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Preface 



In response to the growing interest in building a Neutrino Factory to pro- 
duce high intensity beams of electron- and muon-neutrinos and antineutrinos, 
in October 1999 the Fermilab Directorate initiated two six-month studies. The 
first study, organized by N. Holtkamp and D. Finley, was to investigate the 
technical feasibility of an intense neutrino source based on a muon storage ring. 
This design study has produced a report in which the basic conclusion is that a 
Neutrino Factory is technically feasible, although it requires an aggressive R&D 
program. The second study, which is the subject of this report, was to explore 
the physics potential of a Neutrino Factory as a function of the muon beam 
energy and intensity, and for oscillation physics, the potential as a function of 
baseline. 

The work presented in this report is the result of the enthusiastic contri- 
butions of many people from many institutions. This enthusiasm made the 
organizers job fun. We also want to thank our local sub-group organizers and 
sub-editors for their many effective contributions, ranging from running the 
study groups to editing the report: Bob Bernstein, Debbie Harris, Eric Hawker, 
Stephen Parke, Panagiotis Spentzouris, and Chris Quigg. 

Neutrino Factories seem to have caught the imagination of the community. 
We hope that this report goes some way towards documenting why. 



Steve Geer and Heidi Schellman 
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Executive Summary 



In the Fall of 1999, the Fermilab Directorate chartered a study group to 
investigate the physics motivation for a neutrino factory based on a muon storage 
ring that would operate in the era beyond the current set of neutrino-oscillation 
experiments. Wc were charged to evaluate the prospective physics program as a 
function of the stored muon energy (up to 50 GeV), the number of useful muon 
decays per year (in the range from 10^^ to 10^^ decays per year), and the distance 
from neutrino source to detector. A companion study evaluated the technical 
feasibility of a neutrino factory and identified an R&D program that would lead 
to a detailed design. Our conclusion is that there is a compelling physics case for 
a neutrino factory with a beam energy of about 20 GeV or greater, that initially 
provides at least 0(10^^) muon decays per year. 

The principal motivation for a neutrino factory is to provide the intense, con- 
trolled, high-energy beams that will make possible incisive experiments to pursue 
the mounting evidence for neutrino oscillations. The composition and spectra 
of intense neutrino beams from a muon storage ring will be determined by the 
charge, momentum, and polarization of the stored muons, through the decays 
H~ e'v^Ve or /x"*" e^v^Ve- There is no other comparable source of electron 
neutrinos and antineutrinos. In addition, a neutrino factory would provide well 
coUimated muon neutrino and antineutrino beams. The uncertainties on the 
beam composition and flux are expected to be significantly better than those 
for conventional neutrino beams. If the neutrino factory energy exceeds about 
20 GeV the neutrino beam intensity greatly exceeds the corresponding intensity 
provided by conventional wide band beams. The neutrino factory therefore of- 
fers unprecedented opportunities for precise measurements of nucleon structure 
and of electroweak parameters. The intense muon source needed for the neutrino 
factory would make possible exquisitely sensitive searches for muon-electron con- 
version and other rare processes. 

Experiments carried out at a neutrino factory within the next decade can 
add crucial new information to our understanding of neutrino oscillations. By 
studying the oscillations of i/^, i/g, i?^, and i?e, it will be possible to measure, 
or put stringent limits on, all of the appearance modes ~^ ^t, ~^ J^n, 
and z/^ Ur- This is a necessary step beyond the measurements provided by 
the next generation of neutrino experiments, and will provide a basic test of 
our understanding of neutrino oscillations. It will also be possible to determine 
precisely (or place stringent limits on) all of the leading oscillation parameters, 
including the mixing angle ^13 which appears to be difficult to determine precisely 
with conventional neutrino beams. In addition, a neutrino factory would enable 
us to infer the pattern of neutrino masses; and, under the right circumstances, to 
observe CP violation in the lepton sector. Baselines greater than about 2000 km 
will enable a quantitative study of matter effects and a determination of the mass 
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Wrong-Sign Muon Measurements 
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Figure I: Predicted ratios of Ve i^fj. to z/g —>■ rates at a 20 GeV neutrino 
factory. The upper (lower) band is for 5m\2 < {6771^2 > 0)- The range of 
possible CP violation determines the widths of the bands. The statistical error 
shown corresponds to 10^° muon decays of each sign and a 50 kt detector. Results 
are from Ref. 51. 

hierarchy. If the MiniBooNE experiment confirms the ^ z/g effect reported 
by the LSND experiment, experiments with rather short baselines (a few tens 
of km) could be extremely rewarding, and enable, for example, the search for 
z/g Ur oscillations. 

If the atmospheric neutrino deficit is correctly described by three flavor oscil- 
lations with in the range favored by the SuperKamionkande data, and if the 
parameter sin^ 2^13 is not smaller than ~ 0.01, then exciting cutting-edge long 
baseline oscillation physics could begin with an ~ 50 kt detector at a neutrino 
factory with muon energies as low as 20 GeV delivering as few as 10^^ muon 
decays per year. This "entry-level" facility would be able to measure Ue —>■ 
and Vf. — > z7^ oscillations. For baselines of a few thousand km the ratio of rates 
N{ye — ^ Vp)lN{yf, u^) is sensitive to the sign of 6m'^, and hence to the pat- 
tern of neutrino masses (Fig. I). With 10^^ decays and a 50 kt detector a unique 
and statistically significant measurement of the neutrino mass spectrum could 
be made. In addition, the z/g — >■ t'^ event rate is approximately proportional to 
the parameter sin^2^i3, which could therefore be measured. 

With higher beam intensities and/or higher beam energies the physics po- 
tential of a neutrino factory is enhanced (Fig. II). In particular, as the intensity 
is increased to O(10^°) decays/year z/g —>■ oscillations might be measured, and 
eventually CP violation in the lepton sector observed if the large mixing angle 
MSW solution is the correct description of the solar neutrino deficit. Higher 
beam intensities would also allow smaller values of sin^ 26*13 to be probed (Fig. 
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Figure II: The required number of muon decays needed in a neutrino factory to 
observe —>■ oscillations in a 50 kt detector and determine the sign of (5m^, 
and the number of decays needed to observe oscillations in a few kt 

detector, and ultimately put stringent limits on (or observe) CP violation in the 
lepton sector with a 50 kt detector. Results are from Ref. 51. 

Ill), and higher precision measurements of the oscillation parameters to be made. 
An example of the improvement of measurement precision with neutrino factory 
intensity is shown in Fig. IV for the determinations of sin^ ^23 and sin^ 2^13 . 

The physics program at detectors located close to the neutrino factory is 
also compelling. The neutrino fluxes are four orders of magnitude higher than 
those from existing beams. Such intense beams make experiments with high 
precision detectors and low mass targets feasible for the first time. Using these 
detectors and the unique ability of neutrinos to probe particular flavors of quarks 
will allow a precise measurement of the individual light quark contents of the 
nucleon in both an isolated and nuclear environment. In addition, neutrinos 
provide an elegant tool for probing the spin structure of the nucleon and may 
finally enable resolution of the nucleon spin among its partonic components. 
The high event rates at a neutrino factory would also enable a new generation of 
tagged heavy quark production experiments, precision measurements of electro- 
weak and strong interaction parameters, and searches for exotic phenomena other 
than oscillations. 

Recommendat ions 

The physics program we have explored for a neutrino factory is compelling. We 
recommend a sustained effort to study both the physics opportunities and the 
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Figure III: Limits on sin^ 2^13 that would result from the absence of a z/g — ^ 
signal in a 10 kt detector 7400 km downstream of a 30 GeV neutrino factory 
in which there are 10^° and 10^^/i^ decays, followed by the same number of fi~ 
decays. The limits are shown as a function of The impact of including 

backgrounds in the analysis is shown. Note that the unshaded band shows the 
Sm^ region favored by the SuperK atmospheric neutrino deficit results. Results 
are from Ref. 50. 
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Figure IV: Precision with which the oscillation parameters sin^ ^23 and sin^ 2^13 
can be measured in a 10 kt detector 7400 km downstream of a 30 GeV neutrino 
factory in which there are 10^^, 10^°, and lO^^yU"*" decays. Results are from Ref. 50. 
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machine realities. 

(i) We encourage support for the R&D needed to learn whether a neutrino 

factory can be a real option in the next decade. 

(ii) We propose further studies of detector technologies optimized for a neutrino 

factory, including both novel low mass detectors for near experiments and 
very high mass detectors for long baselines. For long baseline experiments 
detectors should have masses of a few times 10 kt or more that are able 
to detect and measure wrong-sign muons, and detectors of a few kt or 
more able to observe tau-lepton appearance with high efficiency. It is also 
desirable to identify electrons, and if possible measure the sign of their 
charge. Both the detector technologies themselves and the civil engineering 
issues associated with the construction of such massive detectors need to 
be addressed. 

(iii) We recommend continued studies to better compare the physics potential 

of upgraded conventional neutrino beams with the corresponding potential 
at a neutrino factory, and also studies to better understand the benefits of 
muon polarization. 

(iv) The present study concentrated on the muon storage ring as a neutrino 
source and did not cover the additional physics programs which would use 
the proton driver and the high intensity muon beams. We recommend a 
further study directed at these other facets of physics at a muon storage 
ring facility. 
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1 Introduction 



New accelerator technologies offer the possibility of building, not too many years 
in the future, an accelerator complex to accumulate more than 10^^, and per- 
haps more than 10^°, muons per year [1]. It has been proposed [2] to build a 
Neutrino Factory by accelerating the muons from this intense source to energies 
of several GeV or more, injecting the muons into a storage ring having long 
straight sections, and exploiting the intense neutrino beams that are produced 
by muons decaying in the straight sections. If the challenge of producing, cap- 
turing, accelerating, and storing a millimole of unstable muons can be met, the 
decays 

At" e'u^Ue , e^i'^Ue (1) 

offer exciting possibilities for the study of neutrino interactions and neutrino 
properties [2, 3, 4, 5]. In a Neutrino Factory the composition and spectra of 
intense neutrino beams will be determined by the charge, momentum, and po- 
larization of the stored muons. The prospect of intense, controlled, high-energy 
beams of electron neutrinos and antineutrinos — for which we have no other plau- 
sible source — is very intriguing. 

Neutrinos — weakly interacting, nearly massless elementary fermions — have 
long been objects of fascination, as well as reliable probes. One of the most 
dramatic recent developments in particle physics is the growing evidence that 
neutrinos may oscillate from one species to another during propagation, which 
implies that neutrinos have mass. 

If neutrinos z^i, z/2, . . . have different masses mi, m2, . . . , each neutrino flavor 
state may be a mixture of different mass states. Let us consider two species for 
simplicity, and take 

( \ ^ ( cos 9 sm9 \ f ui \ , 
[j^^^ J ~ [ - sine cose )[jy2 J ' 

The probability for a neutrino born as i/^ to oscillate into a i/g, 

P{v. Ue) = sin^ 2esm^ 1.27 ^ • — , 3 

^ " ^' \ 1 eV^ 1 km y ' ^ ' 

depends on two parameters related to experimental conditions: L, the distance 
from the neutrino source to the detector, and E, the neutrino energy. It also de- 
pends on two fundamental neutrino parameters: the difference of masses squared, 
Srn^ = m\ — rn^, and the neutrino mixing parameter, sin^2^. The probability 
that a neutrino born as remain a at distance L is 

PW, ^ = 1 - sm=2»sm= ^1.27 ^ . _ . _j . (4) 

Many experiments have now used natural sources of neutrinos, neutrino radi- 
ation from fission reactors, and neutrino beams generated in particle accelerators 
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to look for evidence of neutrino oscillation. The positive indications for neutrino 
oscillations fall into three classes: [6] 

1. Five solar-neutrino experiments report deficits with respect to the predic- 
tions of the standard solar model: Kamiokande and Supcr-Kamiokande 
(SuperK) using water- Cerenkov techniques, SAGE and GALLEX using 
chemical recovery of germanium produced in neutrino interactions with 
gallium, and Homestake using radiochemical separation of argon produced 
in neutrino interactions with chlorine. These results suggest the oscillation 

^x, with |5m^|soiar ~ 10~^ cV^ and sin^ 2^soiar ~ 1 or a few x 10~^, or 
1 I solar ~ 10"^° eV^ and sin^ 26'soiar ~ 1- 

2. Five atmospheric-neutrino experiments report anomalies in the arrival of 
muon neutrinos: Kamiokande, 1MB, and SuperK using water- Cerenkov 
techniques, and Soudan 2 and MACRO using sampling calorimetry. The 
most striking result is the zenith-angle dependence of the i/^ rate re- 
ported last year by SuperK [7, 8]. These results suggest the oscillation 
^fj. or Vg, with sin^2^atm ~ 1 and |5m^|atm = 10"^ to 10~^ eV^. The 
oscillation i/^ — > i/^ is increasingly the favored interpretation. 

3. The LSND experiment [9] reports the observation of Pg-like events in what 
should be an essentially pure beam produced at the Los Alamos Meson 
Physics Facility, suggesting the oscillation P^ — Pg. This result has not yet 
been reproduced by any other experiment. The favored region lies along 
a band from (sin^ 2^lsnd = 10"^, |(5m^|LSND ~ 1 eV^) to (sin^ 26'lsnd ~ 
1, \Sm\sND ~ 7 X 10-2 eV^). 

A host of other experiments have failed to turn up evidence for neutrino oscilla- 
tions in the regimes of their sensitivity. These results limit neutrino mass-squared 
differences and mixing angles. In more than a few cases, positive and negative 
claims are in conflict, or at least face off against each other. Over the next five 
years, many experiments will seek to verify, further quantify, and extend these 
claims. If all of the current experimental indications of neutrino oscillation sur- 
vive, there are apparently three different mass-squared-difference scales, which 
cannot be accommodated with only three neutrino types. New sterile neutrinos 
may be required. This would be a profound discovery. 

From the era of the celebrated two-neutrino experiment [10] to modern times, 
high-energy neutrino beams have played a decisive role in the development of 
our understanding of the constituents of matter and the fundamental interac- 
tions among them. Major landmarks include the discovery of weak neutral- 
current interactions [11], and incisive studies of the structure of the proton and 
the quantitative verification of perturbative quantum chromo dynamics as the 
theory of the strong interactions [12]. The determinations of the weak mix- 
ing parameter sin^ 9w and the strong coupling constant in deeply inelastic 
neutrino interactions are comparable in precision to the best current measure- 
ments. Though experiments with neutrino beams have a long history, beams 
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of greatly enhanced intensity would bring opportunities for dramatic improve- 
ments. Because weak-interaction cross sections are small, high-statistics studies 
have required massive targets and coarse-grained detectors. Until now, it has 
been impractical to consider precision neutrino experiments using short liquid 
hydrogen targets, or polarized targets, or active semiconductor target-detectors. 
All of these options are opened by a muon storage ring, which would produce 
neutrinos at approximately 10"^ times the flux of existing neutrino beams. 

At the energies best suited for the study of neutrino oscillations — tens of GeV, 
by our current estimates — the muon storage ring is compact. We could build it at 
one laboratory, pitched at a deep angle, to illuminate a laboratory on the other 
side of the globe with a neutrino beam whose properties we can control with 
great precision. By choosing the right combination of energy and destination, 
we can tunc future neutrino-oscillation experiments to the physics questions we 
will need to answer, by specifying the ratio of path length to neutrino energy 
and determining the amount of matter the neutrinos traverse. Although we 
can use each muon decay only once, and we will not be able to select many 
destinations, we may be able to illuminate two or three well-chosen sites from 
a muon-storage-ring neutrino source. That possibility — added to the ability to 
vary the muon charge, polarization, and energy — may give us just the degree 
of experimental control it will take to resolve the outstanding questions about 
neutrino oscillations. Experiments at a Neutrino Factory would seek to verify the 
number of neutrino types participating in the oscillations, precisely determine 
the mixing parameters that relate the flavor states to the mass states, determine 
the pattern of neutrino masses, and look for CP violation in the lepton sector. 

The prodigious flux of neutrinos close to the muon storage ring raises the 
prospect of neutrino-scattering experiments of unprecedented sensitivity and del- 
icacy. Experiments that might be pursued at a Neutrino Factory include precise 
measurements of the nucleon structure (including changes that occur in a nu- 
clear environment), measurements of the spin structure of the nucleon using a 
new and powerful technique, charm measurements with several million tagged 
particles, precise measurements of Standard Model parameters, and searches for 
exotic phenomena. 

We beheve that the physics program at a Neutrino Factory is compeUing and 
encourage support for a vigorous R&D program to make neutrino factories a real 
option for the future. 
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2 Beam properties 



Consider an ensemble of polarized negatively-charged muons. When the muons 
decay they produce muon neutrinos with a distribution of energies and angles in 
the muon rest-frame described by [13]: 

-—^ oc [(3-2a;) + (l-2a;)P^coseem] , (5) 
axa\ Icm 47r 

where x = 2Et,/'m^, 6cm is the angle between the neutrino momentum vector 
and the muon spin direction, and is the average muon polarization along the 
beam direction. The electron antineutrino distribution is given by: 

/72/V- 12r^ 

-f^ oc t^[{l-x) + {l-x)P,cos9cm\. (6) 

dxdi Icm 47r 

and the corresponding distributions for and z/g from /i"*" decay are obtained 
by the replacement "-^m- Oiily neutrinos and antineutrinos emitted in the 

forward direction (cos6';afe ~ 1) are relevant to the neutrino flux for long-baseline 
experiments; in this limit Ey — xE^ax and at high energies the maximum Ey 
in the laboratory frame is given by Emax = 7(1 + /9cos ^cm)?7i^/2, where f3 and 
7 are the usual relativistic factors. The and Vc distributions as a function of 
the laboratory frame variables are then given by: 



d'^N 1 2r2 

Tlur- ^ ^in—R ^i-[(3-2a;) + (l-2a;)P,cos^e,n], 

dxdiliab 7^(1 - /? COS eiabY 47r 



(7) 



and 



d'^N- 1 12r^ 



dxdQiab 7^(1 - /3 cos 6l/„b)^ 47r 

Thus, for a high energy muon beam with no beam divergence, the neutrino 
and antineutrino energy- and angular- distributions depend upon the parent 
muon energy, the decay angle, and the direction of the muon spin vector. With 
the muon beam intensities that could be provided by a muon-coUider type muon 
source [1] the resulting neutrino fluxes at a distant site would be large. For 
example. Fig. 1 shows as a function of muon energy and polarization, the com- 
puted fluxes per 2 x 10^° muon decays at a site on the other side of the Earth 
(L = 10000 km). Note that the (l7g) fluxes are suppressed when the muons 
have P = -1-1 (-1). This can be understood by examining Eq. (8) and noting 
that for P = — 1 the two terms cancel in the forward direction for all x. 



2.1 Interaction rates 

Neutrino charged current (CC) scattering cross-sections are shown as a function 
of energy in Fig. 2. At low energies the neutrino scattering cross section is 
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Neutrino Fluxes at L = 10000 Km 
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Figure 1: Calculated v and v fluxes in the absence of oscillations at a far site 
located 10000 km from a neutrino factory in which 2 x 10^° muons have decayed 
in the beam-forming straight section. The fluxes are shown as a function of 
the energy of the stored muons for negative muons (top two plots) and positive 
muons (bottom two plots), and for three muon polarizations as indicated. The 
calculated fluxes are averaged over a circular area of radius 1 km at the far site. 
Calculation from Ref. 2. 
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Figure 2: The total cross section for charged current neutrino scattering by muon 
and tau neutrinos (top plot), and the ratio of tau to muon neutrino cross sections 
function of neutrino energy (bottom plot). 
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Figure 3: Charged current event spectra at a far detector. The sohd hues indicate 
zero polarization, the dotted hues indicate polarization of ±0.3 and the dashed 
lines indicate full polarization. The P = 1 case for electron neutrinos results in 
no events and is hidden by the x axis. 



dominated by quasi-elastic scattering and resonance production. However, if 
is greater than ~ 10 GeV, the total cross section is dominated by deep inelastic 
scattering and is approximately [14]: 

o{y ^ N ^ t ^ X) ^ 0.67 x 10"^^ cm^ x {E^, GeV) , (9) 
a{V+N -^e+ + X) ^ 0.34 x 10"^^ cm^ x (E^, GeV) . (10) 

The number of u and 17 CC events per incident neutrino observed in an isoscalar 
target is given by: 

N{u + N +X) = 4.0 X lO'^^iE^, GeV) events per gr/cm^ , (11) 
N(T7 + N t + X) = 2.0 X 10'^\E^, GeV) events per gr/cm^ . (12) 

Using this simple form for the energy dependence of the cross section, the pre- 
dicted energy distributions for z/g and interacting in a far detector (cos6' = 1) 
at a neutrino factory are shown in Fig. 3. The interacting energy distribu- 
tion is compared in Fig. 4 with the corresponding distribution arising from the 
high-energy NUMI wide band beam. Note that neutrino beams from a neutrino 
factory can be considered narrow band beams. In practice, CC interactions can 
only be cleanly identified when the final state lepton exceeds a threshold energy. 
The calculated final state lepton distributions are shown in Fig. 5. Integrating 
over the energy distribution, the total u and u interaction rates per muon decay 
are given by: 

. (L, kmy 



= 1.2 X 10 



-14 



xC(z/) events per kt (13) 
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Figure 4: Comparison of interacting z/^ energy distributions for the NUMI high 
energy wide band beam (Ref. 15) with a 20 GeV neutrino factory beam (Ref. 2) at 
L = 730 km and a 30 GeV neutrino factory beam at L = 2900 km. The neutrino 
factory distributions have been calculated based on Eq. (5) (no approximations), 
and include realistic muon beam divergences and energy spreads. 




Figure 5: Lepton energy spectra for CC (top left), (top right), (bottom 
left), and Ve (bottom right) interactions. Note that z is the energy normalized 
to the primary muon energy z = E^/ E^. Calculation from Ref. 16. 
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Table 1: Muon neutrino and electron antineutrino CC interaction rates in the 
absence of oscillations, calculated for baseline length L — 732 km (FNAL — > 
Soudan), for MINOS using the wide band beam and a muon storage ring deliv- 
ering 10^° decays with En = 10, 20, and 50 GeV at 3 baselines. The neutrino 
factory calculation includes a realistic muon beam divergence and energy spread. 



Experiment 




Baseline 
(km) 


(GeV) 


{E.J 
(GeV) 


N(z/^ CC) 
(per kt-yr) 


N(z7e CC) 
(per kt-yr) 


MINOS 


Low energy 


732 


3 




458 


1.3 




Medium energy 


732 


6 




1439 


0.9 




High energy 


732 


12 




3207 


0.9 


Muon ring 


(GeV) 










10 


732 


7.5 


6.6 


1400 


620 




20 


732 


15 


13 


12000 


5000 




50 


732 


38 


33 


1.8x10^ 


7.7x10^ 


Muon ring 


E^, (GeV) 










10 


2900 


7.6 


6.5 


91 


41 




20 


2900 


15 


13 


740 


330 




50 


2900 


38 


33 


11000 


4900 


Muon ring 


(GeV) 










10 


7300 


7.5 


6.4 


14 


6 




20 


7300 


15 


13 


110 


51 




50 


7300 


38 


33 


1900 


770 



and 



Njj = 0.6 X 10 



-14 



where 



{E„GeV) 



3i 



(L, kmy 



xC{u) events per kt , 



C(iy.) = — + P. — , C(i/e) = — -Pu— 



(14) 



(15) 
(16) 



The calculated Ue and CC interaction rates resulting from 10^° muon 
decays in the beam-forming straight-section of a neutrino factory are compared 
in Table 1 with expectations for the corresponding rates at the next generation 
of accelerator-based neutrino experiments. Note that event rates at a neutrino 
factory increase as E^, and are significantly larger than expected for the next 
generation of approved experiments if E^ > 20 GeV. The radial dependence 
of the event rate is shown in Fig. 6 for a 20 GeV neutrino factory and three 
baselines. 



18 




20 40 60 80 100 

Radius in km 



Figure 6: Events/kT of detector as a function of distance from the beam center 
for a 20 GeV muon beam. 

Finally, for an isoscalar target the neutral current (NC) cross sections are 
approximately 0.4 of the CC cross sections[17], and are given by: 

a{v + N^v + X) ^ 0.3 X 10-3^ cm^ X (E^, GeV) , (17) 
a(V + N^V + X) ^ 0.15 X 10"^^ cm^ X (E^,GeV) . (18) 

2.2 Tau neutrino interactions 

Tau neutrino CC interaction rates are substantially less than the corresponding 
i^e and Ufj^ rates, especially near the tau production threshold of ~ 3.3 GeV. The 
NC rates should be the same as those for electron and muon neutrinos. Figure 2 
shows the calculated [18] ratio of i^t/i^^ CC interaction rates as a function of 
the neutrino energy. Near threshold, contributions from quasi-elastic and reso- 
nance production dominate. If the cross sections from Ref. [19] are used, the 
predicted event rates are 5-7% lower. 

2.3 Systematic uncertainties on the muon beam and neu- 
trino flux 

In the neutrino beam-forming straight section the muon beam is expected to 
have an average divergence given by ae = 0(0.1/7). The neutrino beam diver- 
gence will therefore be dominated by muon decay kinematics, and uncertainties 
on the beam direction and divergence will yield only small uncertainties in the 
neutrino flux at a far site. However, if precise knowledge of the flux is required, 
the uncertainties on 6 and ae must be taken into account, along with uncertain- 
ties on the flux arising from uncertainties on the muon energy distribution and 
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Tabic 2: Dependence of predicted charged current event rates on muon beam 
properties at a neutrino factory. The last column lists the required precisions 
with which each beam property must be determined if the uncertainty on the 
neutrino flux at the far site is to be less than ~ 1%. Here A denotes uncertainty 
while a denotes the spread in a variable. Table from Ref. 20. 



Muon Beam 
property 


Beam 
Type 


Rate Target 
Dependence Precision 


Energy {E^,) 


V (no osc) 

7/,, Uf, 


AN/N = 3 AE^jEf,, A{E^) /E^, < 0.003 
AN/N = 2 AE,/E, ME,,)/E, < 0.005 


Direction {M) 


u (no osc) 


AN/N < 0.01 < 0.6 ae 
(for A^ < 0.6 ae) 


Divergence (erg) 


u (no osc) 


AN/N ~ 0.03 Aae/ag Aae/ag < 0.2 
(for ae ~ O.I/7) (for ae ~ O.I/7) 


Momentum spread (cTp) 


u (no osc) 


AN/N ~ 0.06 Aap/ap Aa^/a^ < 0.17 


Polarization (P^) 


z/e (no osc) 
1/^ (no osc) 


ANJN,^ = AP^ AP^ < 0.01 
ANJN,^ = 0.4 AP^ AP^ < 0.025 



polarization. The relationships between the uncertainties on the muon beam 
properties and the resulting uncertainties on the neutrino flux are summarized 
in Table 2. If, for example, we wanted to know the and fluxes at a far site 
with a precision of 1%, we would need to know the beam divergence uq to 20% 
(Fig. 7), and ensure that the beam direction was within 0.6 ug of the nominal 
direction [20] (Fig. 8). 

2.4 Event distributions at a near site 

The event distributions measured in a detector close to the neutrino factory will 
be quite different from the corresponding distributions at a far site. There are 
two main reasons for this difference. First, the near detector accepts neutrinos 
over a large range of muon decay angles 9, not just those neutrinos traveling 
in the extreme forward direction. This results in a broader neutrino energy 
distribution that is sensitive to the radial size of the detector (Fig. 9). Second, if 
the distance of the detector from the end of the beam forming straight section is 
of the order of the straight section length, then the 6 acceptance of the detector 
varies with the position of the muon decay along the straight section. This 
results in a more complicated radial flux distribution than expected for a far 
detector. However, since the dominant effects are decay length and muon decay 
kinematics, it should be modeled quite accurately. (Fig. 10). 

Note that, even in a limited angular range, the event rates in a near detector 
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Beam divergence (mrads) 

Figure 7: Dependence of CC interaction rates on the muon beam divergence 
for a detector located at L = 2800 km from a muon storage ring containing 
30 GeV unpolarized muons. Rates are shown for z/g (boxes) and (circles) 
beams in the absence of oscillations, and for oscillations (triangles) with 

the three-flavor oscillation parameters lAl. The calculation is from Ref. 20. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Beam direction (mrads) 



Figure 8: Dependence of CC interaction rates on the neutrino beam direction. 
Relative rates are shown for a detector at a far site located downstream of a 
storage ring containing 30 GeV unpolarized muons, and a muon beam divergence 
of 0.33 mr. Rates are shown for z/g (triangles) and (circles) beams in the 
absence of oscillations, and for z/g — * oscillations (boxes) with the three-flavor 
oscillation parameters lAl. The calculation is from Ref. 20. 
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Figure 9: Events per gr/cm^ per GeV for a detector 40 m from a muon storage 
ring with a 600 m straight section. The 3 curves show all events and those falling 
within 50 and 20 cm of the beam center. 



Detector 40 m from 600 m decay region 


= 50 GeV 




1 0^° [L decays I 




target 1 gr/cm^ j 









-100 -50 50 100 

X, cm 



Figure 10: Events per gr/cm^ as a function of the transverse coordinate x 50 m 
downstream of a 50 GeV neutrino factory providing 10^° muon decays. The 
central peak is mainly due to decays in the last hundred meters of the decay 
pipe while the large tails are due to upstream decays. 
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Figure 11: Events per gr/cm^ at a near detector as a function of muon beam 
energy. The curves indicate (solid) all events, the dashed and dotted curves show 
the effects of radial position cuts. 

are very high. Figure 11 illustrates the event rates per gram/cm^ as a function 
of energy. Because most of the neutrinos produced forward in the center of mass 
traverse the detector fiducial volume, the factor of 7^ present in the flux for 6* ~ 
is lost and the event rate increases linearly with E^^. For a 50 GeV muon storage 
ring, the interaction rate per 10^° muon decays is 7 million events/gram/cm^. 
Rate calculations are discussed further in the context of specific experiments in 
the section on non-oscillation experiments. Finally, in the absence of special 
magnetized shielding, the high neutrino event rates in any material upstream 
of the detector will cause substantial backgrounds. The event rate in the last 
3 interaction lengths (300 gr/cm^) of the shielding between the detector and 
the storage ring would be 30 interactions per beam spill at a 15 Hz machine 
delivering 2 x 10^° muon decays per year. These high background rates will 
require clever magnetized shielding designs and fast detector readout to avoid 
overly high accidental rates in low mass experiments. 
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3 Oscillation physics 



The recent impressive atmospheric neutrino results from the SuperK experiment 
have gone a long way towards establishing the existence of neutrino oscilla- 
tions [21]. Up to the present era, neutrino oscillation experiments at accelerators 
were searches for a phenomenon that might or might not be within experimental 
reach. The situation now is quite different. The atmospheric neutrino deficit 
defines for us the and oscillation amplitude that future long-baseline os- 
cillation experiments must be sensitive to, namely = O(10~^) eV^ and 
sin^2^ = 0(1)- Experiments that achieve these sensitivities are guaranteed an 
excellent physics program that addresses fundamental physics questions. We can 
hope that future neutrino oscillation experiments will provide the keys we need to 
understand really fundamental questions, for example: the origin of the minute 
neutrino masses and the reason why there are three lepton families. We can- 
not guarantee that these insights will be forthcoming from neutrino oscillation 
measurements, but they might be. For this reason it is important to understand 
how our community can get detailed experimental information on the neutrino 
oscillation scheme, the mass splittings between the neutrino mass eigenstates, 
and the leptonic mixing matrix that controls the oscillation probabilities. A 
neutrino factory would be a new tool, providing a beam of energetic electron 
neutrinos. In the following we address how this new tool might be exploited 
to go well beyond the capabilities of the next generation of neutrino oscillation 
experiments. 

In this section we begin by describing the theoretical basis for neutrino oscil- 
lations, and then define a selection of oscillation parameter sets that can be used 
in assessing the physics program at a neutrino factory. This is followed by a sum- 
mary of the current experimental status and how it can be expected to change 
in the next few years. We then discuss the parameters and the performance 
of candidate detectors at a neutrino factory. The section is completed with a 
survey of the physics measurements that can be performed at a neutrino factory 
as a function of beam energy, intensity, and baseline, and finally, a summary of 
our conclusions. 

3.1 Theoretical framework 

There exist three known flavors of active neutrinos which form left-handed dou- 
blets with their associated charged leptons. The interaction of these active neu- 
trinos with the electroweak gauge bosons is described by the Standard Model 
(SM). In principle there can be additional flavors of neutrino which are sin- 
glets under the electroweak gauge group. These electroweak singlet neutrinos 
do not have electroweak couplings, and their interactions arc not described by 
the SM. Let us denote the flavor vector of the SU(2) x U(l) active neutrinos as 
V — {ve, I'fjL-i i^t) and the vector of electroweak-singlet neutrinos as x = (xi, .., Xn^)- 
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The Dirac and Major ana neutrino mass terms can then be written compactly as 

where Ml is the 3x3 left-handed Majorana mass matrix, Mr is a x right- 
handed Majorana mass matrix, and Md is the 3-row by n^-column Dirac mass 
matrix. In general, all of these are complex, and (Ml)^ = Ml , (Mr)^ = Mr. 
Without further theoretical input, the number Ug of "sterile" electroweak-singlet 
neutrinos is not determined. For example, in the SM, minimal supersymmetric 
standard model (MSSM), or minimal SU(5) grand unified theory (GUT), Us = 0, 
while in the SO(IO) GUT, Ug — 3. (This is true for both the original non- 
super symmetric and the current supersymmetric versions of these GUTs.) Since 
the terms XjRCxkR ai'e electroweak singlets, the elements of the matrix Mr, 
would not be expected to be related to the electroweak symmetry breaking scale, 
but instead, would be expected to be much larger, plausibly of the order of the 
GUT scale. 

Mechanisms involving Ml only for the generation of neutrino masses without 
the presence of electroweak-singlet neutrinos exist. The simplest scenarios, in 
which one or more Higgs triplets are introduced to couple to a pair of left-handed 
neutrinos, arc excluded by measurements of the p parameter. Therefore, other 
extensions of the SM must be considered, for example the addition of one or 
more Higgs singlets, non-renormalizable terms involving a large mass scale such 
as the GUT scale, or R-parity- violating terms in the context of supersymmetry. 

The most natural explanation for the three known ultra-light neutrino masses 

is generally regarded to be the seesaw mechanism [22], which involves Mr, and 
arises from Eq. (19) in the case of = 3 electroweak singlet neutrinos. This 
leads to neutrino masses generically of order 



rriR 

where mD and rriR denote typical elements of the corresponding matrices. With 
rriL, ~ rrit and rriR ~ 10^^ GcV, as suggested in a (supersymmetric) SO(IO) 
grand unified theory framework, a scale of ~ 10^'^ eV is readily obtained. In 
this case the three light neutrino masses are obtained by diagonalization of the 
effective 3x3 light neutrino mass matrix 

M^ = -MdM^^MI (21) 

while the super-heavy neutrinos are determined from the right-handed Majorana 
matrix Mr. 

Additional electroweak-singlet neutrinos may arise in string theory with the 
existence of supersymmetric partners of moduli fields, resulting in the appearance 
of ni light sterile neutrinos. But the presence of these light sterile neutrinos may 
undermine the seesaw mechanism and, for this reason, is not very appealing. 
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(20) 



However, if one tries to fit all of the data from the oscillation experiments, to 
obtain a reasonable it is necessary to include light sterile neutrinos. We shall 
illustrate some of the effects of sterile neutrinos with a toy model in which one 
studies the minimal number, = 1. 



3.1.1 Neutrino Oscillations in Vacuum 



The presence of non-zero masses for the light neutrinos introduces a leptonic 
mixing matrix, U, which is the analogue of the CKM quark mixing matrix, and 
which in general is not expected to be diagonal. The matrix U connects the 
flavor eigenstates with the mass eigenstates: 



K)^^Uai\ui), (22) 

i 

where a denotes one of the active neutrino flavors, e, fioi r or one of the rii light 
sterile flavors, while i runs over the light mass eigenstate labels. The number of 
flavor states considered here is equal to the number of light mass eigenstates, so 
[/ is a square unitary matrix. 

The neutrino mass differences and the mixing parameters can be probed by 
studying oscillations between different flavors of neutrinos, as a function of the 
neutrino energy E and the distance traversed L. The oscillation probability 
P(i/Q — > z/^) is given by the absolute square of the overlap of the observed flavor 
state, 1 1^^), with the time-evolved initially- produced flavor state, \i'o)- In vacuum, 
the evolution operator involves just the Hamiltonian Hq of a free particle, yielding 
the well-known result: 



The CP-even and CP-odd contributions are 

= Sap-^Ei>j Re {UaiU^iU*jUf)j) sin^i 

-PcP-odd(j^a ^p) — —PcP-oddi^a ^ ^p) 

Am^ T 



(23) 



4B ^ 



(24) 



so that 

P{Pa T>p) = P{vp Ua) = PcP-even(l^a ^ i^p) " -PcP-odd(l^a ^ l^p) (25) 

where, by CPT invariance, Pli'a ^p) = P{i^i3 ^a)- In vacuum the CP-even 
and CP-odd contributions are even and odd, respectively, under time reversal: 
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In Eq. (24), Smfj = m,{v,-i)'^ — m^Uj)^, and the combination 6mfjL/{4E) 
in h = c = 1 units can be replaced by 1.2669 • • ■dmjjL/E with dm^j in eV^ and 
(L, E) in (km, GeV). In disappearance experiments (5 — a and no CP- violation 
can appear since the product of the mixing matrix elements is inherently real. 
At distances L large compared to all the individual oscillation lengths, A"!"^ ~ 
E/Smfj, the sine squared terms in Fcp-cvcn average to 0.5 whereas the sine 
terms in Pcp-odd average to zero. Therefore CP violating effects are largest and 
hence easiest to observe at distances between the smallest and largest oscillation 
lengths. 



3.1.2 Three Active Neutrinos Only 

With three neutrinos, the mixing matrix U is the 3x3 unitary Maki-Nagawa- 
Sakata (MNS) matrix [23]. We parameterize U by 



iS 



U 



( C13C12 C13S12 513^' 

~C23Si2 — Sl3S23Cl2e*'' C23C12 — 8138238126*'' C13S23 | , (26) 
V 'S23S12 — -S13C23C126*'' — S23C12 — -§136235126*^ C13C23 



where Cjk = cosOjk and Sjk = sinO jj.. For Majorana neutrinos, U contains 
two further multiplicative phase factors, but these do not enter in oscillation 
phenomena. 

With the plausible hierarchical neutrino mass spectrum mi < 1712 <^ and 
the assumption that the LSND effect is not a neutrino oscillation phenomena, 
we can identify the largest 5m^ scale with the atmospheric neutrino deficit: 
SM"^ = Sm'l^^ = — ^"^31- the approximation that we neglect oscillations 
driven by the small Srri^ scale, the Ue oscillation probabilities can be written as 

P(^e-^e)^ l-4|C/e3P(l-|C/e3nsin^(%^) 

[27) 



1 -sin2(2^i3)sin2( 



AE 



P(z/e-Z^M)^ 4|C/e3nf/,3psin2(%^ 



Sin2(2^i3)sin2(^23)sin2(%^) , ^^^^ 



P{v,^v;)^ 4|f/.3nf/e3psin2(%^ 



sin2(2^i3)cos2(^23)sin^(%^) ^^^^ 



and the v^j^ oscillation probabilities are 



P{v, ^v,)^ 1 - A\U,,\\l - |f/,3n sin2(^!^^) 

= 1 - 4sin2(^23) cos2(^i3)(l - sin2(^23) cos2(^i3)) sin2(%^) , 

(30) 
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V 4E 



Sin2(2^i3)sin2(^23)sin2(%^) , ^^^^ 



sin2(2^^23) cos4(^^i3) sin^^^^^' ^^^^ 



4_B 



The CP-odd contribution to the atmospheric neutrino oscillation probability 
vanishes in the one-mass-scale-dominant approximation. However if we include 
the effects of the small mass scale, Sm^i-, then 



-PcP-odd(i^M 



sm 



4ci2cf3C23Sl2Sl3S23(sin 5) 
L 



2E 



-) sin 



2 

atm ^ 



4E 



+ sm| 



2E 



'-) sin 



4E ) 



(33) 



At distances significantly larger than the atmospheric neutrino oscillation length, 
E/Sm^^j^, the second term in brackets averages to zero whereas the sin squared 
part of the first term averages to one half, leaving 



-PcP-odd(i^M ^r) ^ -2ci2C?3C23Si2Si3S23(sin(5) Sin(^^^) 



(34) 



The Jarlskog factor [24], J, is given by J = ci2Ci3C23Si2Si3S23(sin(5) and is a 
convenient measure of the size of the CP violation. 

If the neutrinos propagate through matter, these expressions must be modi- 
fied. The propagation of neutrinos through matter is described by the evolution 
equation 



.di'a 

l—r- 

dt 



E 

/3 



rrij 
2E. 



+ 



A 
2E„ 



(35) 



where A/{2Ei,) is the amplitude for coherent forward charged- current scattering 
of i/e on electrons. 



A = 2V2GFNeE^ = 1.52 x 10-\V^Yep{g/cm^)E{ GeV) 



(36) 



(for i7e A is replaced with -A). Here is the electron fraction and p{t) is the 
matter density. Density profiles through the earth can be calculated using the 
Earth Model [25], and are shown in Fig. 12. For neutrino trajectories through the 
earth's crust, the density is typically of order 3 gm/cm^, and Ye ~ 0.5. For very 
long baselines a constant density approximation is not sufficient and oscillation 
calculations must explicitly take account of p(t). However the constant density 
approximation is very useful to understand the physics of neutrinos propagating 
through the earth since the variation of the earth's density is not large. 

The propagation Eq. (35) can be re-expressed in terms of mass-squared dif- 
ferences: 



l—r- 

dt 



E 



2E,. 



(37) 
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Figure 12: Density profiles for trajectories through the Earth. Calculation from 
Ref. 16. 



This evolution equation can be solved numerically for given input values of the 
Sm^ and mixing matrix elements. 

In the approximation where we neglect oscillations driven by the small Sm^ 
scale, the evolution equations are: 



d 
I— 
dt 



f 




Srri^ ( 














Vr ) 





A 



^e3 



\2 



It/. 



t3 



v 



(38) 



For propagation through matter of constant density, the flavor eigenstates are 
related to the mass eigenstates i/™ by 



where 



/ 

-C23 
\ S23 



-13 



'13 



-SL3S23 C^S23 
-SSC23 C^^3C23 / 



(39) 



(40) 



and 9"^ is related to ^13 by 

tan26'^^ = sin 2^13/ (^cos2^i3 - 



(41) 



We note that has the form of the vacuum U with the substitutions 

^13 ^^^l3, ^^23^^^23, ^^i2^vr/2, 5 = 0. (42) 
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Equation (41) implies that 



sin^ 2^^ = sin^ 2^13/ ( " cos 2^13) + sin^ 2^13 j . (43) 

Thus there is a resonant enhancement for 

A = 5m^ cos 2^13 (44) 

or equivalently 

The resonance occurs only for positive 5m^ for neutrinos and only for negative 
Sm? for anti-neutrinos.-^ For negative Sm? the oscillation amplitude in Eq. (43) 
is smaller than the vacuum oscillation amplitude. Thus the matter effects give 
us a way in principle to determine the sign of 5m^. 

It is instructive to look at the dependence of the oscillation probabilities on 
the neutrino energy as a function of the oscillation parameters and the baseline. 
Some examples from Ref. [26] are shown in Fig. 13 for v^. i/^ oscillations. Note 
that for parameters corresponding to the large mixing angle MSW solar solution, 
maximal CP violation results in a small but visible effect. Matter effects, which 
have been computed using the density profile from the Earth Model, can have 
substantial effects, and are very sensitive to sin^ 26'i3. 



3.1.3 Three Active Flavor Oscillation Sceneirios 

We now define some representative three-flavor neutrino oscillation parameter 
sets that can be used to establish how well experiments at a neutrino factory 
could determine the oscillation parameters. We begin by considering constraints 
from existing experiments. 

If we assume CPT invariance then the oscillation probability for Pg ~^ is 
equal to that for v^. v^. The CHOOZ results [27] imply: 

sin^ 2^,eac = 4|[/e3|'(l - It/eaH = siu^ 2^i3 < 0.1 (46) 

for the range SM"^ ^ 10~^ eV^. On the other hand, for the solar neutrino 
experiments, with |C/e3|^ <S 1, one finds 

sin^ 2Q,,i,, = 4 1 [/el n C/e2 r = sin^ 2^12 cos^ ^13 ~ sin^ 2^12 (47) 

-'^If the LSND effect is due to neutrino oscillations then 8m} >> 0(fO~^) eV^ and the 
resonance occurs at energies much higher than those of interest at the currently invisioned 
neutrino factory. 
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Figure 13: Dependence of z/g oscillation probability on neutrino energy 

for some representative oscillation parameters. Plots are from Ref. 26 and show 
the effects of varying 5 (top plots), matter effects (middle plots), and sin^ 2^13 
dependence (bottom plots). 
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with sin^ 2^^i2 ~ 0.006 in the case of the small angle MSW solution with 617121 ~ 
6 X 10~® eV^ or sin^ 2^12 ~ 1.0 in the case of the large angle MSW solution with 
5mli ~ 5 X 10~^ eV^, the LOW solution with 5mli ~ 10~^ eV^, or the vacuum 
solutions with Smli ~ 4 x 10~^° eV^ or Smli ~ 8 x 10~^^ eV^. 

The atmospheric neutrino oscillation experiments favor i/^ — > Vt [28], and in 
the one-mass-scale-dominant approximation the best fit from the SuperK exper- 
iment [29] yields 

sin22^„,^ = 4\U^3\^\Ur3\^ = sin2 2e23COS^^i3 = 1.0 (48) 

with Sm1f„^ = 3.5 x 10~^ eV^. Unpublished analyses of a substantially enlarged 
data set by the SuperK experiment have yielded the same central value for 
sin^ 29atm and essentially the same value of Sni^tj^, 2.5 x 10~^ [29]; we shall use 
the pubhshed fits in the following. 

Based on these considerations we define the representative three-fiavor pa- 
rameter sets shown in Table 3. The first three scenarios do not attempt to fit 
the LSND anomaly. These scenarios have the Atmospheric anomaly explained 
by z/^ — Ur oscillation with maximal mixing and the Solar Anomaly explained 
by one of the MSW Solar solutions: 

Scenario lAl - Large Angle MSW 
Scenario IA2 - SmaU Angle MSW 
Scenario IA3 - LOW MSW. 

Alternatively we can keep the LSND anomaly, and either drop the solar 
neutrino deficit, or attempt to find a "fit" (necessarily with a poor x^) that 
explains all three neutrino anomalies [30]: 

Scenario IBl - Atmospheric and LSND 
Scenario ICl - Atmospheric, Solar and LSND 

For scenario ICl the Atmospheric anomaly is a mixture of z/^ z/^ and z/^ 
z/g and the solar electron neutrino flux is reduced by a factor two independent 
of energy. There are large contributions to the fo^' this scenario coming from 
the Atmospheric Neutrino Anomaly as well as the Homestake (Chlorine) Solar 
neutrino experiment. 

Note that the Jarlskog J-factor is small for all scenarios. It is clear that CP 
violation will be very difficult to observe. 
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Table 3: Parameters for the three- flavor oscillation scenarios defined for the 
study. 



T~\ o T* o yy^ dT QT* 

pdjl djIllcLcl 


TAI 


TA9 


T A '? 


1 R1 


1 PI 


0171^2 16 V ) 
0777.21 J 

sin 2 O^oo 


O.O X iU 


O.O X iU 


O.O X iU 


O.O X iU 


u.o 


O A ±u 


u X ±u 


1 V 1 

1 X lU 


u.o 


7 V 1 n^^ 

( X iU 


1 n 


1 


1 


1 


53 


sin^ 2^13 


0.04 


0.04 


0.04 


0.015 


0.036 


sin^ 2^12 


0.8 


0.006 


0.9 


0.015 


0.89 


S 


0,±7r/2 


0,±7r/2 


0,±7r/2 


0,±7r/2 


0,±7r/2 


sin^ 26'„t^ 


0.98 


0.98 


0.98 


0.99 




Sin^ 20reac 


0.04 


0.04 


0.04 


0.03 




Sin^ 26 solar 


0.78 


0.006 


0.88 






SIV? 2eLSND 








0.03 


0.036 


J 


0.02 


0.002 


0.02 


0.002 


0.015 



3.1.4 Three Active and One Sterile Neutrinos 

In order to incorporate the observed z/^ — > and LSND appearance 

results [9] and achieve an acceptable in the fit, it is necessary to introduce 
at least one light sterile neutrino. As discussed earlier, the theoretical case for 
sterile neutrinos is unclear, and various neutrino mass schemes predict anything 
from Tis = to many. To admit just one must be regarded as a rather unnatural 
choice. We consider this case because it allows us to explain the Atmospheric, 
Solar and LSND anomalies with the fewest number of new parameters. 

Scenarios with three nearly degenerate neutrinos (for example m\ < m2 < 
7/23 ^ 7/74 or mi <^ 7/72 < < ^4) are essentially ruled out by a Schwarz 
inequality on the Icptonic mixing elements [31]: \U^iU*^\'^ < |f/^4p|f/e4p < 0.008 
which fails to be satisfied in the allowed LSND region. Of the two scenarios 

with 777 1 < 7772 <S < ''774, the One with 57772^ ~ ^^lolar^ 5m\.^ ~ ^"^atm 

preferred over the other arrangement which is on the verge of being ruled out by 
the Heidelberg-Moscow /9/9oi/ decay experiment [32] giving (777) < 0.2 eV. 

With the three relevant mass scales given by 

5ml^i = bm\^ < 5777^i„ = 5777^3 < ^rri\s^^ = bml^ 

and the fiavors ordered according to {s, e, /i, r}, the 4x4 neutrino mixing 
matrix depends on six angles and three phases and is conveniently chosen to be 
[33] 
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— Ru{&14, 0)-Rl3(6'l3, 0)i?24(6'24, 0)i?23(^23, ^3)-R34(6'34, ^2)-Rl2(6'l2, ^l) 

where, for example, 



(49) 
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C23 


S23e-'^' 










C23 
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l) 



-^23(^23) ^3) 



In the hmit where the mi — m2 and 1713 — rui pairs are considered degenerate, 
-^12(^12, (^1) = -^34(^34, (^34) = I-i and only four angles and one phase appear in 
the mixing matrix 



/ C14C13 -Ci4Si3S23e"^^ - S14S24C23 C14S13C23 - Sl4S24S23e~"^^ S14C24 \ 

jj ^ C24C23 C24S23e~"^3 S24 

-Sl3 -Ci3S23e*''^ C23C13 

V — S14C13 5145135236*^^ — C14S24C23 ~5i4Si3C23 — 014524523^ C14C24 / 

(50) 

with the same angle and phase rotation convention adopted as before. 

In this one-mass-scale-dominant approximation with the large mass gap la- 
beled 5M'^ — 5m\gj^jj, the oscillations are again CP-conserving, and a short 
baseline experiment is needed to determine the extra relevant mixing angles and 
phase. The oscillation probabilities of interest are: 

P(^e - ^m) = ^(^M - ^e) = 4cf3ci,ci3.i3 sin^ (^) , 
P(^e ~^ ^t) — 4CI3CI4 [(5i3S^4s|3 -|- 014(^3824) 

-2C14S14C23S23S13S24 COS ^3] sin^ (^^) , (51) 
Piu, ^u,)^ 1 - 4c?3ci3(5i3 + sl,43) sin^ (^) , 

P{^H ~^ ^t) = 4C13C23 [(S13S14C23 + C]^4S23'^24) 

-F2Ci4Si4C23S23Sl3S24 COS 53] sin^ (^^^ 

If the neutrinos propagate through matter, these expressions must be mod- 
ified. Matter effects for the three active and one sterile neutrino scenario are 
similar in nature to those for the three active neutrino case, Eq. (35). However 
in Eq. (35) a flavor diagonal term that only contributes to an overall phase has 



34 



been discarded. This term comes from the coherent forward scattering amph- 
tude for the active flavors scattering from the electrons, protons and neutrons 
in matter via the exchange of a virtual Z-boson. Since the sterile neutrino does 
not interact with the Z-boson this term must be added to the diagonal terms for 
the active neutrinos (or cquivalently subtracted from the diagonal part for the 
sterile neutrino). That is in Eq. (35) 

AAA' 

-SaeSpe ^ TTTT^ae^/Je " TTFT^asSps (52) 



O ZT* ^ O IT* ^ O IT* 

Z£/^ ZH/i, ^l-Jv 

where A is given by Eq. (36) with replaced by — 1(1 — Ye) for electrically 
neutral matter. 

In order to search for CP violation, at least two mass scales must be relevant. 
For simplicity consider 

5m2i = 0, 5m|3 = 5m^, 

(5mi2= (5mii = 5M2, (53) 

with five angles and two phases present, since U\i{d\i^b\) = I. The CP-odd 
parts of the relevant probabilities are: 

PcP-odd(i^e ^ ^m) = 8c?3ci3C24C34S24S34 Sin(52 + Ss) (^) siu^ (^) 
-PcP-odd(i^e T^t) = 4C23C24 {2Ci4Si4C23S23Si3S24(s?3Sm " ^^4) Sn\{52 + 5^) 

+C14C34S13S14S34 [(si3 - sl^) s\u52 + s^3s|4sin(52 + 253)] 

+ Ci4C24Sl3Sl4S23S24(ci4 - SinS^} 

x(*Si).in=(-fii) 

-PcP-odd('^,. I^t) = 8C^3C23C24C34S34 [C14C23S13S14 siu 5^ + C14S23S24 sin(52 + S^)] 

(54) 

where only the leading order term in Sm^ has been kept. The CP-even expres- 
sions also have such additional small corrections. 

The present atmospheric neutrino data favors the — > i/t oscillation over the 
i/^ — > oscillation. On the other hand, if a solar neutrino oscillates significantly 
into a sterile neutrino, only the small angle MSW solution is viable since the 
large angle solutions fail to provide enough u + e~ —>■ u + e~ elastic scattering to 
be consistent with SuperK measurements [29] . Hence if it turns out that one of 
the large angle mixing solutions is the correct solution to the solar anomaly then 
something other than a single light sterile neutrino will be needed to explaining 
the solar, atmospheric and LSND results. 
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Table 4: Parameters for the four-flavor oscillation scenarios defined for the study. 
Note that for these parameter sets Sml^ ~ ^n^li ~ Sml2 ~ = ^M'^, and 

sin^ 2^14 = sin^ 2^13 = sin^ 2^24 = sin^ 2^23 



parameter 


llAl 


llBl 


Sml^ (eV^) 
Sml^ (eV2) 
SM^ (eV2) 
sin^ 2^34 


3.5 X 10^^ 


3.5 X 10-3 


6 X 10-6 
0.3 


6 X 10-6 
1.0 


1.0 


1.0 


sin^ 2^12 


0.006 


0.006 


sin^ 2^14 


0.03 


0.003 


Si 








S3 


0,±7r/2 



0,±7r/2 




3.1.5 Scenctrios with Three Active plus One Sterile Neutrino 

We now consider some representative four-flavor neutrino oscillation parameter 
sets that can be used to establish how well experiments at a neutrino factory 
could determine the oscillation parameters. As was noted earlier, the only viable 
solutions with one sterile and three active neutrinos require that there be two 
sets of almost degenerate neutrinos separated by the largest 5m^. We begin 
by considering the constraints from CHOOZ and LSND. Note that the effective 
two-component atmospheric and solar mixing angles are: 

sm^29atm= 4|[/^3p|^7^4|' = 44sin'2^^34 ^ ^ 

(55) 

sin^ 29soi = MUei\^\Ue2\^ = 44 sin' 2^12 
The CHOOZ constraint [27] from P(Pe ^ ^e) is: 

4 sin' 2^24 + C24 sin' 2^23 < 0.2 (56) 
while the LSND constraint [9] from P(t'j^ i^e) is: 

10~3 < 44 sin' 2^23 < 10"' . (57) 

With this in mind, the parameter sets we have deflned are summarized in Table 4. 
They are: 

Scenario llAl - Low Mass LSND 
Scenario IIBl - High Mass LSND 
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3.2 Where will we be in 5-10 years ? 



In this section, we briefly discuss the prospects for currently operating, planned, 
or proposed experiments exploring neutrino oscillations. The discussion will be 
broken down according to the various oscillation modes. The current limits and 
the expected reach of some of the future experiments are summarized in Fig. 14, 
and Tables 5 and 6. 

3.2.1 i/^ ^ i/^, Us 

The evidence for z/^ disappearance in atmospheric neutrinos at SupcrK is con- 
vincing [34]. The preferred region of parameter space is given [35] by 10^^ < 
Srri^ < 10~^ eV^ at near maximal mixing (sin^ 2^ ~ 1). The inter- 
pretation of the atmospheric neutrino deficit is disfavored by the SuperK data 
and is ruled out by the CHOOZ [36] and Palo Verde [37] experiments. The two 
immediate issues are (1) the precise determination of Sm^ and sin^ 29, and (2) 
discrimination between z/^ z/^ and u^ ^ Ug- 

Future SuperK data will probably not shrink the currently preferred region of 
parameter space by very much. The precision with which 6m^ can be extracted 
from the observed event distributions depends on the precisions with which the 
event-by-event values of L and E are determined. The greatest sensitivity to 
Sm^ comes from the sample of events with values oi L/E corresponding to the 
region of the first oscillation maximum, which in practice are those events from 
neutrinos coming from approximately the horizontal direction. However, for 
these events the L precision is limited by the angular resolution of the detector. 

SuperK can discriminate between i/^ — > Vg and — > Ur by looking for matter 
effects and by measuring the number of NC interactions. The lack of evidence 
for matter effects in up-going muons (in both partially-contained events and in 
the NC-enriched multi-ring event sample) already disfavors at the 99% 

confidence level. This result is expected to become firmer in the future. In ad- 
dition, if there is a significant fraction of incident sterile neutrinos, there will be 
fewer NC events detected. The cleanest sample of NC events is the sample of 
events with a single detected vr''. By comparing the ratio of 7r° ("two-electron") 
events to z/g CC (single-electron) events, SuperK already has a statistically sig- 
nificant handle on the NC/CC ratio (±6% for 848.3 days livetime). However, 
the measurement is currently limited by large uncertainties on the NC single 7r° 
cross section (23 — 25% total systematic). K2K will measure this cross section in 
its near detector, and over the next few years this new information may produce 
the most dramatic improvement in Vt/us discrimination [28]. 

The next generation of long-baseline experiments have been designed to be 
sensitive to oscillations with parameters that correspond to the SuperK favored 
region. The currently running K2K experiment [39] will cover Srn^ > 2xl0~^ eV^ 
after 3-5 years of running, and MINOS [15] will cover 8m^ > 0.6 x 10~^ eV^ 
(both at 90% CL and for maximal mixing). These experiments are expected to 
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Figure 14: The current and expected limits at some of the future neutrino oscil- 
lation experiments. Note that different oscillation modes are shown together. 
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confirm the neutrino oscillation interpretation of the atmospheric neutrino data 
by about 2005. In addition to searching for z/^ disappearance, K2K can also look 
for a distortion in the neutrino energy spectrum using quasi-elastic events, and 
MINOS can study NC/CC event energy distributions. In 2005 two experiments, 
OPERA [40] and ICANOE [41], are expected to begin taking data at the CNGS 
beam. Both OPERA and ICANOE aim primarily at r-appearance and will cover 
drn? > 2 X 10~^ eV^ after about 5 years of running. 

Unless < 2 x 10"^ eV^ (allowed at 99% CL at SuperK), we expect to have 
a complete accelerator based experimental confirmation of atmospheric neutrino 
oscillations by 2010, and measurements of Srn^ and sin^2^ at the 0(10%) level 
(see Table 6). If Sm^ < 2 x 10"'^ eV^ some additional experiments will be 
necessary, such as MONOLITH (30 kt calorimeter) [42] or AQUARICH (novel 
1 Mt Water Ccrcnkov) [43] . These experiments may study atmospheric neutrinos 
and exploit good angular resolution to search for dips in the zenith angle (or 
L/E) distribution. 

3.2.2 i/^ <-> I^e 

For large drri^ ~ 1 cV^ suggested by the LSND experiment [9], Mini-BooNE [44] 
is expected to cover the entire preferred region of LSND parameter space with a 
wide safety margin. In the event of a positive signal, they plan to build another 
detector (BooNE) that will be able to measure the parameters with a precision 
0(10%) (see Table 6). Independent confirmation from ICANOE would also 
be expected. Should all of the experimental indications for oscillations (LSND, 
atmospheric, and solar) be confirmed we may be seeing evidence for the existence 
of sterile neutrinos. This would be a very exciting discovery, would raise many 
questions, and would require a new round of experiments. 

For Sm^ ~ 10~^-10~^ eV^ we expect some <-> Vf. mixing if the heavier of the 
two mass cigcnstatcs involved in the atmospheric neutrino oscillation contains 
any admixture of (i-e., if Ues 7^ 0). Current limits from CHOOZ [36] and 
Palo Verde [37] require l^/gsl < 0.1. SuperK by itself is unlikely to improve on 
this sensitivity. K2K can look for Ue appearance and improve the sensitivity 
to a finite jC/esI in some Srn^ range. MINOS and ICANOE are expected to be 
sensitive to sin^ 2^13 > O{10^^) in the Sni^ region of interest by searching for z/g 
appearance in their predominantly beams. At this time it is not clear what 
is the interesting range for sin^ 2^13. If this mixing angle is not too small then 
K2K/MIN0S/ICAN0E can make a first measurement. The baselines for these 
experiments are too short, and statistics will be too limited, to observe matter 
effects. For very small mixing angles, comparable with the Small Mixing Angle 
MSW solution for the solar neutrino deficit (see [45]), an order of magnitude 
improvement in sensitivity beyond these experiments is required to make a first 
observation of — > i/g oscillations. 
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3.2.3 i/e Vfj,, Vr, T^s 

Reactor and solar neutrino experiments can only look for these oscillations in 
the Ve disappearance mode. 

The SNO [46] detector should discriminate between z/g z/^,r and Pe — > 
solutions to the solar neutrino deficit by studying the distortion in the v^. 
energy spectrum and by measuring the NC/CC ratio. The spectral distortion 
should occur for the SMA solution and for some regions of the VAC solution. 
Borexino [47] (or possibly KamLAND) will study '''Be solar neutrinos, and should 
see day/night effects for the LOW scenario and seasonal effects for the VAC 
solution. The absence of the ^Be electron neutrino ffux would strongly suggest 
the SMA solution. There arc additional experiments proposed to study lower 
energy neutrinos (esp. pp)\ HELLAZ, HERON, LENS, etc (see [48] for a recent 
overview). KamLAND [49] will look for the disappearance of Pg from reactors 
with sensitivity down to 5m^ > 10~^ eV^ for large mixing angles. With all of this 
data in the next 5-10 years we should have convincingly tested whether or not 
any of the current neutrino oscillation solutions to the solar neutrino problem 
are correct. 

None of the solar neutrino experiments, however, discriminate between v^. 
Vfj^ and i^e ^T- MINOS, OPERA, and ICANOE can look for r appearance but 
cannot separate — > i^r from ^ v-j.. 

3.2.4 Summciry 

To summarize, in 5-10 years: 

(i) — > Vr-i^s- If the 5m^ associated with the atmospheric v^^ deficit exceeds 

~ 2 X 10'^ eV^ accelerator experiments will measure 5rn? and sin^2^ with 
precisions 0(10%). If 5m^ is less than ~ 2 x 10~'^ eV^ new experiments 
will be required to accomplish this in the 2010 era. 

(ii) ^ Ve- If the LSND oscillations are confirmed BooNE would measure the 
associated 5rn? and sin^ 26 with precisions 0(10%). However the oscillation 
framework (sterile neutrinos ?) might be complicated. If LSND is not 
confirmed and if sin^2^i3 > 10~^, the first evidence for a finite value of 
sin^ 2^13 would be expected at long baseline accelerator experiments. If 
sin^ 2^^i3 < 10~^ then z/^ z/g will not be observed in the accelerator 
experiments and new experiments with at least an order of magnitude 
improved sensitivity will be needed. 

(iii) Vf, u^, Ut-, Vg. Either one or none of the current solar neutrino deficit 
solutions will be remaining. If one survives, we will know whether the solar 
neutrino deficit is due to Uf. ^ ^s- If the ^ u^, i/^ mode is favored we will 
not be able to distinguish between i/e — > I'lx or i/g — > i^r. In addition, — > i^t 
will not have been observed at long baseline accelerator experiments. 
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(iv) Sterile neutrinos. If the LSND, atmospheric, and solar neutrino oscillation 
results are all confirmed we may be seeing evidence for the existence of 
sterile neutrinos. This would be a very exciting discovery ! Many new 
questions will arise requiring new experimental input. 

Finally, it is worthwhile considering the possibility that a conventional neu- 
trino beam and the corresponding detectors undergo significant upgrades within 
the coming decade. For example, a Fermilab proton driver upgrade might enable 
the acceleration of up to about a factor of four more beam in the Main Injec- 
tor, resulting in a corresponding increase of the NUMI beam intensity. With 
an additional factor of 2 - 3 increase in detector mass, the event samples might 
be increased by an order of magnitude. However, systematic uncertainties must 
also be considered. For example, there will be limiting systematic uncertain- 
ties on the measurements of 5777,32 and sin^ 2^23 with a MINOS-type experiment 
that arise from the uncertainties on the near/far detector CC reconstruction ef- 
ficiencies, backgrounds to CC events from NC interactions, and an assumed 2% 
flux uncertainty from the near/far detector extrapolation. These uncertainties 
would prevent the precise determination of the oscillation parameters, even in 
the limit of infinite statistics. The ultimate (infinite statistics) precision that 
could be achieved with a MINOS-type experiment is shown in Fig. 15. A vcry- 
long-basclinc neutrino factory experiment would be able to make very significant 
improvements to the precision with which Sml2 and sin^ 26*23 are determined. 
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Table 5: Experimental neutrino oscillation observations expected in the next 
5-10 years at accelerator based experiments. 



Scenario 


Experiment 


i/^ Disap. 
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I/e ^ Vr 
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Table 6: Neutrino oscillation mixing angle and leading Sw? measurements ex- 
pected in the next 5-10 years at accelerator based experiments. 



Scenario Experiment 


Parameter 

sin2 2^i2 sin2 2^23 sin2 2^i3 6 Sm'^ (eV^) 


ivZrv 

MINOS 

ICANOE 
OPERA 
BooNE 


ou/o ou/o 

io%t io%t 

13% 60% 11% 
20% 14% 


TAO ICOIC 

MINOS 

ICANOE 
OPERA 
BooNE 


ou/o ou/o 

io%t io%t 

13% 60% 11% 
20% 14% 


TAQ 'WOV 
lAo IVZIV 

MINOS 
ICANOE 
OPERA 
BooNE 


oO/o oU/o 

io%t io%t 

13% 60% 11% 
20% 14% 


TR1 KOK 

MINOS 
ICANOE 
OPERA 
BooNE 


ou/o ou/o 
10% 15% 
13% 11% 

20% 14% 
10% 10% 


ICl K2K 
MINOS 
ICANOE 
OPERA 
BooNE 


100% 100% 
10% 15% 
25% 5% 7% 
5% 7% 
10% 15% 10% 




sin2 2^23 sin2 2^34 dml^ (eV^) drnj^ (eV^) 


TTAl ICOK 

MINOS 
ICANOE 
OPERA 
BooNE 


ou /o ou /o 
10% 6% 
10% 13% 7% 11% 
30% 20% 30% 14% 
10% 10% 


IIBl K2K 
MINOS 

ICANOE 
OPERA 
BooNE 


30% 50% 
10% 6% 

50% 13% 50% 11% 
20% 14% 

10% 10% 



^ With sin^ 2^23 constraint from Super K. 
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Figure 15: The expected precision that could be achieved by a MINOS-hke ex- 
periment (low energy, baseline 732 km) in the limit of infinite statistics but with 
conservative estimates of systematic errors. The calculated sensitivities are based 
only on disappearance measurements. The oscillation parameters correspond to 
scenario lAl, and the regions of sensitivity shown are at 90% CL. 
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3.3 The neutrino factory oscillation physics program 

We now consider the program of neutrino oscillation measurements at a neu- 
trino factory in the era beyond the next generation of long baseline experiments. 
The main goals in this era are likely to be to precisely establish the oscillation 
framework, determine the pattern of neutrino masses, measure matter effects to 
confirm the MSW phenomenon, make precise measurements or place stringent 
limits on all of the mixing-matrix elements (and hence mixing-angles), and ob- 
serve or place stringent limits on CP violation in the lepton sector. A neutrino 
factory can address each of these goals: 

(i) Establishing the oscillation framework. This requires measuring as a func- 

tion of L/E, or putting stringent hmits on, all of the oscillation prob- 
abilities P(fe I'x) and Pii^fj. — > t^x)- The oscillation framework can 
be established by summing the probabilities (a) P(z/e ~^ ^e) + P{^e 
v^) + P(z/e ^ Vr), and (b) P(z/^ z/g) + P(z/^ v^) + P{v^ Vr)- In a 
three-fiavor mixing framework, both sums should be unity for all L/ E. If 
there are sterile neutrinos participating in the oscillations one or both of 
the sums will be less than unity. Part (b) of the test will almost certainly 
be made with conventional neutrino beams, although with a precision that 
will be limited by the statistics and by the uncertainty on the 

P{^^x ~^ ^e) measurement arising from the 0(1%) v^. contamination in the 
beam. Part (a) of the test, which includes the first observation of (or strin- 
gent limits on) z/g oscillations, can only be made with an energetic 
[Ey > 10 GeV) z/g (or z/g) beam, and will therefore be a unique part of the 
neutrino factory physics program. 

(ii) Determining the pattern of neutrino masses. The present experimental data 

suggests that, within a three-fiavor mixing framework, there are two neu- 
trino mass eigenstates separated by a small mass difference, and a third 
state separated from the pair by a "large" mass difference SM'^. What is 
unknown is whether there is one low state plus two high states, or two 
low states plus one high state. This can be determined by measuring the 
sign of 5M^. The only way we know of making this measurement is to 
exploit matter effects which, in a very long baseline experiment, alter the 
probabilities for oscillations that involve electron neutrinos; the modifica- 
tion being dependent on the sign of 5M^. In principle the measurement 
could be made using a conventional neutrino beam and measuring z/^ — > 
and z7^ — j> z/g transitions over a baseline of several thousand km. However, 
the 0(1%) z/g (z7g) contamination in the beam will introduce an irreducible 
background that is comparable to, or larger than, the z/g signal. In con- 
trast, at a neutrino factory it appears that the measurement can be done 
with great precision. Hence, determining the sign of 5M'^ and the pattern 
of neutrino masses would be a key measurement at a neutrino factory. 
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(iii) Measuring matter effects to confirm the MSW phenomenon. The same 
technique used to determine the sign of 5m|2 can, with sufficient statistics, 
provide a precise quantitative confirmation of the MSW effect for neutrinos 
passing through the Earth. The modification to P{i'e — ^ ^/i), for example, 
depends upon the matter parameter A (Eq. (36)). Global fits to appear- 
ance and disappearance spectra that arc used to determine the oscillation 
parameters can include A as a free parameter. The quantitative MSW test 
would be to recover the expected value for A. This measurement exploits 
the clean i/g — > i/^ signal at a neutrino factory, and would be a unique part 
of the neutrino factory physics program. 

(iv) Making precise measurements or placing stringent limits on all of the 
mixing-matrix elements. In practice the measured oscillation probability 
amplitudes are used to determine the mixing angles. If any of the angles 
are unmeasured or poorly constrained the relevant entries in the mixing 
matrix will also be poorly determined. At present there is only an upper 
limit on ^13, the angle that essentially determines the — > i/^ oscillation 
amplitude. A neutrino factory would provide a precise measurement of, or 
stringent limit on, this difficult angle. In fact, because all of the 

and fx oscillation amplitudes can be measured at a neutrino factory, 

global fits can be made to the measured spectra to provide a very precise 
determination of the mixing angles. This exploits the v^. component in the 
beam. Finally, it should be noted that it is important to test the overall 
consistency of the oscillation framework by determining the mixing angles 
in more than one way, i.e. by using more than one independent set of 
measurements. Clearly the beam is an asset for this check. 

(v) Placing stringent limits on, or observing, CP violation in the lepton sector. 

Most of the oscillation scenarios defined for the study predict very small 
CP violating amplitudes. An important test of these scenarios would be 
to place stringent experimental limits on CP violation in the lepton sector. 
The LMA scenario lAl might result in sufficiently large CP violating effects 
to be observable at a neutrino factory. The CP test involves comparing 
T^e with Ve Vn oscillatiou rates, possible at a neutrino factory be- 

cause backgrounds are very small. A search for CP violation in the lepton 
sector with the required precision cannot be done with a conventional neu- 
trino beam, and is therefore a unique part of the neutrino factory physics 
program. 

Note that it is the i^g (Ve) component in the neutrino factory beam that 
drives the oscillation physics program. A i/e beam would (a) enable a basic 
test of the oscillation framework that cannot be made with a i/f^ beam, (b) 
enable the first observation of (or stringent limits on) z/g Vr oscillations, 
(c) make a convincing determination of the pattern of neutrino masses that 
would be difficult or impossible with a conventional neutrino beam, (d) make a 
quantitative check of the MSW effect only possible with a neutrino factory beam. 
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(e) enable measurements or stringent limits on all of the (three-flavor) mixing 
angles with a precision that requires both v^. and beams, and (f) measure 
or put meaningful limits on CP violation in the lepton sector, which requires a 
signal purity only available at a neutrino factory. 

A neutrino factory operating in the next decade, after the next generation 
of long baseline experiments, would appear to be the right tool at the right 
time. However, before we can quantitatively assess how well a neutrino factory 
might realize the physics program we have listed, we must first understand the 
capabilities of neutrino detectors in the neutrino factory era. 

3.4 Detector considerations 

We would like to measure the oscillation probabilities Piv^ — > Vf)) as a function 
of the baseline L and neutrino energy E (and hence LjE) for all possible initial 
and final fiavors a and /5. This requires a beam with a well known initial fiavor 
content, and a detector that can identify the flavor of the interacting neutrino. 
The neutrinos interact in the detector via charged current (CC) and neutral 
current (NC) interactions to produce a lepton accompanied by a hadronic shower 
arising from the remnants of the struck nucleon. In CC interactions the flnal 
state lepton tags the flavor (/3) of the interacting neutrino. 

At a neutrino factory in which, for example, positive muons are stored, the 
initial beam consists of 50% z/g and 50% V^. In the absence of oscillations, the Vf. 
CC interactions produce electrons and the CC interactions produce positive 
muons. Note that the charge of the final state lepton tags the flavor (a) of the 
initial neutrino or antineutrino. In the presence of z/g z/^ oscillations the z/^ 
CC interactions produce negative muons (i.e. wrong-sign muons). Similarly, 
Vf. oscillations produce wrong-sign electrons, 17^ — > oscillations produce 
events tagged by a r"*" and i/g — > oscillations produce events tagged by a 
T~ . Hence, there is a variety of information that can be used to measure or 
constrain neutrino oscillations at a neutrino factory, namely the rates and energy 
distributions of events tagged by (a) right-sign muons, (b) wrong-sign muons, 
(c) electrons or positrons (their charge is difficult to determine in a massive 
detector), (d) positive r-leptons, (e) negative r-leptons, and (f) no charged 
lepton. If these measurements are made when there are alternately positive and 
negative muons decaying in the storage ring, there are a total of 12 spectra 
that can be used to extract information about the oscillations. Some examples 
of the predicted measured spectra are shown as a function of the oscillation 
parameters in Figs. 16 and 17 for a 10 kt detector sited 7400 km downstream of 
a 30 GeV neutrino factory. Clearly, the high intensity i/g, '^^i ^-nd beams 
at a neutrino factory would provide a wealth of precision oscillation data. 

The detectors required at a neutrino factory will have many similarities to 
the detectors that have been designed for the next generation of experiments at 
conventional neutrino beams. However, there are some important differences. 
First, we can anticipate more massive detectors. The sensitivity of a neutrino 
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Figure 16: Visible energy spectra for four event classes when 10^^/i~ decay in 
a 30 GeV neutrino factory at L = 7400 km. Black histogram: no oscillations. 
Blue dotted histogram: 5ml^ = 3.5 X 10"=^ eVVc^ sin^ ^23 = 1- Red dashed 
histogram: = 7 x 10"'^ eV^/c^, sin^ 623 = 1. The distributions in this figure 

and the following figure are for an ICANOE-type detector, and are from Ref. 50. 




Figure 17: Same as previous figure, but with positive muons circulating in the 
storage ring. The difference between the two figures is due to the different cross 
section for neutrinos and antineutrinos, and to matter effects. 
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factory oscillation experiment is proportional to the product of the detector mass 
and beam intensity. It is likely that the cost of increasing the MINOS detector 
fiducial mass (for example) by a factor of a few is smaller than the cost of 
increasing the neutrino factory beam intensity by a factor of a few. Therefore, 
we believe that it is reasonable to assume that detectors at a neutrino factory 
would be a factor of a few to a factor of 10 more massive than the generation of 
neutrino detectors presently under construction. Second, the presence of both 
neutrinos and antineutrinos in the same beam at a neutrino factory places a 
premium on measuring the sign of the charge of the lepton produced in CC 
interactions. Charge-sign determination may not be practical for electrons, but 
is mandatory for muons and highly desirable for r-leptons. Finally, a relatively 
low energy threshold for the detection and measurement of wrong-sign muons 
is very desirable. This is because high muon detection thresholds require high 
energy interacting neutrinos, and hence a high energy neutrino factory. Since 
the muon acceleration system at a neutrino factory is likely to be expensive, low 
energies are preferable. 

In the following sections we begin by considering general detector issues for 
the measurement of final state muons and r-leptons, and then consider some 
specific candidate detectors for a neutrino factory. Some of these detector types 
are quite new and are just beginning to be studied; for the more mature detectors 
the "neutrino" energy resolution, the signal efficiency, background rejection, and 
fiducial mass are discussed. 

3.4.1 Muon identification and measurement 

The detection and measurement of muons (especially those of opposite sign to 
the muons in the storage ring) is crucial for many of the key oscillation physics 
measurements at a neutrino factory. Before considering some specific neutrino 
factory detectors it is useful to consider more generally muon backgrounds and 
related issues. Background muons can be produced in NC and CC interactions 
by: 

(i) Pions or kaons from the hadronic shower that decay to produce a muon. 

(ii) Non- interacting pions which fake a muon signature (punch-through) . 

(iii) Charm meson production and muonic decay. 

A background muon event can be produced when a background "muon" of 
the appropriate sign is recorded in (a) a NC event or (b) a CC event in which the 
primary lepton has been lost. If the background muon has the same charge sign 

as that in the storage ring the resulting event will be a background for disap- 
pearance measurements, but more importantly, if it has the opposite sign then 
the event will be a background for wrong-sign muon appearance measurements. 

The integrated wrong-sign background fraction from the hadronic shower 
is shown in Fig. 18 as a function of the minimum muon energy accepted for 
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Steel/Scintillator and water detectors downstream of 20 GeV and 50 GeV neu- 
trino factories. The cliarm background comes from i/^ CC events where the 
primary muon was less than 2 GeV. The peak at low muon energies is from the 
hadron shower itself and from punch through, while the long tail is from shower 
particles decaying to muons. 




Figure 18: Background levels from punch through, pion/kaon decay, and charm 
backgrounds for 20 GeV (left) and 50 GeV (right) neutrino factories. The frac- 
tion of neutrino interactions that produce a wrong-sign muon background event 
is shown as a function of the minimum muon energy accepted. 

In general there are two different standards for background levels which are 
relevant: that of a disappearance experiment and that of an appearance ex- 
periment. Background estimates are not trivial, but if the backgrounds for a 
disappearance measurement arc at the one per cent level, then the uncertain- 
ties on those backgrounds can be expected to be small compared to the flux 
uncertainty. On the other hand, wrong-sign muon appearance measurement 
uncertainties are expected to be dominated by the statistics. An extremely ag- 
gressive background level requirement would be to have less than of the order of 
one background event. If there are several thousand CC events expected, then 
this would require a minimum background rejection factor of 10^. 

Backgrounds can be suppressed by imposing a minimum energy requirement 
on the measured muon. Figure 19 shows the effect of several different minimum 
muon energy cuts on a simulated oscillation signal observed in a steel-scintillator 
type detector at a 20 GeV muon storage ring, at a baseline length of 2800km 
[51]. A muon threshold energy of 4 GeV for example depletes the low energy 
part of observed measured "neutrino energy" distribution, degrading but not 
completely removing the information about the neutrino oscillation parameters 
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Figure 19: Reconstructed neutrino energy distribution for several different min- 
imum muon energy cuts for a 20 GeV ring. Result is from Ref. 51. 

that is encoded in the shape of the distribution. A 4 GeV threshold at a 20 GeV 
neutrino factory is probably tolerable. If higher thresholds are needed to reject 
backgrounds, then a higher energy neutrino factory is desirable. If a lower energy 
neutrino factory is to be viable, then lower muon thresholds are desirable. 

As is shown in Fig. 18, to get to a background level of 10~^ one would need 
a 5 (6.5) GeV muon momentum cut in Steel/scintillator (Water) for a 20 GeV 
muon storage ring, and a 10 (12) GeV muon momentum cut in Steel/Scintillator 
(Water) for a 50 GeV muon storage ring. Clearly more background rejection 
is desirable. Fortunately muons from hadron decay in the hadronic shower are 
likely to be more aligned with the shower direction than muons from the leptonic 
vertex of the CC interaction. This provides another handle on the background. 
A useful variable to cut on is the momentum of the muon in the direction trans- 
verse to the hadronic shower (pt). Figure 20 shows the generated pf distribution 
for background and signal events, with no cut on the final state muon momen- 
tum. Note that requiring p'f > 1 the background is extremely low, while the 
signal efficiency is high. The resolution with which pf is determined is detec- 
tor dependent, and for detectors with reasonable transverse and longitudinal 
segmentation is dominated by the hadronic energy resolution. 

3.4.2 r— lepton identification and measurement 

The detection and measurement of r-leptons is crucial for z/^ Vr and z/g 
measurements at a neutrino factory. Note that v^. oscillations will be of 

special interest since they will not have been previously observed. The Ve 
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Figure 20: Distributions of the square of the muon momentum component trans- 
verse to the hadronic shower direction for charged current events compared 
to background muons for a 20 and 50 GeV muon storage ring. 

signal can be separated from Vj- "background" if the sign of the r-lepton 

charge is measured. The majority of r-lepton decays produce either one charged 
track (electron, muon, of hadron) or three charged tracks (hadrons). There are 
two general techniques that can be used to identify r-leptons. The first technique 
exploits the one-prong and three-prong topologies, and uses kinematic cuts to 
suppress backgrounds. The second technique uses a detector with a high spatial 
resolution to look for the displaced vertex or kink resulting from r-lepton decay. 

The advantage of the displaced vertex or kink detection r-lepton technique 
is that the detailed r-lcpton decay is measured and background suppression is 
therefore large. The disadvantage is that detectors that have sufficient spatial 
resolution are necessarily less massive than coarse-grained detectors. 

The advantage of the kinematic technique is that a very massive detector 
can be used. If the r-leptons decay muonically (BR = 17%) a measurement of 
the muon charge-sign determines the sign of the r charge. However, there are 
substantial backgrounds that must be reduced. In the case of muonic r decays, 
the backgrounds are from (a) z/^ (or z/^) CC interactions which typically produce 
muons at high momentum and high pj, and (b) meson decays (discussed earlier) 
which are at low momentum and low p|. For r — > e decays, the main background 
comes from Ug and Pg CC interactions. Fortunately the undetected neutrinos 
from r decays result in a larger missing transverse momentum than expected for 
backgroTind events. Exploiting these kinematic characteristics the backgrounds 
can be reduced by a large factor. For example, for an ICANOE-type detector 
a background rejection factor of 200 has been estimated, with a corresponding 
signal efficiency of 30%. In the electron channel background can also come 
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Figure 21: Example of a Ue Charged current event from the full simulation of 
the ICANOE detector. 

from NC interactions which produce photon conversions or Dalitz 7r° decays. 
These backgrounds can be suppressed in detectors with good pattern recognition 
allowing conversions, for example, to be identified and rejected. The analysis of 
hadronic r decays requires the identification of the r decay product inside a jet. 
This can only be done with a detector having good pattern recognition. It has 
been demonstrated that with an ICANOE-type detector a background rejection 
factor of 200 can be expected for r — > 1 prong, r — p, and r — Svr decays, with 
a signal efficiency of 8%. 

3.4.3 A Liquid Argon neutrino detector 

We have studied the performance of a large Liquid Argon neutrino detector 
at a neutrino factory using the ICANOE monte carlo program. One ICANOE 
detector unit consists of a liquid argon TPC followed by a magnetic spectrometer. 

The Liquid Argon TPC has extremely fine granularity, producing bubble 
chamber like event images. Figure 21 shows an example of an electron neutrino 
charged current event-note the separation between the electromagnetic shower 
and the hadronic shower of the nucleon remnant. The TPC is instrumented with 
3 mm pitch wires which allow tracking, dE/dx measurements, electromagnetic, 
and hadronic calorimetry. Electrons and photons can be identified and their 
energies are measured with a resolution given by aE / E = 0.03/V^ © 0.01. The 
hadron energy resolution is given by cte/E = 0.2/a/E © 0.05. The magnetic 
spectrometer is primarily needed to measure muon energy and charge, but it is 
assumed that it will also be instrumented as a calorimeter to allow the hadron 
energy of showers which leak into the spectrometer to be correctly measured 
(albeit with worse resolution). The muon momentum resolution is expected to 
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be 20%. 

In the design we have simulated, the hquid Argon module is 18 m deep with 
a cross section of 11.3m x 11.3m. The active (total) mass of one Liquid Argon 
module is 1.4 kt (1.9 kt). The magnetized calorimeter module is 2.6 m deep 
with a cross section of 9m x 9m, and has a mass of 0.8 kt. It consists of 2 m of 
steel, corresponding to 7.4 Xint and 59 Xq, interleaved with tracking chambers. 
Four Super-Modules are assumed, yielding a total detector length of 82.5 m and 
a total active mass of 9.3 kt that is fully instrumented. 

ICANOE can reconstruct neutrino (and antincutrino) events of all active 
flavors, and with an energy ranging from tens of MeV to tens of GeV, for the 
relevant physics analyses. The unique imaging capabilities of the liquid argon 
TPC allow one to cleanly determine whether a given event is a CC event, a 
Ve CC event, or a NC event. 

For our studies the ICANOE fast simulation was used. Neutrino interac- 
tion events are generated, with a proper treatment of quasi-elastic interactions, 
resonance and deep-inelastic processes. The 4- vectors for all the particles gener- 
ated are smeared, according to the resolutions derived from the full simulation. 
Muonic decays of pions and kaons are also considered, for a proper wrong- and 
right-sign muon background treatment. Once a 2-GeV cut is placed on the out- 
going muon momentum, the background levels tend to be about 10~^ times the 
actual charged current event rate, and are dominated by meson decay in the 
hadronic shower. 

Examples of simulated oscillation signals in an ICANOE-type detector at 
a neutrino factory are shown in Figs. 16 and 17. More detailed results from a 
study of the sensitivity that might be achieved using an ICANOE-type detector 
are discussed in the oscillation measurements section of this report. 

3.4.4 A magnetized Steel/Scintillator neutrino detector 

Steel/Scintillator calorimeters have been used extensively in past neutrino ex- 
periments. Their performance is well understood and well simulated. Typically 
a magnetized Steel/Scintillator (MINOS-like) neutrino detector consists of iron 
plates interspersed with scintillator planes. To obtain transverse position in- 
formation the scintillator can be segmented transversely, or a separate detector 
system (e.g. drift chambers) used. Penetrating charged particles (muon candi- 
dates) can then be reconstructed. With a reasonable transverse segmentation, 
the transverse position resolution is dominated by multiple coulomb scattering. 
The detector performance depends primarily on its longitudinal segmentation. 
The segmentation needs to be fine enough to determine whether a charged track 
has penetrated beyond the region of the accompanying hadronic shower. If it 
has, then the penetrating track is a muon candidate. The muon momentum 
resolution is determined by the magnetic field and the thickness of the steel 
plates. 
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Neutrino CC and NC interactions have well defined signatures. In a MINOS- 
like detector NC interactions produce a hadronic shower reconstructed as a large 
energy deposition in a small number of scintillator units. A i/^ or 17^ CC inter- 
action will produce a muon in the final state, characterized by a long track in 
addition to the hadronic shower. These events can be identified provided the 
muon penetrates well beyond the hadronic shower. This imposes a minimum 
track-length, and hence minimum energy, requirement on muons that can be 
identified. If the muon is not identified the CC interaction will look like a NC 
event. A Vf, or z7g CC interaction, will produce an electron in the final state which 
cannot be resolved, so these events look similar to NC interactions. A z/^ or 
CC interaction will also look like a NC interaction unless the r-lepton decays 
muonically. 

To study the performance of a magnetized Steel/Scintillator detector at a 
neutrino factory we have considered a detector geometry similar to the CCFR/NuTeV 
calorimeter [52], but with the addition of a toroidal magnetic field of IT. The 
detector is constructed from 3 x 3 x 0.3 m^ modules (see Fig. 22). The 0.7 kt 
CCFR detector consists of 42 modules. A neutrino factory detector with a mass 
of 50 kt (10 X the MINOS detector) would require 3000 of these modules. The 
ultimate transverse size (and hence module mass) that is practical is probably 
determined by the largest size over which a large magnetic field can be generated. 
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Figure 22: Schematic of a CCFR/NuTeV calorimeter module. 



In the following we consider how well a magnetized Steel/Scintillator detector 
can identify and measure wrong-sign muon events at a neutrino factory. For our 
simulations, we used the parameterized Monte Carlo developed by the NuTeV 
collaboration,modified to include particle tracking in the magnetic field. The 
hadron energy resolution of this detector is described in detail in [52], and is 
approximately given by (Je/E = 0.85/a/E. The muon momentum resolution 
depends on the track length in the steel, and whether the muon is contained 
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within the detector. For muons which range out in the detector the effective 
momentum resolution is ap/P = 0.05, while for tracks which leave the fiducial 
volume of the detector the resolution is described by ap/P ~ OMcs^BdL-, where 
the angles 9mcs and 9BdL describe respectively the change in direction due to 
multiple scattering and curvature in the magnetic field. 

The simulation includes a detailed parameterization of the hadron-shower 
development, with the inclusion of charm production and tt, K decays (the data 
set on which the decay probability parameterization was tuned contained only 
muons with momentum higher than 4 GeV/c). Note that tt punchthrough was 
not included in the parameterization, but is expected to make only a small 
contribution to background muons above 4 GeV. 
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Figure 23: Reconstructed yU~ with respect to the shower direction for 20 GeV 
and 50 GeV /i^ decaying in a neutrino factory. The muons are required to have 
energies exceeding 4 GeV. 

To be conservative, and reduce the dependence of our study on low energy 
processes that may not be adequately described by the Monte Carlo program, in 
our analysis all muons with generated energy below 4 GeV are considered lost. 
Muons with track length in steel less than 50 cm past the hadronic shower are 
also considered lost. All other muons are assumed to be identified with 100% 
efficiency, and measured sufficiently well to determine their charge sign. For the 
background events we considered (i) all the vr, K decay events producing "wrong- 
sign" muons in NC interactions, and (ii) all the charm production and tt, K decay 
events producing "wrong-sign" muons in CC events where the primary muon was 
considered lost. To reduce the backgrounds, we cut on Pf. The reconstructed 
Pf distribution is shown in Fig. 23 for signal and background muons in a 10 kt 
detector 2800 km downstream of 20 GeV and 50 GeV neutrino factories which 
provide 10^° fi^ decays. The oscillation parameters corresponding to the LMA 
scenario lAl. As expected, background wrong-sign muons, tend to have smaller 
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Pi than genuine wrong-sign muons from the leptonic vertex. The reconstructed 
wrong-sign muon spectrum is shown in Fig. 24 for a 20 GeV storage ring before 
(top plot) and after (bottom plot) muon energy, track length and > 2 GeV^ 
cut were applied. Signal and background rates are summarized in Table 7. After 
the cuts the signal/background ratio is above 10 to 1 in scenario lAl for a 
detector 2800 km away, while 40 — 50% of the — > signal events are retained. 
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Figure 24: Reconstructed wrong-sign muons as a function of the muon energy 
for a fi+ 20 GeV ring. Top plot accepted events for the signal (z/g z/^-stars- 
and Ue ^ i't ^ IJ' + X-crosses) and the potential backgrounds (x) . The bottom 
plot shows the signal and the background after cuts. 



3.4.5 A Water Cerenkov detector 

Preliminary studies have explored the possibility of using a large water Cerenkov 
detector as a distant target for a neutrino factory beam. Traditionally this type of 
detector has been used for measuring much lower energy neutrinos than expected 
at a muon storage ring, but to date water Cerenkov neutrino detectors are the 
only existing neutrino detectors with masses already in the 50 kt range. Water 
is of course the lightest target material under consideration in this report, but 
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Table 7: Wrong-sign muon rates after all cuts for a 10 kt steel-scintillator detector 
downstream of a neutrino factory providing 10^° muon decays. The oscillation 
parameters correspond to scenario lAl. The loss of signal acceptance and the 
background rejection are due solely to the kinematic and reconstruction cuts. 

H Ring ^ ^ix 
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0.50±0.02 


4.5 X 10-* 


50 




55.2 
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0.7 


0.41±0.02 
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this type of detector has several advantages when extrapolating to large masses, 
namely (i) low cost target material, (ii) only the surface of a very large volume 
needs to be instrumented, and (iii) good calorimctry. A large volume guarantees 
containment of hadronic and electromagnetic showers (as well as muons up to a 
certain energy). The low density of the target and good angular resolution from 
the Cerenkov cone might yield an overall hadron angle resolution that is as good 
as or better than the corresponding resolution obtained with steel-scintillator 
calorimeters. 

Water Cerenkov devices as large as 50 kt (SuperK) are already in operation 
and is expected to continue data-taking for ten years or more. Therefore, the 
response of the existing SuperK detector at a baseline distance of 9100 km has 
been studied as a test case. Next generation detectors, up to 1 Mton in mass, 
are technically feasible and arc currently under consideration for proton decay 
and neutrino measurements, sited perhaps at the Kamioka mine or elsewhere. 

For this initial study, the primary question is the suitability of a water 
Cerenkov detector for the higher energy neutrino beam produced by a 10-50 
GeV muon storage ring. At these energies, the multiplicity of hadrons is greater 
than for typical atmospheric neutrino interactions, and event topologies are cor- 
respondingly more complex. Figure 25 shows the Cerenkov light produced in a 
typical neutrino event from a 50 GcV muon storage ring at the SuperK detec- 
tor: the circles in the display are estimates of the outgoing angles of different 
charged particles produced in the hadronic shower. Some particle identification 
is possible from the pulse-height information. Reconstruction software from the 
SuperK experiment must be further optimized to study the detector response to 
neutrinos from 10 GeV and 50 GeV muon storage rings. It is worth noting that 
neutrinos produced by a 50 GeV muon beam induce a large number of events 
in the material (rock) surrounding any detector (producing an entering muon), 
and for a SuperK sized device these events outnumber the those produced in 
the detector's water volume. Both contained and entering events have therefore 
been studied. 

The response of a detector the size of SuperK changes drastically as a func- 
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Figure 25: Simulated neutrino event from a 50 GeV muon storage ring in the Su- 
perKamiokande detector. The rings indicate where the reconstruction software 
found charged particles in the hadronic shower, as well as the exiting muon. 



tion of beam energy. At a 10 GeV neutrino factory, 57% of the muon CC events 
are fully contained in the inner water volume, whereas only 11% are fully con- 
tained at a 50 GeV neutrino factory. This large difference only exists for events 
containing penetrating muons; at 10 GeV (50 GeV) both v^. and NC events are 
contained greater than 98% (90%) of the time. The existing /i-like particle iden- 
tification algorithm works to produce a reasonably pure (89%) z/^-CC sample for 
fully contained events in the 10 GeV beam, but e-like events are a mixture of 
of z/g-CC, NC and t'^-CC contamination. Exiting and entering events are pure 
samples of z/^-CC simply because of their penetrating nature. The muon angular 
resolution (3°) is much less than the muon- neutrino angular correlation (15°). 
With 2 X 10^° decays at a 50 GeV neutrino factory and a baseline of 9100 km, 
approximately 200, 000 v^, CC events would be observed entering or exiting the 
current SuperK detector. Combined with muon charge identification this sample 
should be able to provide good oscillation measurements. 

Implementing charge identification in a water Cerenkov detector is not triv- 
ial. Two possibilities have been proposed: (i) several large water targets, each 
one followed by a thin external muon spectrometer, and (ii) a magnetic field 
introduced into the water volume itself. Although the first design would have 
lower geometrical acceptance and a higher muon energy threshold, it would pose 
much less of a problem for the phototubes since the magnetic field would pre- 
sumably be well-contained in the spectrometer. The second proposal could in 
principle have good low energy muon momentum acceptance, but the resolution 
on the muon and hadron shower angles might be compromised. 

For a magnetic field internal to the target, 0.5-1 kC is sufficient to visually 
determine the charge of a several meter- long (> 1 GeV) muon, but no automated 
algorithms have yet been developed. A number of conceptual magnet designs 
have been studied: solenoidal, toroidal, and concentric current loops in the center 
or at the ends of the detector. A detailed study of one particular design has shown 
that one can immerse the central volume of a SuperK sized detector in a 0.5 kC 
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magnetic field while leaving only a 0.5 G fringe field (which may be acceptable 
with shielding and/or local compensation) in the region of the PMTs. Many 
of the difficulties inherent in placing a field inside a water detector would be 
avoided if an alternative hght collector (insensitive to the field) were used. Work 
on magnet design and alternative light collection is ongoing, but the internal 
magnetic field option must be considered speculative at this point. 

The results we will describe in the remainder of this section are for a wa- 
ter Cerenkov detector with an external magnetic field, because neutrino event 
reconstruction is more straightforward to simulate and the spectrometer tech- 
nology is well-understood. Although the studies of this detector are very pre- 
liminary, they look promising and warrant further investigation. We have used 
a LUND/GEANT Monte Carlo program which uses as its geometry a 40 x 40 x 
100 m^ box of water, followed by a 1 m long muon spectrometer. This simula- 
tion can be used to study acceptance issues and background contamination for 
a range of geometries and storage ring energies. 

Figure 3.4.5 shows the geometrical acceptance for the box-hke water Cerenkov 
detector as a function of distance of the neutrino interaction vertex from the 
spectrometer, for CC z/^ events from 20 and 50 GeV storage rings. The loss in 
acceptance close to the spectrometer is due to rejection of events where there is 
more than one muon which traverses the spectrometer (where the extra muon 
comes from background processes). It is clear that for a 20 GeV muon storage 
ring one would want a muon spectrometer much more frequently than once every 
100 m. Of course, noting that steel has a density of 8 times that of water, the 
smaller the ratio of water thickness to steel thickness the more it approximates 
a magnetized steel/scintillator target interspersed with water volumes with fine 
granularity. 

Clearly more work is needed to optimize the design for this kind of detector, 
but it might be an inexpensive compromise between a coarse-grained sampling 
calorimeter and a very fine-grained liquid argon TPC. 

3.4.6 Specialized r-lepton detectors 

The measurement of r-lepton appearance in large mass neutrino detectors is 
challenging. There are several ideas that might lead to viable new r-appearance 
detectors within the next 5-10 years, and that might be suitable for use at 
a neutrino factory. We briefly describe three examples in the following: (i) a 
perfluorohexane Cerenkov detector, (ii) a hybrid emulsion detector, and (iii) a 
very fine-grained micro-strip gas chamber target. 

Consider first a Cerenkov detector filled with perfiuorohexane (6*6-^14), which 
has a density 1.7 times that of water. This has been proposed [53] for use in the 
CERN to Gran Sasso beamline. The detector geometry consists of several target 
volumes followed by short muon spectrometer modules. A 1 Ton perfluorohex- 
ane detector (with a very different geometry) exists at DELPHI. The r-lepton 
signature in this type of detector consists of a sparsely populated Cerenkov ring 
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Spectrometer Acceptance for Neutrinos (solid) and Anti — Neutrions (dashed) 
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Figure 26: Acceptance in a water target for charged current (solid) and 
(dashed) events in a 20 and 50 GeV storage ring, as a function of distance of the 
neutrino interaction vertex from the muon spectrometer. 



from the r before it decays, together with a more densely populated ring from 
the daughter muon. The two rings would have offset centers. Figure 27 shows a 
simulated quasi-elastic Ut event (no hadron energy) from this kind of detector. 
This technique would probably not work for events with high energy hadron 
showers because of the large number of charged particles that would result in 
overlapping rings near the initial r-lepton ring. 

Next consider a hybrid emulsion detector consisting of, for example, thin 
(~ 100/im) sheets of emulsion combined with low-density (~ 300yum) spacers. 
The signature for a r-lepton 1-prong decay would be a change in direction of 
the track measured before and after the spacer [54] . For charge identification the 
detector could be within a large magnetic field volume. With an emulsion track 
angular resolution of 2 mrad, a 5a charge-sign determination of a 10 GcV/c 
charged particle could be achieved with a 2 T field and a 1.2 mm thick spacer 
[55] . An ~ 20 kt hybrid emulsion detector of this type might consist of 20 kt of 
steel segmented into 1 mm thick sheets, and an equal volume of thin emulsion 
layers plus low density spacers. The resulting detector would fit into the ATLAS 
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Figure 27: Quasi-elastic z/^ event in a perfluoroHexane Cerenkov detector: the 
ring described by about eight hits on the left is from the tau before it decays. 

barrel toroid magnet, which has a magnetic field ranging from 2 to 5 Tesla [56]. 
A hybrid emulsion detector with an external downstream muon spectrometer will 
be used by the OPERA experiment, which is to be put in the CERN to Gran 
Sasso beam. The muon spectrometer will determine the charge sign for r ^ yU 
decays provided the muon reaches the spectrometer. According to the OPERA 
studies [57], with an average neutrino energy of 20 GeV the total efficiency for 
seeing the r decays is 29% (including the branching ratios). The efficiency is 
largely geometric and should not be compromised by the addition of a magnetic 
field, provided the bend in the spacer due to the magnetic field is much less than 
the "apparent bend" due to the r-lepton decay. 

Finally, consider a target consisting of a tracking chamber constructed from 
micro-strip gas chambers (MSGCs) and a low Z material (for example, nylon) in 
a large magnetic field volume. This would be a NOMAD [58]-like detector with 
a much larger 0(1 kt) fiducial mass and an improved spatial resolution. Because 
of the low Z of the material electrons can travel a long distance in the detector 
before showering, and with a high enough field their charge can therefore be 
measured. Although a kink is not seen, the tau decay could be distinguished 
kinematically. For example, nylon has a radiation length of 37 cm. With a B 
field of 1 Tesla and MSGC's every 10 cm one would have an 8cr measurement 
of a 50 GeV electron's charge. This idea is worthy of further consideration, 
particularly if the LSND signal is confirmed and lower-mass tau detectors are 
warranted. 
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3.4.7 Detector summciry 

In our initial studies we have simulated the performance of steel/scintillator and 
liquid Argon detectors at a neutrino factory. Results arc encouraging. These 
technologies could provide detectors with masses of order 10 kt (liquid Argon) to 
a few xlO kt (steel/scintillator) that yield good wrong-sign muon identification 
and adequate background rejection. Our simulations of the capabilities of water 
Cerenkov detectors at a neutrino factory are less advanced, but initial results 
are encouraging, and this detector technology might permit very large detector 
masses to be realized. Some relevant characteristics of steel/scintillator, liquid 
Argon, and water Cerenkov detectors are listed in Table 8. It is premature to 
choose between detector types at this early stage. However, some general points 
are worth noting: 

(i) We believe that a cost optimization of detector mass (cost) versus neutrino 

factory beam intensity (cost) will probably favor detectors that are at 
least a factor of a few to a factor of 10 more massive than, for example, 
ICANOE or MINOS. A detector mass in the range 10 kt to 50 kt does not 
seem unreasonable. 

(ii) The minimum energy a muon must have for good identification and mea- 
surement may well determine the minimum viable muon storage ring en- 
ergy. This threshold is a few GcV, and is detector technology dependent. 
With a steel/scintillator detector and a threshold of 4 GcV, for exam- 
ple, the minimum acceptable neutrino factory energy appears to be in the 
neighborhood of 20 GeV. 

In this initial study we have not comprehensively considered to what extent mas- 
sive detectors at a neutrino factory need to be deep imderground. It seems very 
likely that detectors with low detection thresholds (water Gerenkov and liquid 
argon) will need to be well protected from cosmic ray backgrounds, regardless 
of the neutrino factory energy. For the steel/scintillator detector, the cosmic 
ray backgrounds for charged current events with muons in them are likely to 
be small for a detector at the surface of the earth, but there will be substantial 
background to neutral current or z/g charged current interactions. Finally, we 
note that the development of a new generation of very massive detectors capable 
of identifying and measuring the charge-sign of muons and r-leptons, would be 
of great benefit to a neutrino factory. There is a possible area of mutual inter- 
est with the nucleon decay community in developing the technology for a really 
massive 1 Mton scale water Gerenkov detector. This possibility deserves further 
investigation. 
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Table 8: Comparison of detector parameters for candidate detectors at a neutrino 
factory. 



Characteristic 




Detector Technology 




Steel/Scint 


Liquid Argon TPC 


Water Cerenkov 


Resolutions of: 








Electron Energy 


50%/^ 


3%/v^ © 1% 


0.6©2.6%//E 


Hadron Energy 




20%//E © 5% 


20-30% 


Muon Energy 


5% 


20% 


20%t 


Hadron Shower Angle 




rad 

(each hadron) 




Muon Angle 


5% for 
50cm track 


.02© .21/^ 


3° 


Maximum mass 


50 kton 


30 kton 


IMton? 


What limits size? 




safety, tunnel 


tunnel 


Required Overburden** 


m 


50 m 


50-lOOm 


Analysis Cuts 


> 4 GcV 
P2 > 2 GeV^ 


> 2 GeV 




Background level 


10-^ 


2 X 10-5 





** The overburden required for all technologies depends on the neutrino factory 
duty factor. The overburden required for a steel-scintillator calorimeter also 
depends on the energy of the muon storage ring; but in the past this type of 
detector has been used at ground level with minimal contamination in the 
charged current sample above a neutrino energy of 5 GeV.^ The muon momentum 
resolution would be comparable to that of an ICANOE detector if the muon 
spectrometer were separated from the water tank volume. 
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3.5 Oscillation measurements 

Using the oscillation scenarios described in section 3.1 as examples, we can now 
assess how well the neutrino oscillation physics program outlined in section 3.3 
can be pursued at a neutrino factory with the detectors described in section 3.4. 
In the following sub-sections the oscillation measurements that can be made at 
a neutrino factory are discussed as a function of baseline, muon beam energy, 
and muon beam intensity. In particular we consider: 

(i) The first observation of z/g — > oscillations, the measurement of the sign of 

Sm^ and hence the pattern of neutrino masses (section 3.5.1). 

(ii) The first observation of Ue — > oscillations (section 3.5.2). 

(iii) The measurement of oscillations (section 3.5.3). 

(iv) Precision measurements of the oscillation parameters (section 3.5.4). 

(v) The search for CP violation in the lepton sector (section 3.5.5). 

The results are based on the calculations described in more detail in Refs. [59, 
51, 60, 50]. The calculations from Ref. [50] are for an ICANOE type detec- 
tor, and include realistic resolution functions, analysis cuts, and background 
modeling, but use a constant average matter density to compute matter effects. 
The calculations from Refs. [60, 51, 59] all use resolution functions typical of 
steel/scintillator detectors and, unless explicitly stated, reasonable thresholds 
on the detected muon energies. The calculation from Ref. [51] docs not include 
backgrounds but covers a broad range of scenarios, and uses the explicit trans- 
Earth density profile to compute matter effects. In contrast, the calculations 
from Refs. [60, 59] have been used to look at only a few oscillation scenarios, 
but include backgrounds and use respectively the average Earth density and the 
explicit density profile in computing matter effects. It should be noted that 
although there are significant differences in the details implemented in the cal- 
culations, in general all the four groups arrive at similar assessments for the 
measurement sensitivity at a neutrino factory as a function of energy, intensity, 
and baseline. 

3.5.1 Observation of v^. — > oscillations and the pattern of neutrino 
nicisses 

At a neutrino factory z/g — ^ oscillations would be signaled by the appearance 
of CC interactions tagged by a wrong-sign muon [2]. Within the framework 
of three-flavor oscillations the oscillation amplitude is approximately 

proportional to sin^2^i3. At the present time only an upper limit exists on 
sin^ 2^13. The next generation long-baseline oscillation experiments are expected 
to be able to improve the sensitivity to sin^^ia fa 10~^, i.e. about one order of 
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magnitude below the present bound. If sin^ 26*13 is in this range we would expect 
to observe Ve oscillations at a relatively low intensity neutrino factory, 

measure matter effects, and determine the pattern of neutrino masses [63] . This 
is discussed further in the remainder of this sub-section. 

sin^ 2-;?,3 yielding 1 /i" evts/ 1 O" jx* Decays 
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Figure 28: Reach in sin^ 2^13 for the observation of 10 /i~ events from i/g — 
oscillations, shown versus baseline for three 6171^2 spanning the favored SuperK 
range. The other oscillation parameters correspond to the LMA scenario lAl. 
The curves correspond to lO^^/i'^ decays in a 20 GeV neutrino factory with a 
50 kt detector, and a minimum muon detection threshold of 4 GeV. Results are 
from Ref. 51. 

It is useful to define [51] the sin^ 26*13 "reach" for an experiment as the value of 
sin^ 2^13 for which a given physics goal would be met. We take as our initial goal 
the observation of 10 z/g — > events tagged by a wrong-sign muon. Consider 
first the sin^ 26*13 reach for a 50 kt detector sited a distance L from a 20 GeV 
neutrino factory in which there are 10^^/i"'" decays in the beam-forming straight 
section. The baseline-dependent sin^ 2^13 reach is shown in Fig. 28 for a three- 
flavor oscillation scenario in which 5m\^, sin^ 2^12, and sin^ 2^23 correspond to the 
LMA scenario lAl, and the value of 5m|2 is varied over the favored SuperK range. 
Backgrounds are expected to be less than one event for L > 2800 km (Table 9), 
and are not included in the calculation shown in the figure. If ^mgg is in the center 
of the SuperK range, the sin^ 26*13 reach is about an order of magnitude below 
the currently excluded region, improving slowly with decreasing L. However, 
at short baselines (L < 2800 km) backgrounds may degrade the sin^ 2^13 reach. 
The reach improves (degrades) by a about a factor of 2 (3) if 6171^2 is at the 
upper (lower) end of the current SuperK range. If the oscillation parameters 
correspond to the LMA scenario lAl (sin^ 2^13 = 0.04), then only 2 x 10^^ muon 
decays are required at a 20 GeV neutrino factory to observe 10 signal events 
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in a 50 kt detector at L = 2800 km. The calculation [51] assumes that CC 
events producing muons with energy less than (greater than) 4 GeV are detected 
with an efheiency of (1). The number of muon decays needed to observe 10 
events is shown in Fig. 29 as a function of for the LMA scenario lAl, 
the SMA scenario IA2, and the LOW scenario IA3. The required muon beam 
intensities decrease with increasing E^, and are approximately proportional to 
E~^-^. Compared with the SMA and LOW scenarios, slightly less intensity is 
needed for the LMA scenario, showing the small but finite contribution to the 
signal rate from the sub-leading Srn^ scale. In all three scenarios (LMA, SMA, 
LOW) a 20 GeV neutrino factory providing 10^^ decays in the beam-forming 
straight section would enable the first observation of z/g oscillations in 

a 50 kt detector provided sin^ 26*13 > 0.01. It should be noted that although 
sin^ 2^13 could be very small, there are models [61] that predict sin^ 2^^i3 ~ 0.01. 

Having established Ve ~^ oscillations, further data taking would facilitate 
the measurement of matter effects and the determination of the sign of 5m^, 
and hence the pattern of neutrino masses. To illustrate the effect of matter on 
the z/g oscillation probability, the predicted measured energy distributions 

2800 km downstream of a 30 GeV neutrino factory are shown in Figs. 30 and 31 
for respectively — > and — > I^^ wrong-sign muon events. The distributions 
are shown for a range of positive and negative values of (5m|2- Note that for a 
given 1 577132!, if Sml2 < we would expect to observe a lower wrong-sign muon 
event rate and a harder associated spectrum when positive muons arc stored in 
the neutrino factory than when negative muons are stored. On the other hand, 
if 57TT.32 > we would expect to observe a higher wrong-sign muon event rate 
and a softer associated spectrum when positive muons are stored in the neutrino 
factory than when negative muons are stored. Hence, measuring the differential 
spectra when positive and negative muons are alternately stored in the neutrino 
factory can enable the sign of ^mgg to be unambiguously determined [63]. 

The expected number of wrong-sign muon events are listed in Table 9 for 
the LMA scenario lAl, and a 50 kt detector downstream of a neutrino factory 
providing 10^^//"*" decays and the same number of decays. The event rates are 

shown for both signs of Smf^2j for various storage ring energies and baselines. 
Even at a 20 GeV neutrino factory the signal rates at L = 7332 and 2800 km 
are large enough to permit the sign of Sml2 to be determined with a few years 
of data taking. 

We conclude that for the LMA, SMA, and LOW three-flavor mixing sce- 
narios we have considered, a 20 GeV neutrino factory providing 10^^ decays 
in the beam-forming straight section would be a viable entry-level facility. In 
particular, with a 50 kt detector and a few years of data taking cither z/g — > z/^ os- 
cillations would be observed and the sign of 6ml2 determined or a very stringent 
upper limit on sin^ 26*13 will have been obtained (discussed later). Long baselines 
(> 2000 km) are preferred. The longest basehne we have considered (7332 km) 
has the advantage of lower total event rates and hence lower background rates. 
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Figure 29: The required number of muon decays needed in the beam-forming 
straight section of a neutrino factory to achieve the physics goals described in the 
text, shown as a function of storage ring energy for the LMA scenario lAl, SMA 
scenario IA2, LOW scenario IA3, and a bimaximal mixing scenario BIMAX. 
The basehne is taken to be 2800 km, and the detector is assumed to be a 50 kt 
wrong-sign muon appearance device with a muon detection threshold of 4 GeV 
or, for z/g — >■ appearance, a 5 kt detector. Results are from Ref. 51. 
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Figure 30: Predicted measured energy distributions for CC events tagged by a 
wrong-sign (negative) muon from i/g — ^ i^ti oscillations (no cuts or backgrounds), 
shown for various 5m|2, as labeled. The predictions correspond to 2 x 10^° 
decays, Ef^ = 30 GeV, L = 2800 km, with the values for (5m^2) sin^ 26^13, sin^ 2^23, 
sin^ 2^12, and 6 corresponding to the LMA scenario lAl. Results are from Ref. 63. 
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Figure 31: Same as previous figure, for CC events tagged by a wrong-sign (pos- 
itive) muon from Ue —>■ i^n oscillations. 



69 



Table 9: Wrong-sign muon rates for a 50 kt detector (with a muon threshold of 
4 GeV) a distance L downstream of a neutrino factory (energy E^) providing 10^^ 
muon decays. Rates are shown for LMA scenario lAl with both signs of 8171^2 
considered separately. The background rates listed correspond to an assumed 
background level of 10~^ times the total CC rates (see section 3.4) with no 
energy dependence. Energy dependent cuts might suppress backgrounds further. 
Results are from Ref. 63. 





L 




ji^ stored 






^ stored 




GeV 


km 


(5m32 > 


5777.32 < 


Backg 


577732 > 


577732 < 


Backg 


20 


732 


52. 


36. 


7.3 


32. 


26. 


6.5 




2800 


46. 


9.2 


0.43 


7.1 


26. 


0.36 




7332 


33. 


0.97 


0.063 


0.55 


19. 


0.05 


30 


732 


100. 


72. 


25. 


58. 


45. 


24. 




2800 


90. 


26. 


1.6 


19. 


43. 


1.4 




7332 


43. 


3.3 


0.19 


2.1 


33. 


0.17 


40 


732 


150. 


110. 


60. 


83. 


65. 


58. 




2800 


140. 


48. 


4.0 


36. 


64. 


3.8 




7332 


54. 


5.6 


0.49 


3.1 


28. 


0.43 


50 


732 


200. 


140. 


120. 


110. 


84. 


120. 




2800 


180. 


71. 


7.9 


53. 


82. 


7.7 




7332 


56. 


8.0 


1.1 


5.0 


34. 


1.0 



3.5.2 Observation of oscillations 

We begin by considering the LMA scenario lAl, and ask: What beam intensity 
is needed to make the first observation of — > Vt oscillations in a detector that 

is 2800 km downstream of a 20 GeV neutrino factory ? The z/g and the 

accompanying event rates are shown in Fig. 32 as a function of the 

oscillation parameters sin^ 26'i3 and 577732. The v^. signal rate is sensitive 

to both of these parameters, and hence provides an important consistency check 
for three-flavor mixing: the observation or non-observation of a ^ v^. signal 
must be consistent with the oscillation parameters measured from, for example, 
~^ ^11, and disappearance measurements. For the LMA scenario 

lAl the observation of 10 signal events in a 5 kt detector (with 30% r-lepton 
efficiency) would require 3 years with 7 x 10^^//''" decays per year in the beam 
forming straight section. Very similar beam intensities arc required for the SMA 
and LOW scenarios (IA2 and IA3). Note that, over the sin^ 2^13 range shown in 
Fig. 32, the — >■ Vr rates are one to two orders of magnitude higher than the 
Vf. —>■ Vt rates. Hence, we will need a detector that can determine the sign of the 
tau-lepton charge at the 2(7 — 3(7 level, or better. 

Let us define the sin^ 2^13 "reach" for an experiment as the value of sin^ 2^13 
for which we would observe 10 — > events when there are 10^*^ muon decays 
in the beam-forming straight section of a neutrino factory. The sin^ 2^13 reach 
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Figure 32: Vr CC appearance rates in a 5 kt detector 2800 km downstream of a 
20 GeV neutrino factory in which there are lO^'^/x"'" decays in the beam-forming 
straight section. The rates are shown as a function of sin^ 2^13 and 5ml2 with 
the other oscillation parameters corresponding to the LMA scenario lAl. The 
top 3 curves are the predictions for events and the lower curves are for 

z/g events. Results are from Ref. 51. 
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Figure 33: Reach in sin^ 2^13 for the observation of 10 Ur oscillation events, 
shown as a function of baseline for four storage ring energies. The oscillation 
parameters correspond to the LMA scenario lAl. The curves correspond to 10^° 
/i"*" decays in a 20 GeV neutrino factory with a 5 kt detector. Results are from 
Ref. 51. 



is shown as a function of the baseline and storage ring energy in Fig. 33 for a 
5 kt detector and an oscillation scenario in which all of the parameters except 
sin^2^i3 correspond to scenario lAl. The reach improves with energy (approxi- 
mately ~ E~^'^ [51]) and is almost independent of baseline except for the highest 
energies and baselines considered. However, backgrounds to a i/g — ^ J^t oscillation 
search have not been studied in detail, and are not included in the calculation. 
Background considerations will favor longer baselines. The number of muon de- 
cays needed to observe 10 Ve ~^ events is shown as a function of muon beam 
energy in Fig. 29 for the LMA scenario lAl, the SMA scenario IA2 , and the 
LOW scenario I A3. We conclude that within these three-flavor mixing scenar- 
ios, a 20 GeV storage ring in which there are O(10^°) muon decays per year 
would begin to permit an observation of, or meaningful limits on, Ue Ur os- 
cillations provided a multi-kt detector with good tau-lepton identification and 
charge discrimination is practical. 

Next, consider the oscillation scenarios IBl (atmospheric + LSND scales) and 
ICl (three-flavor with atmospheric, solar, and LSND data stretched). In these 
cases the leading is large (0.3 eV^/c'^) and medium baseline experiments 
(L = 10-100 km) become interesting. As an example, consider a medium baseline 
experiment a few xlO km downstream of a 20 GeV neutrino factory in which 
there are lO^^yU"*" decays. The z/g and accompanying z7^ Vr event rates 
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Figure 34: Ur CC appearance rates in a 1 kt detector downstream of a 20 GeV 
neutrino factory in which there are lO^^/i"*" decays. Rates are shown as a function 
of the basehne L and phase 5, with the other oscillation parameters correspond- 
ing to the LSND + Atmospheric scenario IBl. Predictions for — > Vt and 
17^ are shown separately, as labeled. Results are from Ref. 62. 



are shown in Fig. 34 as a function of the baseline and the phase 5 with the other 
oscillation parameters corresponding to scenario IBl. In contrast to the — >■ 
rates, the z/g — >■ Vr rates are very sensitive to 5, and for \5\ > 20° can be very 
large, yielding thousands of events per year in a 1 kt detector at L = 60 km, 
for example. Note that the corresponding z7^ Vj. rate is of order 100 events. 
For small \5\ the Vj. rate will dominate the r appearance event sample. 

For larger \5\ the Vf. — rate dominates. Good r charge determination will 
therefore be important to measure both z/g — > Vr and — »• Vt oscillations. 

Now consider the r appearance rates in scenario ICl. In this case the rates 
are not sensitive to 5 and, for a 1 kt detector at L = 60 km, there are about 8000 
z/g — > events and 93000 — > Vr events [62]. A detector with 3cr (or better) 
r-lepton charge discrimination would enable these two rates to be separately 
measured. 

We conclude that measurements of the Ve — *■ i^r oscillation rate at a neutrino 
factory would provide an important test of the oscillation scenario. In LMA, 
SMA, and LOW three-flavor oscillation scenarios, a 20 GeV neutrino factory 
providing O(10^°) muon decays could permit an observation of, or meaningful 
limits on, z/g — > z/^ oscillations. In LSND-type scenarios where the leading dm? 
scale is large, a 20 GeV neutrino factory providing 0(10^^) muon decays might 
already permit hundreds of z/g — events to be measured. It should be noted 
that the feasibility of a multi-kt detector with good r identification efficiency 
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(for example 30%) and good charge sign determination has not been explored 
in detail at this stage, and further work is required to identify the best detector 
technology for this, and determine the expected resolutions and efficiencies. 

3.5.3 Measurement of — > oscillations 

The present SuperK data suggests that the atmospheric neutrino deficit is due 
to — > oscillations. If this is correct the next generation of accelerator based 
long baseline experiments are expected to measure these oscillations. Neverthe- 
less, for a fixed neutrino factory energy and baseline, it is important to measure 
or put stringent constraints on all of the appearance channels so that the sum of 
the appearance modes can be compared with the disappearance measurements. 
Hence, we briefly consider v-j- rates at a neutrino factory. Note that at 

a 20 GeV neutrino factory the average interacting neutrino energy is of order 
15 GeV, and for 5m^ within the favored SuperK range, the first oscillation max- 
imum occurs at baselines of 7000 ± 3000 km. At shorter baselines the oscillation 
probabilities are lower and hence the signal/background ratio is lower, although 
the signal rate can be higher. 

Consider first a 5 kt detector 2800 km downstream of a 20 GeV neutrino 
factory in which there are 10^° muon decays. The — > Vt event rates are 
shown in Fig. 32 as a function of sin^2^i3 and 5m\2, with the other oscillation 
parameters corresponding to the LMA scenario lAl. If negative muons arc stored 
in the neutrino factory, the resulting event rates would be about a factor 

of two higher than the 17^ Vr rates shown in the figure. A neutrino factory 
providing O(10^°) muon decays would enable i/^ — >• appearance data samples 
of a few hundred to a few thousand events to be obtained. Similar rates are 
expected in SMA and LOW three-flavor mixing scenarios. 

Next consider a longer baseline example in which a 10 kt ICANOE-type 
detector is 7400 km downstream of a 30 GeV neutrino factory which provides 
10^° muon decays in the beam forming straight section [50]. The main advantage 
of a longer baseline is that the total interaction rate, and hence the r-lepton 
background, is reduced. The energy distribution for events in which there is no 
charged lepton can directly reflect the presence of a signal (see Fig. 16). The 
non-r events in this event sample can be suppressed using topology-dependent 
kinematic cuts. It is desirable that the r charge-sign also be determined which, 
with an external muon spectrometer, will be possible for the r — > subsample. 
We conclude that the measurement of — > v-,- oscillations with high statistical 
precision will be possible at a neutrino factory in long and very long baseline 
experiments. A more complete study is warranted. 

3.5.4 Determination of sin^2^i3, sin^ 2^235 bmk^^ 

Consider first the determination of sin^ 2^13. The most sensitive way to measure 
sin^ 2^13 at a neutrino factory is to measure the — > oscillation amplitude. 
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which is approximately proportional to sin^26'i3. More explicitly, the value of 
sin^ 26*13 is extracted from a fit to the spectrum of CC interactions tagged by a 
wrong-sign muon. Background contributions from, for example, muonic decays 
of charged mesons must be kept small, which favors small total event samples 
and hence long baselines. 
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Figure 35: Visible energy spectrum for events tagged by wrong-sign muons 
in an ICANOE-type detector (full histogram). The oscillation parameters are 
6ml2 = 3.5 X 10"^ eV^/c^, sin^ 6'23 = 1, and sin^ 26^13 = 0.05. Also shown are 
the contributions from oscillations (black dashed curve), v^. v^, with 

a subsequent muonic decay of the r lepton (red curve), and background from 
muonic decays of pions or kaons in neutral current or charged current events 
(blue dot-dashed curve). Results are from Ref. 50. 

As an example, consider a 10 kt ICANOE-like detector that is 7400 km 
downstream of a 30 GeV neutrino factory [50]. The simulated energy spectrum 
of wrong-sign muon events is shown in Fig. 35 for three-flavor oscillations with 
the parameters dm^^ = 3.5 x 10"^ eV^/c"^, sin^2^23 = 1, and sin^2^i3 = 0.05. 
Note that the backgrounds predominantly contribute to the low energy part of 
the spectrum. To fit the observed spectrum and extract sin^ 26*13 niatter effects 
must be taken into account. The modification of the oscillation probability due 
to matter effects is a function of the profile of the matter density p between the 
neutrino source and the detector. The density profile is known from geophysical 
measurements, and this knowledge can either be used in the fit, or alternatively p 
can be left as a free parameter. It has been shown that both methods give consis- 
tent results [50] , and that the uncertainties on the fitted values of p and sin^ 26*13 
are not strongly correlated. However, the fitted value for sin^ 26*13 does depend 
on the assumed values for sin^ 6*23 and 5m\2- The measured right-sign muon (z/^ 
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Figure 36: Results from a global fit to the visible energy distributions for various 
event classes recorded in a 10 kt ICANOE-type detector 7400 km downstream 
of a 30 GeV neutrino factory. The 68% CL contours correspond to experiments 
in which there are 10^^, 10^°, and 10^^/i"'' decays in the neutrino factory (as 
labeled) followed by the same number of /i~ decays. Upper plot: density fixed 
to its true value. Lower plot: density is a free parameter of the fit. Results are 
from Ref. 50. 



disappearance) distribution, together with the distributions of events tagged by 
electrons, r-leptons, or the absence of a lepton, can be used to constrain these 
additional oscillation parameters. Hence, the best way to extract sin^ 26'i3 is 
from a global fit to all of the observed event distributions, with the oscillation 
parameters (and optionally p) left as free parameters. If the density profile is 
left as a free parameter, the fit determines its value with an uncertainty of about 
10% [50, 59]. This provides a quantitative test of the MSW effect ! Examples of 
fit results in the (sin^ 2^13, sin^ ^23)^plane are shown in Fig. 36 for 10^^, 10^°, and 
10^^ muon decays in the neutrino factory. As the beam intensity increases the 
measurements become more precise. With lO^^/x"'" and /i~ decays sin^ 2^13 and 
sin^ 26'23 are determined with precisions of 40% and 20% respectively. With 10^^ 
decays these precisions have improved to ~ 5%. If the baseline is decreased from 
7400 km to 2900 km the oscillation parameters are determined with compara- 
ble (although slightly worse) precisions (Fig. 37). We conclude that within the 
framework of three-flavor mixing, provided sin^ 2^13 is not too small, a global 
flt to the observed oscillation distributions would enable sin^2^i3, sin^6'23, and 
5m32 to be simultaneously determined, and the MSW effect to be measured. 

Consider as a second example a 20 kt MINOS-type detector 2800 km down- 
stream of a 30 GeV neutrino factory providing 10^° muon decays. Some prelim- 
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Figure 37: Results from a global fit to the visible energy distributions for various 
event classes recorded in a 10 kt ICANOE-type detector downstream of a 30 GeV 
neutrino factory in which there are 10^°//"'" decays in the neutrino factory followed 
by the same number of /i~ decays. The 68% CL contours correspond to baselines 
of 7400 km and 2900 km, as labeled. Results are from Ref. 50. 

inary fit results [60] for a LMA type scenario are shown in Fig. 43b. For this 
example the analysis required < 4 GeV/c, but did not use the full set of cuts 
described in section 3.4, and therefore tolerated a background level a factor of 
a few greater than shown in Fig. 24 and Table 7. Nevertheless, the fits to the 
measured distribution of energies for events tagged by wrong-sign muons were 
able to determine sin^ 2^13 and Sml2 with precisions of 14% and 10% respectively 
for scenario lAl. 

To illustrate the ultimate sensitivity to the oscillation parameters that might 
be achievable at a high intensity neutrino factory, consider next a 40 kt Fe- 
scintillator detector downstream of a 50 GeV neutrino factory in which there 
are lO^^/i"*" decays followed by 10^^/i~ decays [59]. Fit results in the (matter 
density, sin^ 2^i3)-plane are shown in Fig. 38 for three baselines. The precision 
on the sin^ 2^13 determination is a few percent. Note that the analysis described 
in Ref. [59] suggests that backgrounds can be suppressed to less than 10~^ of 
the total CC rate in the detector. This impressive level of background rejection 
deserves further study. At the shortest baselines (732 km) matter effects are too 
small to obtain a good determination of the matter density parameter. 

Consider next the precision with which the oscillation parameters can be 
determined if sin^ 26'i3 is very small, and hence no u^, —> oscillation signal is 
observed. The resulting limits on sin^ 2^13 are shown as a function of 5m\2 in 
Fig. 39 for a 10 kt ICANOE type detector 7400 km downstream of a 30 GeV 
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Figure 38: Fit results in the (A, ^^i3)-plane for a simulated experiment in which 
a 40 kt Fe-scintillator detector is a distance L km downstream of a 50 GeV 
neutrino factory in which there are lO^^/j, decays. The density parameter A is 
defined in Eq. (36). The curves are 68.5, 90, and 99% CL contours. Results are 
from Ref. 59. 

neutrino factory in which there are 10^°;U+ decays followed by 10^'^ /i~ decays [50]. 
The resulting upper limit on sin^ 2^13 would be 0(10^^ — 10^'^), about three orders 
of magnitude below the present experimental bound, and one to two orders of 
magnitude below the bound that would be expected at the next generation of 
long-baseline experiments. 

The limit would become even more stringent at a higher intensity neutrino 
factory. As an example of the ultimate sensitivity that might be achievable, in 
Fig. 40 the limits on sin^ 29i3 are shown as a function of ^^'^ baseline for a 
40 kt Fe-scintillator detector downstream of a 50 GeV neutrino factory in which 
there are 10^^/x decays [59]. The non-observation of Ue — > oscillations could 
result in an upper hmit on sin^ 2^13 below 10~^ ! With this level of sensitivity 
oscillations driven by the sub-leading 5m^ scale might be observed [51]. 
For example, the number of muon decays required to produce 10 z/g events 
in a 50 kt detector 2800 km downstream of a neutrino factory is shown for a 
bimaximal mixing scenario (sin^ 2^13 = 0) in Fig. 29 as a function of the stored 
muon energy. Approaching 10^^ muon decays might be sufficient to observe 
oscillations driven by the sub-leading scale, but would require background levels 
of the order of 10~^ of the total CC rate, or better. 

With a vanishing or very small sin^ 2^^i3 only the z/^ — > z/^ oscillations will 
have a significant rate, and the oscillation parameters sin^ 2623 and 6171^2 can be 
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Figure 39: Allowed regions in oscillation parameter space calculated for a simu- 
lated experiment in which there are decays followed by N fi~ decays in a 
30 GeV neutrino factory that is 7400 km from a 10 kt ICANOE-type detector. 
The contours correspond to A^ = 10^° and 10^^ with and without backgrounds 
included in the calculation. Results are from Ref. 50. 

determined by fitting the right-sign muon {u^^ disappearance) spectrum. Good 
sensitivity can be obtained provided the baseline is chosen such that the first 
oscillation maximum occurs in the middle of the visible energy spectrum. As a 
first example, spectra of events tagged by right-sign muons are shown in Fig. 41 
as a function of 6ml2 for a 10 kt MINOS-type detector 2800 km downstream of 
a 30 GeV neutrino factory in which there are 2 x 10^°/i~ decays in the beam- 
forming straight section [51]. The position of the oscillation maximum (resulting 
in a dip in the observed distributions) is clearly sensitive to 5m|2. The depth of 
the observed dip is sensitive to the oscillation amplitude, and hence to sin^ 26*23. 
The visible energy spectrum of right-sign muon events can be fit to obtain 
sin^ 26*23 and We begin by considering the statistical precision that could 

be obtained with a perfect detector having MINOS-type resolution functions, no 
backgrounds, no selection requirements, and no systematic uncertainty on the 
neutrino flux. Fit results are shown in Fig. 42. For 6171^2 = 3.5 x 10^^ eV^/c^ the 
fit yields statistical precisions of a few percent on the the values of the oscillation 
parameters. If L is increased to 7332 km, the statistical precision improves to 
about 1%. With this level of precision it is likely that systematic uncertainties 
will be significant [60]. To illustrate this in Fig. 43a the la contours are shown 
in the ((5^32, sin^ 26*23) from fits which include backgrounds together with 0% 
and 2% systematic uncertainties on the beam flux. With a 2% flux uncertainty 
the precision on 6m'^2 sin^ 26*23 are respectively 11% and 14%. 
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Figure 40: Allowed regions in oscillation parameter space calculated for a sim- 
ulated experiment in which a 40 kt Fe-scintillator detector is a distance L km 
downstream of a 50 GeV neutrino factory in which there are 10^^// decays. The 
curves are 90% CL contours for L — 732 km (dashed), 3500 km (sohd), and 
7332 km (dotted). Results are from Ref. 59. 

As a second example, consider a 10 kt ICANOE-type detector that is down- 
stream of a 30 GeV neutrino factory in which there are lO^^/x"*" decays in the 
beam-forming straight section followed by 10^° /x~ decays [50]. The sensitivity 
to the oscillation parameters has been studied by fitting simulated visible energy 
distributions for events tagged by a right-sign muon. The analysis includes a 2% 
bin-to-bin uncorrelated systematic error on the number of neutrino interactions 
which takes into account the uncorrelated uncertainties on neutrino flux, the 
cross section, and the selection efficiency. To reduce background from charged 
meson decays, the events entering the fit are those with muons having momenta 
> 2 GeV/c. Figures 44-47 show fit results in the (sin^26'23, 5m|2)-plane as a 
function of the oscillation parameters and baseline. Note that for the "short" 
baseline (L = 732 km) the first oscillation maximum for the reference vahic of 
5m\2 = 3.5 X 10~^ eV^/c^ occurs at a neutrino energy of about 2 GeV. This is 
too low to produce a clear dip in the visible energy spectrum, and as a result 
sin^ 2^23 and 5m|2 can only be determined with relatively low precision (Fig. 44), 
and the fit results are sensitive to systematic uncertainties on the neutrino flux. 
At the longer baselines (L = 2900 km and 7400 km) the oscillation dip is visible, 
and the oscillation parameters can be measured with a precision that is mostly 
determined by the statistical uncertainty (Fig. 45). For a 30 GeV neutrino fac- 
tory and Sml2 = 3.5 x 10~^ eV^/c^ the longer baseline (7400 km) yields the 
most precise result. Specifically, for 10^°/x decays the statistical precisions on 
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Figure 41: Visible energy distributions for events tagged by a right-sign muon 
in a MINOS-type detector 2800 km downstream of a 20 GeV neutrino factory 
in which there are 2 x 10^°/i~ decays. Predicted distributions are shown for four 
values of with the other parameters corresponding to the LMA scenario 

lAl. For each panel, the points with statistical error bars show an example 
of a simulated experiment. The light shaded histograms show the predicted 
distributions in the absence of oscillations. Results are from Ref. 51. 
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Figure 42: Fit results for simulated disappearance measurements with a 10 kt 
MINOS-type detector 2800 km downstream of a 30 GeV neutrino factory in 
which there are 2 x lO^'^/i" decays. For each trial point the la, 2a, and 3a 
contours are shown for a perfect detector (no backgrounds) and no systematic 
uncertainty on the beam flux. The 68%, 90% and 95% SuperK regions are 
indicated. Results are from Ref. 51. 
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Figure 43: Fit results (1 a contours) for (a) simulated z/^ disappearance mea- 
surements with a 10 kt MINOS-type detector 2800 km downstream of a 30 GeV 
neutrino factory in which there are 2 x 10^°/x~ decays, with and without an ad 
hoc 2% systematic uncertainty on the neutrino flux, and (b) wrong-sign muon 
appearance measurements including an ad hoc 2% systematic uncertainty on the 
flux. The acceptance for a muon is zero for < A GeV and unity for > 4: 
GeV. Backgrounds are included but no p^ cut has been used. Results are from 
Ref. 60. 
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Figure 44: Fit results for simulated disappearance measurements with a 10 kt 
ICANOE type detector 732 km downstream of a 30 GeV neutrino factory in 
which there are lO^'^/x decays. The effect of a systematic uncertainty on the 
neutrino flux is shown. Results are from Ref. 50. 

sin^2^^23 and 5m\2 are respectively about 10% and 1%. With 10^"*^ /i decays the 
sin^ 26'23 precision improves by about a factor of 2. It should be noted that the 
best basehne choice depends on 5m\2 (Figs. 46-47), or more specifically 8771^2/ E. 

We conclude that, within the framework of three-flavor mixing, the oscilla- 
tion parameters sin^2^i3, sin^2^23, and 5m\2 can be determined at a neutrino 
factory by fitting the observed visible energy distributions for various event types. 
A comprehensive study of the expected precisions of the measurements as a func- 
tion of the oscillation parameters, baseline, and neutrino factory parameters has 
not yet been undertaken. However, detailed studies have been made for some ex- 
amples in which there are 10^°//+ decays followed by 10^° decays in a 30 GeV 
neutrino factory. For these examples we find that (i) if sin^ 2^13 > 0(10^^) global 
fits can be used to determine its value, (ii) if sin^ 2^^i3 is too small to observe 
~^ oscillations then we would expect to place the very stringent upper limit 
on its value of 10~^ or better, and (iii) the values of sin^ 2^23, and Sml2 could 
be determined with precisions of respectively better than or of order 10% and 
of order 1%, provided the baseline is chosen so that the dip corresponding to 
the first oscillation maximum is in the middle of the visible energy distribution. 
At a high-intensity neutrino factory (for example with 10^^ decays of 50 GeV 
muons) the mixing angles could be measured with a precision of a few percent, 
and if sin^ 2^13 is vanishingly small, the resulting upper limit could be at the 
O(10-^)-level. 
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Figure 45: Fit results for simulated disappearance measurements with a 10 kt 
ICANOE type detector 2900 km (top plot) and 7400 km (bottom plot) down- 
stream of a 30 GeV neutrino factory in which there are (a) 10^°/i decays and (b) 
10^^/i decays. Results are from Ref. 50. 

3.5.5 Search for CP violation 

In the majority of the three-flavor oscillation scenarios described in section 3.1 
the CP violating amplitude is too small to produce an observable effect. Nev- 
ertheless, in these cases stringent limits on CP violation would provide an im- 
portant check of the overall interpretation of the oscillation data. If however 
the LMA scenario provides the correct description of neutrino oscillations, CP 
violating effects might be sufficiently large to be observable at a high-intensity 
neutrino factory [59, 51]. This is illustrated in Fig. 48 which shows, as a function 
of baseline at a 20 GeV neutrino factory, the ratio R for 5 = and ±7r/2, where 
R is defined as the Vf. event rate divided by the z/g event rate. The 

upper group of curves is for 5m\2 < 0, the lower group is for 5m\2 > 0, and 
the statistical errors correspond to 10^^ muon decays of each sign with a 50 kt 
detector. If L is a few thousand km a non-zero S can produce a modification to 
R that is sufficiently large to be measured ! 

Since the Ue —>■ oscillation rates are to a good approximation proportional 
to sin^ 26'i3, it is useful to define the sin^ 26'i3 reach as that value of sin^ 26i^ that 
will produce a 3(T change in the predicted ratio R when 5 is changed from to 
±7r/2. The sin^ 2^13 reach is shown as a function of baseline and stored muon 
energy in Fig. 49 for a 50 kt detector at a neutrino factory in which there are 
lO^^/i"*" decays followed by lO^^yU" decays. With an optimum baseline of about 
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Figure 46: Fit results for simulated disappearance measurements with a 10 kt 
ICANOE type detector 2900 km downstream of a 30 GeV neutrino factory in 
which there arc 10^°/i~ decays followed by 10^°/i"'' decays. Results are shown for 
3 values of 5m§2, and are from Ref. 50. 

3000 km (for 5m\2 = 3.5 x 10"^ eV^/c^) the sin^2^i3 reach is approximately 
10~^, an order of magnitude below the current experimental bound. Thus, in 
a LMA scenario, CP violation in the lepton sector might be measurable at a 
neutrino factory providing 0(10^^) muon decays. 

As an example, consider a 40 kt Fe-scintillator detector downstream of a 

50 GeV neutrino factory providing lO^^/i"*" decays followed by 10^^//~ decays 
in the beam-forming straight section [59]. The results of fits to the simulated 
wrong-sign muon event distributions, with 6 and sin^ 29is left as free parameters, 
are shown in Fig. 50 for various baselines, with the sub-leading scale Smli ~ 
1 X 10~^ eV^/c^. The analysis includes the detector resolutions, reasonable 
event selection criteria, and backgrounds. As might be expected from Fig. 48 
at L = 7332 km there is little sensitivity to S, and at the "short" baseline L = 
732 km the fit has difficulty untangling 6 from sin^ ^Oi^. However, at a baseline 
of L = 3500 km for the example shown 6 and sin^ 2^13 can be determined with 
precisions of respectively about 15° and a few percent. Note that a combination 
of baselines can yiekl a modest improvement in the precision of the measurement. 
The sensitivity to CP violation decreases with decreasing 617121- Figure 51 shows 
as a function of sin^ 2^^i3 the lowest value of 6ml-^ for which the maximal CP 
phase 5 = 7r/2 can be distinguished from a vanishing phase at L = 3500 km. This 
limiting t^rnl^ is below the current central value for the LMA parameter space 
suggested by solar neutrino deficit, and is about 2 x 10"^ eV^/c^, independent 



85 



<^ 0.009 



I 0.008 



= 30 GeV, L = 7400 km, 2 x 10 |x decays 



a) = 3,5 x 10"' eV" 

b) ^m\:> = 5x10'^ eV 

c) = 7x10'' eV^ 



0.007 



c) 



0.006 



0.005 



b) 



0.004 



a) 



0.003 



90% C.L 



0.002 



0.2 



0.4 



0.6 



0.8 



1 

^23 



Figure 47: Fit results for simulated disappearance measurements with a 10 kt 
ICANOE type detector 7400 km downstream of a 30 GeV neutrino factory in 
which there are 10^°//+ decays followed by lO^^/x"*" decays. Results are shown for 
3 values of 5m|2, and are from Ref. 50. 
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Figure 48: The ratio i? of Pe — ^ ^ix to Ve — ^ix event rates at a 20 GeV neutrino 
factory for 5 = and ±7r/2. The upper group of curves is for 5m|2 < 0, the lower 
group is for 5m\2 > 0. The statistical errors correspond to 10^^ muon decays 
of each sign and a 50 kt detector. The oscillation parameters correspond to the 
LMA scenario lAl. With no matter or CP effects R ~ 0.5 for all baselines. 
Results are from Ref. 51. 
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Figure 49: Reach in sin^ 29i3 that yields a 3a discrimination between (a) 5 = 
and 7r/2 with 6ml2 > 0, (b) 5 = and 7r/2 with (5m|2 < 0, (c) 5 = and — 7r/2 
with 5m32 > 0, and (d) 5 = and — 7r/2 with (5m|2 < 0. The discrimination is 
based on a comparison of wrong-sign muon CC event rates in a 50 kt detector 
when 10^^ positive and negative muons ahernately decay in the neutrino factory. 
The reach is shown versus basehne for four storage ring energies. The oscillation 
parameters correspond to the LMA scenario lAl. Results are from Ref. 51. 
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Figure 50: Fit results in the CP phase 6 versus ^13 plane for a LMA scenario 
with 6m2i = 1 X 10~^ eV^/c^. The 68.5, 90, and 99% CL contours are shown for 
a 40 kt detector a distance L km downstream of a 50 GeV neutrino factory in 
which there are lO^^/i^ and 10^^/i~ decays. Results are from Ref. 59. 
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Figure 51: The lowest value of Smli, shown as a function of ^13, for which the 
maximal CP phase 5 = 7r/2 can be distinguished from a vanishing phase in a 
LMA oscillation scenario. The curve corresponds to a 40 kt detector 3500 km 
downstream of a 50 GeV neutrino factory in which there are lO^^/i"*" and 10^^/i~ 
decays. Results are from Ref. 59. 
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of sin^ 2^13. 

Finally we note that the sensitivity of short and medium baseline experiments 
to CP violation in a three-active plus one sterile neutrino scenario has been 
considered in Ref. [64]. They concluded that a 1 kt detector and a 100 km 

baseline could provide a clean test of CP violation, particularly in the r lepton 

appearance channel. 

3.6 Summary 

The oscillation physics that could be pursued at a neutrino factory is compelling. 
In particular, experiments at a neutrino factory would be able to simultane- 
ously measure, or put stringent limits on, all of the appearance modes i/g — > i^r, 
J^e J^ix, and Ufj, — > Ur- Comparing the sum of the appearance modes with the 
disappearance measurements would provide a unique basic check of candidate 
oscillation scenarios that cannot be made with a conventional neutrino beam. 
In addition, for all of the specific oscillation scenarios we have studied, the z/g 
component in the beam can be exploited to enable crucial physics questions to 
be addressed. These include (i) the pattern of neutrino masses (sign of 5m^) and 
a quantitative test of the MSW effect, (ii) the precise determination of (or strin- 
gent limits on) all of the leading oscillation parameters, which in a three-flavor 
mixing scenario would be sin^ 26*13, sin^ '2023, and 5m|2, and (iii) the observation 
of, or stringent limits on, CP violation in the lepton sector. 

To be more quantitative in assessing the beam energy, intensity, and baseline 
required to accomplish a given set of physics goals it is necessary to consider 
two very different experimental possibilities: (a) the LSND oscillation results 
are not confirmed by the MiniBooNE experiment, or (b) the LSND results are 
confirmed. 

(a) LSND not confirmed. Fairly extensive neutrino factory studies have been 
made within the framework of three-flavor oscillation scenarios in which 
there is one "large" dni^ scale identified with the atmospheric neutrino 
deficit results, and one small 6m'^ identified with the solar neutrino deficit 
results. A summary of the energy dependent beam intensities required to 
cross a variety of "thresholds of interest" is provided by Fig. 29. A 20 GeV 
neutrino factory providing 10^^ muon decays per year is a good candidate 
"entry-level" facility which would enable cither (i) the first observation of 
i^e i^ix oscillations, the first direct measurement of matter effects, and 
a determination of the sign of 5m\2 and hence the pattern of neutrino 
masses, or (ii) a very stringent limit on sin^ 2^13 and a first comparison of 
the sum of all appearance modes with the disappearance measurements. 
The optimum baselines for this entry-level physics program appears to be 
of the order of 3000 km or greater, for which matter effects are substan- 
tial. Longer baselines also favor the precise determination of sin^2^i3. A 
20 GeV neutrino factory providing 10^° muon decays per year is a good 
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candidate upgraded neutrino factory (or alternatively a higher energy fa- 
cility providing a few xlO^^ decays per year). This would enable the first 
observation of, or meaningful limits on, i/g — > Vj- oscillations, and precision 
measurements of the leading oscillation parameters. In the more distant 
future, a candidate for a second (third ?) generation neutrino factory 
might be a facility that provides 0(10^^) decays per year and enables the 
measurement of, or stringent limits on, CP violation in the lepton sector. 

(b) LSND confirmed. Less extensive studies have been made for the class of 
scenarios that become of interest if the LSND oscillation results are con- 
firmed. However, in the scenarios we have looked at (IBl and ICl) we 
find that the Vf, — > Vr rate is sensitive to the oscillation parameters and 
can be substantial. With a large leading Srn^ scale medium basehnes (for 
example a few xlO km) are of interest, and the neutrino factory inten- 
sity required to effectively exploit the v^. beam component might be quite 
modest (< 10^^ decays per year). 

The neutrino factory oscillation physics study we have pursued goes beyond 
previous studies. In particular we have explored the physics capabilities as a 
function of the muon beam energy and intensity, and the baseline. Based on the 
representative oscillation scenarios and parameter sets defined for the study, it 
would appear that a 20 GeV neutrino factory providing 0(10^^) decays per year 
would be a viable entry-level facility for experiments at baselines of ~ 3000 km or 
greater. There are still some basic open questions that deserve further study: (1) 
We have sampled, but not fully explored, the beam energy and intensity required 
to explore the scenarios that become relevant if the LSND oscillation results are 
confirmed. (2) Possible technologies for a very massive neutrino factory detector 
have been considered, but these considerations deserve to be pursued further. 
The chosen detector technology will determine whether it is necessary to go deep 
underground. (3) We have developed tools that can explore the utility of having 
polarized muon beams. The physics payoff with polarization is a detailed issue. 
It deserves to be studied in the coming months. 

Based on our study, we believe that a neutrino factory in 5-10 years from 
now would be the right tool at the right time. 

4 Non— Oscillation Physics 

Due to the theoretically clean nature of weak interactions, conventional neutrino 
scattering experiments have always provided precise measurements of fundamen- 
tal parameters. These include: a crucial role in the extraction of parton distribu- 
tion functions, measurements of the Weinberg angle[65], and the strong coupling 
constant [66] a^, which are competitive with any other methods. Perhaps be- 
cause of this success, we forget how crude existing neutrino experiments are. 
The high statistics experiments such as CDHSW[67] and CCFR/NuTeV[66, 65], 
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in order to obtain samples of more than 10^ events, rely on coarsely segmented 
massive iron/scintillator calorimeters weighing close to 1000 tons. Measurements 
on proton targets and detailed studies of the final state have been confined to 
very low statistics bubble chamber experiments. As a result we have virtually 
no precise measurements of neutrino-proton scattering and no measurements on 
polarized targets which could offer new insights into the spin structure of the 
nucleon. 

The advent of a neutrino factory, with neutrino fluxes of 10^°/year instead 
of the 10^^~^^ at existing facilities would open a new era in conventional neu- 
trino physics. We would be able to use low mass targets and high resolution 
detection technologies and still achieve better statistical power than present-day 
experiments. For example a 50 GeV muon storage ring at the above rate would 
produce around 18 M neutrino charge-current interactions per year in a 10 kg 
hydrogen target. This is 5-10 times the statistics of the CCFR and NuTeV 
experiments with 600 ton detectors. Better understanding of neutrino fluxes 
from the decay of monochromatic muons will also reduce many of the dominant 
systematic errors. 

In this study we have concentrated on measurements that are only possible 
with higher fluxes rather than repeating older measurements with thousands 
of times the statistics. As a result, the statistical errors shown are often not 
negligible, but without the high flux at a neutrino factory the measurements 
themselves would be impossible. 

Outline 

Due to the breadth of the field we are unable to give a complete survey and 
instead highlight a few of the areas where the high flux beam at a neutrino 
factory allows new measurements: 

• A description of a low mass target/detector and typical rates in such a 
detector. 

• Nucleon deep inelastic scattering measurements and a proposed detector 
design. 

• Neutrino cross section measurements, a topic of great interest to the nuclear 
physics community and also needed to understand normalization at a far 
neutrino oscillation detector. 

• Spin structure functions, which have never been measured in neutrino 
beams. 

• The potential of the neutrino factory as a clean source of single tagged 
charm mesons and baryons. 
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Figure 52: Kinematics of neutrino scattering in the parton model. The energy- 
momentum tranfer from the leptons to the proton is and the fraction of the 
proton momentum carried by the struck quark is approximately x. 

• Elcctroweak measurements in both the hadronic and purely leptonic sec- 
tors. 

• Use of the very clean initial state to search for exotic interactions. 

• Searches for anomalous neutrino interactions with electromagnetic fields. 

4.1 Possible detector configurations and statistics 

For studies of charged current deep-inelastic scattering on proton targets, the 
optimal detector system is probably a target followed by precision magnetic 
tracking sytems, an electromagnetic calorimeter and a muon detection system. 
Such detectors have been used in muon scattering experiments at CERN and 
FNAL and in the new generation neutrino scattering experiments CHORUS [68] 
and NOMAD [58]. A low mass target followed by tracking and electromagnetic 
calorimetry makes the electron anti-neutrinos in the beam a source of additional 
statistics rather than backround, except in the case of neutral current studies. 

The target itself should be thin enough that particles produced within it have 
a small probability of interacting before they reach the tracking systems. In this 
study we considered liquid hydrogen and deuterium targets - both polarized and 
unpolarizcd - and heavier solid nuclear targets. The hydrogen and deuterium 
targets are Im long while the polarized target is 50 cm long. All targets arc 20 
cm in radius, to fit the central beam spot at 50 GeV. For lower beam energies 
the beam spot grows in size as ~ 1/E. Nuclear targets are scaled so that the 
interaction length in the material is constant at 14%. The charged current muon 
neutrino interaction rates are summarized in table 10. 

The numerical estimates in this stTidy use, unless otherwise noted, 10^° 50 
GeV muon decays in a 600 m straight section. 
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Table 10: Charged current muon-neutrino scattering rates in a small target 
located near a muon storage ring. Rates are per 10^° muon decays. The detector 
is located {Ix E^, GeV)meters away from the ring to assure that primary muons 
have ranged out before the detector. 



Machine 


Target 


Thickness, cm 


Events 


50 GeV neutrino factory 


Liquid H2 


100 


12. IM 




Liquid D2 


100 


29.0M 




sohd HD 


50 


9.3M 




C 


5.3 


20.7M 




Si 


6.3 


25.4M 




Fe 


2.3 


31.6M 




Sn 


3.1 


39. IM 




W 


1.3 


44.3M 




Pb 


2.4 


46.5M 


CCFR/NuTeV 


Fe 


600 


~ 2M 



These are the total event rates for charged current muon-neutrino scattering. 
The anti-neutrino rates arc half as large. Kinematic cuts reduce the statistics by 
less than a factor of two. We have only considered muon-neutrino charge current 
scattering for structure function measurements, although for such thin targets, 
electron neutrino scatters should also be reconstructable with high precision. 



4.2 Neutrino Scattering Kinematics 

The kinematic variables for neutrino deep inelastic scattering are illustrated in 
figure 52: 

q = - kf , = -q^ = 2EeE^ - - 2E^pe cos Oiab, (58) 

u = {ppq)/Mc^Ee-E'„ (59) 

X = Q^/2Mu, (60) 

y = Mu/iKPp) = (1 + cosecM)/^ ~ u/Ee, (61) 

= 2Mu + M^ -Q^, (62) 

where the k are the neutrino and final state four vectors, Pp is the proton four- 
vector, M is the target nucleon mass, E^ is the incoming neutrino energy Ei,pi 
are the outgoing Icpton energy and momentum 6iab is the lepton angle with 
respect to the incoming beam, q is the four-momentum transfer to the target, u 
is the energy transfer, x is the Bjorken x variable, y is the scaled energy transfer 
and is the invariant mass of the final state hadronic system squared. 

Fig. 53 shows the kinematic region for a neutrino factory as compared to 
other deep-inelastic scattering experiments. 
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Figure 53: Comparison of kinematic ranges for present DIS experiments with a 
50 GeV Neutrino factory. 
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For Q « E and s « Mw the the unpolarized neutrino (anti-neutrino) 
scattering cross section is: 



'dxdy " — 2^ — ^[F^' '{x,Q )±xF^' '[x,Q )]+ (63) 

[F,^^^\x,Q')TxF^^''\x,Q')]{l-yr- 
2y'FL{x,Q'), 

where the Fi are Structure Functions. Fl = F2 — 2xFi is a purely longitudinal 
structure function. The xF^ contribution changes sign for anti-neutrino scatter- 
ing. There are additional structure functions F4 and F5 which are suppressed by 

factors of the Icpton mass over the proton mass squared. For z/^ and scattering 
at very low energies, these terms can become quite important. 

If the target is longitudinally polarized with respect to the neutrino polar- 
ization, then the cross section difference [80]: 

d^y = - 2 - W^"^^ ^ -^^^ + 

9 ,^ xM . . xyM..,^ xM . 

y^x{l + -^)^3 + (1 - y - ^) [(1 + ^)^4 + ^5]}, 

is described by two parity conserving Polarized Structure Functions gi and g2, 
and by three parity violating Polarized Structure Functions g^.g^ and g^y. How- 
ever, if the nucleon is transversely polarized, the cross section difference is: 



^'(^r ^r)'^^'^^ _ GlM ^ l^yM[2{l-y)E-xyM]{ t 2xy{^g, + g^) (65) 



dxdy 167r^ ^ 2' 

2 ,^ xyM y 

+xy g3 + {l-y- -^)94 - ^95}- 

The transverse cross section is suppressed by M/ Q with respect to the longitu- 
dinal cross section. 



4.3 Total cross section Measurements 

A measurement of the total CC neutrino scattering cross section is both of 
intrinsic interest and essential to precision measurements at a neutrino factory. 

Wc currently know the cross sections for neutrino scattering at the 2-3% level 
[14] at energies above 30 GeV but at energies approaching the resonance region 
(2 GeV) the uncertainty increases considerably. Because muon decay is so well- 
understood, the flux and hence the total cross section should be measureable 
across the full energy spectrum to the 1% level. 
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The yield Y of neutrino interactions observed in any detector can be written 

as: 

Y = no{E,r)a{E)e{E,r)N 

where no is the flux of incident neutrinos as a funtion of the neutrino energy E 
and distance from the center of the detector, r. a is the cross section, e{E, r) is 
the detector efficiency, and N is the number of target particles. 

For a far detector, nQ{E,r) ~ no{E) is mainly determined by the beam 
divergence and the muon decay kinematics and can probably be estimated from 
the machine parameters and decay model with precisions at the 1% level. For 
a near detector, no{E,r) depends mainly on the muon decay kinematics and 
geometry, with contributions from beam size and divergence at the few percent 
level. 

Without a measurement of the absolute number of neutrinos, the best way 
to determine the flux is to normalize to a very well-understood process. By 
comparing the yield of the normalization process and the total event rate, one 
then has a measurement of the total cross section. Inverse muon decay (z/^+e~ 

^ l-Ue and Ve + — > fi hX^/^) provides a clean channel for mapping the 

beam flux no{E, r) at a near detector and for normalizing the total cross section 
measurement. It has the limitation of an energy threshold of ~ 11 GeV and 
no corresponding channels for the opposite sign beam. Quasi-elastic scattering 
is an additional normalization mode since it has a much lower energy threshold 
and occurs for beams made with muons of either charge. Finally, scattering from 
atomic electrons is suppressed by a factor of order me/rrip relative to the normal 
neutrino nucleon interactions, but still yields an event rate of ~ 10^ interactions 
per gr/cm^ for 10^° 50 GeV jjT decays. 

If the ratio of flux shapes at far and near detectors can be understood at the 
1% level, then measurements of nQ{E,r)a{E) at a near detector can be used to 
precisely predict the number of events expected in the absence of oscillations at 
a far detector. Such precise flux measurements are also important for the suite 
of measurements described in the remainder of this chapter. 

4.4 Structure function measurements 

In principle, the structure functions can be extracted by fits to the y dependence 
of the cross section. To date this has proven very difficult as the data must be 
binned in x, y and and no experiment has had sufficient statistics to perform 
such an analysis with high accuracy[67, 70]. 

Instead, high statistics experiments[14] such as CHARMII, CCFR and CDHSW 
have relied on massive targets (Iron, Calcium) which are approximately iso-scalar 
and have combined neutrino and anti-neutrino information in order to extract 
average structure functions. The structure functions averages have naive parton 
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model interpretations: 

F^{x,Q') ~ J2i^q{x,Q')+xq{x,Q')), (66) 
F^ix,Q') ~ Y.i^q{x,Q')-xq{x,Q')), (67) 

where F2{x,Q^) and F^i^x^Q"^) arc the average of neutrino and antineutrino 
structure functions measured on a target whicli is an average of neutron and 
proton and q(a;, Q"^) and q(a;, Q"^) represent the parton distribution functions or 
total probability of finding a quark or antiquark in the proton: 

q(x,g2) = u(a;,g2) +d(a;,g2) +s(x,g2) +c(a;,g2)... (68) 
q(a;,g2) = u(x, g^) + d(a;, g^) + s(x, g^) + c(a:, g^)... (69) 



Given the expectation of 12 M (24 M) events/year in a 1 m hydrogen (deu- 
terium) target at a 50 GeV muon storage ring we can do a complete analysis of 
each channel i/p, un^ Vp, Vn without averaging. Such an analysis allows a unique 
extraction of individual quark flavor parton distribution functions. 

For example, in the case of i/p scattering, a is exchanged and the reaction 
is only sensitive to negatively charged quarks. Due to the helicity dependence 
of the interaction only left-handed d type and right handed u quarks will be 
involved. 

The leading order parton model cross section is simply 

^ ~ M^a;[(d4a;,g2)+s^(a-,g2)) + (u«(x,g2)-t- (70) 

cn{x,Q')){l-y)\ (71) 



and the different contributions can be extracted from the y dependence of this 
cross section and the corresponding anti-neutrino cross section. The relative s 
and d quark contributions can be measured in charm production. 

For an unpolarized target qL{,x) = qR^x) = \q{x). For a polarized quark 
<iL{x) = \{q{x) + Sq{x)) and qR{x) = ^{q{x) — Sq{x)) where 5q{x) is the degree 
to which the quark spin is aligned with the proton spin^. Thus a a^p measurement 
on an unpolarized target can determine d -|- s and u -|- c by averaging over the 
proton spin, while by measuring the polarization asymmetry one can measure 
Sd + Ss and Su + Sc. 

Scattering on neutrons can be related to scattering on protons by an isospin 
transformation which exchanges u and d quarks and anti-quarks. Differences 

^The traditional Aq spin distributions from electron and ninon scattering measure the sum 
Aq = Sq + 6q as photon probes cannot tell quarks and anti-quarks apart. 
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of neutron and proton cross sections can then be used to cancel the u and d 
components leaving observables sensitive only to s and c distributions. 

4.5 Perturbative QCD 

Neutrinos do not couple directly to gluons. As a result, QCD effects appear 
in neutrino scattering as higher order corrections to the leading order parton 
model. Measurements of the dependence of neutrino cross sections are some 
of the most sensitive measurements of the strong coupling constant as[66] and 
information on the gluon distribution can be obtained from its coupling to the 
structure functions via the DGLAP[69] evolution equations. The neutrino struc- 
ture functions can be divided into two types, singlet and non-singlet, depending 
on their sensitivity to gluon effects in their evolution. 

The structure functions 2xFi, F2 and gi are singlet functions and are directly 
coupled to the gluon distribution via the evolution equations. The structure 
functions xF^ + xF^, 2xg3, and g^ averaged over neutrino and anti-neutrino are 
non-singlet functions and their evolution is independent of the gluon distribution. 
The combination — also cancels the gluon contributions and is thus non- 
singlet in nature. 

To date, extractions of as from non-singlet distributions have been statistics 
limited and strongly affected by flux uncertainties. The additional factor of 10- 
100 in statistics and improved flux understanding at a neutrino factory should 
allow vastly improved measurements of strong interaction parameters in this 
very clean channel. 

Once the quark distributions and strong interaction effects have been thor- 
oughly studied in the non-singlet structure function, that knowledge can be used 
for improved constraints on the gluon distributions via the evolution of the sin- 
glet structure functions. 

4.6 Nuclear Effects 

Experiments at a neutrino factory of nuclear effects in the distribution of par- 
tons within nuclei relative to protons and deuterons are of interest to both the 
nuclear and high energy communities. These nuclear effects have been stud- 
ied extensively using muon and electron beams but have only been observed in 
low-statistics bubble chamber experiments[72] using neutrinos. If we consider 
the behavior of the structure functions F2{x, Q^) measured on a nucleus (A) to 
F2{x, Q"^) measured on a nucleon as a function of x we pass through four distinct 
regions in going from x — Q io x — l.Q: 



98 




o.oo O.Ol 0.02 0.03 0.04 0.05 



Figure 54: Rca:D2 fo^' both F2 and xF^ as measured with 14 M events on each 
target. 



Shadowing Region x < 0.1 

In the shadowing region {x < 0.1) there are several effects that should yield 
a different ratio Ra = -^2{A)/-^2(Af) when using neutrinos as the probe. In the 
limit — > 0, the vector current is conserved and goes to 0. The axial- vector 
part of the weak current is only partially conserved (PCAC) and F2{x,Q'^) — » 
a non-zero constant as — > 0. According to the Adler theorem [73] the cross 
section of — N can be related to the cross section for vr — at = 0. This 
relation can be studied in both proton and in heavy nucleii. 

As we increase from but keep it under 10 GeV^ in the shadowing region 
we enter the region of vector meson dominance (VMD) in fi/e — A scattering. 
The physics concept of VMD is the dissociation of the virtual boson into a 
quark/antiquark pair, one of which interacts strongly with the 'surface' nucleons 
of the target nucleus (thus the 'surface' nucleons 'shadow' interior nucleons). In 
u — A scattering there is an additional contribution from axial- vector mesons that 
is not present in fj,/e — A scattering. Boros et al. [74] predict that the resulting 
shadowing effects in u—A scattering will be roughly 1/2 that measured in fi/e—A 
scattering. In a more quantitative analysis, Kulagin [75] used a non-perturbative 
parton model to predict shadowing effects in u — A scattering. At 5 GeV^ he 
predicts equal or slightly more shadowing in u — A scattering than in /i/e - A 
scattering. He also attempts to determine quark flavor dependence of shadowing 
effects by separately predicting the shadowing observed in F2{x, Q"^) (sum of all 
quarks) and xF^^x, Q^) (valence quarks only). Fig. 54 shows the results of a run 
with 14M events/target using predictions of Kulagin's model for F2 and xF^. As 
can be seen, the predicted difference between the shadowing on sea and valence 
quarks is clearly visible down to x ~ 0.03. 
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Anti-shadowing Region (0.1 < x < 0.2) 

Drell-Yan experiments have also measured nuclear effects. Their results are 
quite similar to DIS experiments in the shadowing region. However, in the anti- 
shadowing region, where makes a brief but statistically significant excursion 
above 1.0 in DIS, Drell-Yan experiments see no effect. This could be an indication 
of difference in nuclear effects between valence and sea quarks. Eskola ct al [76] 
quantified this difference by using a leading order/leading twist DGLAP model. 
Taking the work of Kulagin and Eskola together implies that nuclear effects 
in xF-i{x,Q'^) should be quite dramatic with more shadowing than F2{x,Q'^) 
at lower x and then i?^ rising fairly rapidly to yield significant antishadowing 
around a; = 0.1. With 14 M events on each target we should be able to measure 
antishadowing effects and the difference between shadowing effects in F2{x, Q^) 
and xF3(x, Q^) to the 6 a statistical level. 

EMC-efFect Region(0.2 <x< 0.7) 

Determination of individual quark contributions to the EMC-effect will be chal- 
lenging since the participation of sea quarks, and thus the difference between 
F2{x,Q'^) and xFs{x,Q'^), shrinks rapidly with increasing x. However, Eskola's 

— (2) 

predictions for this region indicate that the contribution of u and d to R\ in 
the range of this experiment remains well below unity so that the quantity 
{Rj^ - Ra ) should remain negative well into the EMC-effect region. 

Behavior of F2{x,Q'^) as a; 1 in a Nuclear Environment 

When working in the fermi-motion region it has been shown that we need to add 
more than the Fermi gas model to a simple nucleon to reproduce the behavior 
of F2{x,Q^) at high x. Few-nuclcon-correlation models and multi-quark cluster 
models allow quarks to have a higher momentum which translates into a high-x 
tail. In this region F2{x, Q"^) should behave as e~"^. There have been analyses of 
this behavior in similar kinematic domains using fi+C and z/+ Fe interactions. 
The BCDMS [77] muon experiment finds a = 16.5 ± 0.5 while the CCFR[78] 
neutrino experiment finds a = 8.3 ± 0.7 ± 0.7 (systematic). Is the value of a 
dependent on the nucleus? One would expect any few nucleon correlation or 
multi-quark effects to have already saturated by Carbon. Is a dependent on the 
probe? 

Summciry 

There is a rich program of studying nuclear effects with a neutrino probe in a 
high statistics neutrino factory experiment. The effects could be measured to 
statistically significant accuracy in a 2 year exposure to the beam in the near- 
detector experiment described above. The data gathered would allow separate 
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measurements of the effects on valence quarks and sea quarks across mucli of tfie 
X range. 

The nuclear community would surely be excited by this valuable tool for 
nuclear research at a neutrino factory. 

4.7 Spin Structure 

An intense neutrino beam at a neutrino factory would create significant event 
rates in compact detectors. This opens the possibility of using a polarized target, 
and hence a completely new class of neutrino measurements becomes possible. At 
present we know very little about the spin structure functions Qi — introduced 
in Equations 65 and 66. In particular, the parity violating functions have only 
been explored via weak-interference measurements of proton form factors by 
the SAMPLE collaboration [79] with low statistics. A neutrino factory would 
allow direct high-statistics measurements of all of these structure functions and 
should be able to answer many unresolved questions about the spin structure of 
the nucleon. 

Formalism 

The nucleon spin (|) can decomposed in terms of quark and gluon contributions: 



where AE = Ait -|- Ad -I- As -|- Ac is the net quark hehcity and Ag is the net 

gluon helicity along the nucleon spin direction, while Li are their relative orbital 
angular momentum. ( We use Aq as a shorthand for the integral / Aq{x)dx.) 

To date, the only experiments which have studied the spin structure of the 
nucleon are low energy charged lepton polarized deep-inelastic scattering exper- 
iments (PDIS) where only the parity conserving polarized structure functions g{ 
and g2 can be measured. 

g{ can be written in the leading order parton model as a sum of a nonsinglet 
and singlet part [81]: 



The first non-singlet term evolves independently of the gluonic spin contri- 
bution while the second is coupled to, and thus depends on, the gluon spin 
contribution Ag. 




(72) 
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The integral structure functions have the following relation to the parton 
spin contributions: 



j dxg[{x,Q^) 

KNsiQ') + KsiQ' 



1 1 

^as + -as 



Cf 



1 1 

-aa H — as 

2 6 



(75) 
(76) 

(77) 



Where the Ci are coefficient functions and the axial charge matrix elements 



03 = F + D~Au-Ad (78) 
as = 3F-D~ Au + Ad-2As (79) 
ao = Au + Ad + As = AE (80) 

(81) 

can be expressed in terms of coupling constants F and D obtained from neutron 
and hyperon beta decays [83]. Because the interaction between A^f and AE in 
the evolution of the singlet (oq) component, interpretation of in terms of the 
quark spin is problematic. Fig. 55 shows NLO QCD predictions for AE as a 
function of /\g. The data in the NLO fit are from [82]. 

Neutrino beams introduce both additional parity violating spin structure 
functions QziQi and and new combinations based on sums and differences of 
neutrino and anti-neutrino scattering. 

For example [80] the sums 



rr + rr = j dx[g7{x,Q'')+g7{x,Q')] (82) 
rr + Ff = j dx[gT{x,Q^)+gl^{x,Q^)] (83) 

for both proton and neutron targets are only sensitive to the singlet ao term 
and no input from beta decay is necessary. 

The parton model interpretation of these new structure functions is: 

OZM') = 2xgrix,Q') (84) 
= -x[6d{x, Q^) + 5s{x, Q2) - du{x, Q^) - 5c{x, Q^)] , 

glU^^Q") = '2^97{^,Q') _ (85) 

= -x[6u{x, g2) + 6c{x, g2) - 6d{x, Q^) - 5s{x, Q'^)] . 
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Figure 55: Model dependent decomposition of singlet term into quarks and gluon 
based on PDIS data, qq ^1 — S^A^f. The QPM expectation and the results 
from a NLO fit of the evolution of most of the available data on g\ are also 
shown. From the fit it was found that Ag = 1.6 ± 0.3 ± 1.0, where the error is 
dominated by theoretical uncertainties. 



103 



On a deuterium target, the u and d contributions to can be cancelled 
leading to a direct measurement of the strange sea contribution to the nucleon 
spin [80] 

^uinp) _ ginp) ^ _2(^g ^ ^ 2((5c + 5c), (86) 

which can also be studied via polarization asymmetries in charm production 
from strange quarks [84]. 

The structure functions xq^iQa and g^, like F3 are non-singlet functions in 
which contribution from gluons cancel. Comparison of the non-singlet functions 
with the single functions gi and F2 is an indirect way of measuring the contri- 
bution of gluons A^i. 

Experimental Setup 

A promising target technology is the 'ICE' target [85] , a sohd hydrogen-deuterium 
compound in which the protons or the deuterons can be polarized independently. 
The expected polarization and dilution are Pf/=80% and /// = 1/3 for H, and 
Pd=50% and fo = 2/3 for deuteron. A 7 kg (pt=l.l gr/cm^) polarized target 
with the qualities mentioned above would be 20 cm in radius and 50 cm thick, 
similar to the other light targets proposed for structure function studies. Raw 
event rates in such a detector would be around 20M per 10^*^ muon decays. 

If such a data sample is analyzed in 10 in x bins, the error in each x bin 
would be 5gi ~ (fPrVN)-^ ~ 1%. 

If the neutrino beam intensities and polarized target described above are fea- 
sible, the physics motivations would be very strong. We will be able to do high 
precision measurements where we can cleanly separate singlet (gluon-dependent) 
from non-singlet (gluon-free) terms. Furthermore, due to the nature of the neu- 
trino charged current interactions it will be possibility to perform a measurement 
of the polarization of the proton's quarks by flavor, with sea and valence contri- 
butions separated. 

4.8 Charm Production and — Mixing 

Neutrino interactions are a very good source of clean, sign-tagged charm parti- 
cles. Single charm quarks are produced via the processes 



vs —>■ i^c Cabbibo favored (87) 

ud i'c Cabbibo suppressed (88) 

us — > i'^c Cabbibo favored (89) 

ud — > i'^c Cabbibo suppressed (90) 
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Figure 56: Charm and bottom quark production as a fraction of the total cross- 
section as a function of E,^. 
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E,, =50 GeV 
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Figure 57: Charmed hadron spectra from neutrino interactions in a near detector 
from a 50 GeV muon storage ring. 

The fraction of heavy flavor produced as a function of E^, is shown in Fig. 56. 
An experiment at a 50 GeV muon storage ring with 10^° muon decays and a a one 
ton (fiducial) target made up of silicon strip detectors interleaved with heavier 
material would observe ^ 3 x 10^ muon-neutrino charged- current interactions 
and around 1.2 x 10* charm hadrons with energies above 4 GeV/year. All of 
these charmed hadrons are flavor tagged at the point of production by the charge 
of the outgoing primary lepton (c production with i~ and c production with 

There are several interesting physics motivations for charm studies at muon 
storage rings, including measurements of the strange contribution to proton 
structure and spin; however, the primary motivation for producing charm by 
this method is the cleanliness of the final state relative to hadroproduction and 
the flavor tagging in production. This experimental fact compliments the the- 
oretically "clean laboratory" of charm in searches for FCNC, CP asymmetries 
and D'^ oscillations, all of which are small in the standard model because 

of the lack of coupling of charm to the heavy top quark. 

Although this study has concentrated on neutrino energies below 50 GeV, we 
note that similar arguments hold for bottom production and that for machines 
with energies above 100 GeV, single B meson production rates can reach 100 
per gr/cm^ of target. Because nuclei mainly consist of u quarks rather than c 
quarks, the u — > b rate will be enhanced and a clean measurement of Vub without 
final state effects may be possible. 

As an example of the physics reach of a neutrino charm factory, consider the 
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example of — mixing measurements. The most sensitive current searches 
for time- integrated mixing place limits on the process of ~ 5 x 10~^ [90, 91]. 
BaBar expects to have sensitivity to mixing at the ~ 5 x 10~^ level after several 
years at design luminosity [92]. These measurements are ultimately limited by 
tagging mistakes and backgrounds to final state or D° identification from 
doubly- Cabibbo suppressed decays, such as K^'k~ which occur at the few 

part per thousand level. 

At a 50 GeV muon storage ring, with a high mass detector, one could probe 
£)0 _ £)0 mixing via 

vN ^ cl-X 

£-X, 

and its charge conjugates. The appearance of like-signed leptons would indicate 
mixing, where opposite-signed leptons are expected. Assuming 50% of the charm 

produces hadronizes as a or D^, this would result in the observation of 2 x 10^ 
tagged neutral charm meson semi-leptonic decays in either the muon or electron 
channel. 

4.9 Precision Electroweak Measurements 

Precision measurements of electroweak parameters from neutrino experiments 

have played an important role in testing the Standard Model and in searching 
for new physics. Even with the wealth of on-shell W and Z bosons today at 
colliders, the neutrino data remains important and is complimentary to collider 
studies. The intense fiux of neutrinos from a neutrino factory opens up a new 
era of precision electroweak measurements previously limited by statistics. 

Three interesting precision measurements can be contemplated at a neutrino 
factory, and each can be cast as a precision measurement of the weak mixing 
angle, sin^ 9w- The experiment that most dramatically highlights new capabil- 
ities at a neutrino factory is the study of neutrino-electron cross sections. The 
second possible experiment is the extraction of sin^ 9iy from the neutral-current 
to charged-current deep inelastic scattering (DIS) cross-sections, a measurement 
which currently is the most precise test of weak interactions in neutrinos and 
which is currently limited by statistics [86], although difficult theoretical system- 
atics are not far behind. A final possibility is the study of weak boson scattering 
from photons in the so-called neutrino "trident" process. 

Only the first two possibihties are discussed here in detail. Trident processes 
are not considered here because of the difficult theoretical systematics that will 
ultimately limit the measurements [87]. 
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Neutrino-electron scattering 



Neutrino-electron elastic scattering, 

pe~ i/e~, (91) 

is perhaps the most promising reaction for the precise probe of electroweak un- 
ification from neutrino interactions. Because the target particle is point-like, 
its structure does not introduce uncertainties in extracting parameters of the 
electroweak interaction from observed cross-sections. These measurements will 
therefore be limited only by statistical and experimental uncertainties. The best 
previous measurements of i/e scattering come from the CHARM II experiment 
which observed approximately 5000 events [88] . 

Two distinct measurements are possible at a muon storage ring because 
neutrino-electron scattering includes different diagrams for the beams from pos- 
itive and negative muons. The processes in the i^ijyl and VJ^Ve beams are 

i/^e~ i/^e~. Pee" Pge", (92) 

where the electron-neutrino and electron-antineutrino scattering process includes 
a charged-current (VF^ exchange) in the t-channel and s-channel, respectively. 
The differential cross-section with respect to y = E^/ E^, is given by 



daive — >• ve ) 2G%meEi, 
dy TT 



gl + dU^-yf , (93) 



where terms of 0{me/ Ey) are neglected and expressions for the left-handed and 
right-handed coupling constants, ql and qr, for each of the processes shown in 
Eq. (92) are given in Table 11. The numerical values for the cross-sections after 
integrating over y are: 



a{ve' ve~) = 1.72 x 10"^' x Ey[GeV] x 



2 I -'- 2 
9l + g^R 



(94) 



where the values for the final term are also given in Table 1 1 . Note that neutrino- 
electron scattering has a much lower cross-section than DIS, roughly down by the 
ratio of mg to nip. Experimentally, it should be noted that the observed neutrino- 
electron scattering rate will be summed over both beams since the observed final 
states are identical. 

The experimental signature for z/ — e scattering is a single negatively charged 
electron with a small transverse momentum relative to the incoming neutrino, 

pf-""^ < pm,E,. (95) 

The normalization mode is the appearance of a single muon with similarly low 
. Of course, the neutrino beam itself has a characteristic divergence from 
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Table 11: ql and QRhy v — e scattering process 



Reaction 


9l 


9r 






— 1 + sin^ 611^ 


sin^ 611^ 


0.0925 


Vf,e~ —>■ t'eC" 


i + sin^ 6w 


sin^ 9\v 


0.5425 




sin^ 


-i + sin^ 


0.0758 


i/ge" — >• z/ge^ 


sin^ 9w 


1 + sin^ 


0.2258 



i/e Scattering Signal/Background at low p„ /j,'^ Beam, perfect resolution 
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Figure 58: Signal to noise in the low pt region (p| < + ""2^" ) as a function 
of E,. 



mean beam direction is given by 



decay kinematics of and therefore the observed lepton relative to the 



(e) 
Pt 



+ (96) 



For a 50 GeV storage ring, this factor is dominated by the fundamental pt of the 
interaction and is typically ~ 90 MeV. For a lower energy storage ring of about 
15 GeV, these factors become equal. 

The primary background to this measurement is from quasi-elastic Uf, — N or 

Ue — N scattering events which occur at p[^^ up to ^/rrii^^E^. Fig. 58 shows the 
estimated signal to background ratios expected in the low pt region. 

Because of the exceptionally low cross section, the target must be very mas- 
sive. The detector must therefore be capable of resolving the pt with much 
better resolution than the background spread. This favors the use of a fully 
active, high resolution tracking detector with sub-radiation length sampling in 
order to resolve the pt of the single electron before it is significantly broadened 
by shower development. A liquid Argon TPC, such as the one proposed for 
the ICANOE experiment [41] might be ideal for such a measurement. Another 
possibility would be a scintillating fiber/tungsten calorimeter. 
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Sensitivity vs y^u,, Ej,= 30 GeV 
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Figure 59: Statistical uncertainty in sin^ 9w for a luminosity of 10 
function of t/cut for a 30 GeV neutrino beam. Note that the ^~ produced beam 
is much less sensitive to sin^ 9w due to nearly exact cancellation in the sin^ 9w 
dependence of the two neutrino species in the beam. 



The largest experimental challenge for such a measurement is likely to be 
the normalization of the absolute neutrino flux. Despite the precise knowledge 
of muon decays, it would be extremely difficult to predict the precise neutrino 
flux at the 10""^ level merely from monitoring the parent muon beam. Instead, 
the signal processes will probably have to be normalized to the theoretically 
predictable processes of inverse muon decay, z/^e~ — > z^e/^~; and muon production 
through annihilation, Vee~ — > T'^^~ ^ both of which occur only in the v^JT^ beam. 
Normalization of the Ve beam may be possible through comparison to neutrino- 
nucleon scattering, vN — > l^N' , in the V^Ve and Vfjp^ beam. 

For a 20 ton detector at a 50 GeV muon storage ring, with 2 x 10^*^ jjT 
decays/year there will be approximately 1.5 x 10^° DIS charged current events 
and 8.5M z/^/z7e-clcctron scatters per year. This leads to an estimated sensitivity 
of (5sin2^gJ^*) ~ 0.0002 

Neutrino-nucleon scattering 

There were two dominant systematic uncertainties in present-day meaurements 
of the weak mixing angle in neutrino nucleon scattering, contamination in the 
beam and the kinematic suppression of scattering from strange quarks in the 

charged current channel. 

For an isoscalar target, the neutral current rates can be related to the charged 
current rates via [17]: 
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lsm''ew[r-rAR,] (97) 

y 

Rj,-AR, = {^-sm''ew + lsm^ew)[l-AR,] + 

^smHw[---AR,] (98) 
9 J- ^ 

where Ru/v is the ratio of neutral to charged current cross sections, r ~ 0.5 is 
the ratio of charged current anti-neutrino to neutrino cross sections, and ARs 
and Ai?c are small corrections for the kinematic suppression of s — c in charged 
current scattering where the neutral current process s — > s has no suppression. 
The charm corrections can be eliminated by a judicious subtraction of neutrino 
and anti-neutrino rates [97] but with a consequent reduction in statistical power. 

Present-day experiments [98, 88, 65] have had integrated fluxes of 10^^-10^^ 
neutrinos and have relied on dense nuclear targets. In such targets neutral 
current events are distinguished from charged current events by the presence or 
absence of a muon in the final state. In a dense calorimeter, electron neutrino 
charged current induced events look similar to neutral current events as the 
electron is lost in the hadronic shower. They are a significant background for 
precision measurements with conventional beams produced by pion and kaon 
decay and would be even more significant at a neutrino factory. 

The most precise measurement to date is from the NuTeV collaboration [65] 

of 

sin^ Ow = 0.2253 ± 0.0019(stoi) ± O.OOlO(sysi). (99) 

At a neutrino factory, the neutrino fiux will be several orders of magnitude 
higher but the beam will consist of approximately equal numbers of and Vf.. 
This makes a detector capable of distinguishing electron charged current events 
from neutral current events desirable and implies a low density detector such as 
those considered for the deep-inelastic scattering studies. 

We have considered several possible observables for a neutrino factory mea- 
surement and propose: 



a{v^, CC) - a(i/e, CC) 1 - gr 



or 



(7(i/^, CC) 

for the VpiPi beam, and 



a{V^,NC) + a{v^,NC) _ R" + g-^R"" 
a{ue,CC)-a{v^,CC) ~ 1 - g^'r 



(102) 
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or 



a{V^, NC) + (7(l/e, NC) + a(i/e), CC) 
a{v,,CC) 



(103) 



r r 



for the I^t'e beam, where R^/u is the ratio of neutral to charged current cross 
sections from Eq. (97). The observable i?^ requires electron identification while 
R'^ requires only muon identification. 

The variable g is the energy-weighted flux ratio between and or, equiv- 
alently, between 17^ and Ve' 

The flux ratio for neutrinos and anti-neutrinos g is: 



and is well determined by the muon decay kinematics. However, the relative 
detection efficiencies for muons and electrons must be known at the 2 x 10~^ 
level in order to determine sin^ 9\y to 10^'^ by the first method. In addition, the 
charm contributions are not cancelled in this observable and must be measured 
directly in the same experiment. 

For the Re measurement, which requires electron identification, an active 
target of 20 cm radius, 10 gr/cm^ thick consisting of either CCD's or silicon 
strip detectors (~ 140 300-/xm detectors) spaced over a meter and followed by 
the tracking, electromagnetic and hadron calorimetry and muon identification 
proposed above for structure function measurements would yield 15M muon and 
8M electron charged current deep- inelastic scattering events/ 10^° n~ decays and 
would yield a statistical precision of 0.0004 in sin^ 9w- The charm corrections 
partially cancel in this observable and would also be measured directly via the 
2M charm events/year produced in such a detector. 

The R measurement, which relies only on muon identification can be done 
with a much denser target, perhaps an iron/silicon sandwich calorimeter. Such a 
calorimeter 200 gr/cm^ thick would have a statistical sensitivity of Asin^ 9w ~ 
0.0001 per year at a 50 GeV machine. This method is quite similar to the method 
used in the NuTeV [65] measurement and would be dominated by systematic 
errors. 

4.10 Heavy Lepton Mixing 

A muon storage ring offers ample opportunities to search for new phenomena in 
yet unexplored physical regions. One such opportunity is the ability to search 



9 = 



mE^jE^JE,^ _ 6 

! ^{E^^)E,^dE,^ r 



(104) 
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for the existence of neutral heavy leptons. Several models describe heavy isospin 
singlets that interact and decay by mixing with their lighter neutrino counter- 
parts [108, 109]. The high intensity neutrino beam created by the muon storage 
ring provides an ideal setting to search for neutral heavy leptons with a mass 
below the muon mass, 105.6 MeV/c^. 

It is postulated that neutral heavy leptons (Lq) could be produced from muon 
decay when one of the neutrinos mixes with its heavy, isospin singlet partner. 
Neutral heavy leptons can be produced via one of two channels: 

fx- ^ Lo + Ve + e- (105) 

/x" ^ i/^ + Lo + e" (106) 
The branching ratio for each of these reactions is given by: 

BRiii ^ Lo/xe) = \Ui\\l - 8x1 + + 12a;^lna;^) (107) 

Here Xm = raL^/m^ and jC/jp is the mixing constant between the specific type of 
neutrino and the neutral heavy lepton: C/j = (Loli^j). Note that \U^\'^ and jC/gp 
need not be identical. 

Once produced, a neutral heavy lepton of such low mass can either decay via 
Lq — > vvv, Lq — > i/ee, or Lq — > 71/. The most viable mode for detection is the 
two-electron channel. For this particular decay mode, the Lq can decay either 
via charged current or charged and neutral current interactions. The branching 
ratio for this decay process has been previously calculated [110]. Since the decay 
width is proportional to f/|, the number of Lq's detectable is proportional to 
Uf ■ Uf in the limit where the distance from the source to the detector is short 
compared to the hfetime of the Lq. 

Using the above model, one can estimate the number of neutral heavy leptons 
produced at the muon storage ring which later decay within a given detector: 

Nl^ = N^* BR{ii Loue) * e * e"^/^^" * BR{Lo detectable) * (1 - e"'"'/^"") 

(108) 

Here N^, is the number of neutrinos produced from muon decay, BR{ii — > Loue) is 
the branching ratio of muons decaying into neutral heavy leptons versus ordinary 
muon decay, L is the distance from the beamline to the detector, 51 is the length 

of the detector, e is the combined detector and geometric efficiency, r is the Lq 
lifetime, and BR{Lq — > detectable) is the branching ratio for the neutral heavy 
lepton decaying via a detectable channel (presumably Lq uee). 

In estimating the sensitivity to Lq production at the muon storage ring, we 
make a few underlying assumptions. We assume that the storage ring utilizes 
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a pure, unpolarized muon beam with straight sections such that 25 percent of 
the muons will decay to neutrinos pointing towards the detector. We assume 
that the fiducial volume is 3 meters in diameter and 30 meters in length, (which 
is probably compatible with the need for empty space before a conventional 
detector) and that the detector has sufficient tracking resolution to detect the 
e^e~ vertex from the Lq decay. We assume for now that the background is 
negligible. These parameters correspond to the fiducial volume of the decay 
channel used for the Lq search at E815 (NuTeV) [111, 112]. 

The sensitivity of the detector has been calculated for a number of different 
muon energies and beam intensities. Fig. 60 shows limits on the Lq-v^j^ mixing 
as a function of Lq mass. One achieves the best limits from using relatively low 
energy/high intensity muon beams. This is a major improvement over previous 
neutral heavy lepton searches, where limits do not reach below 6.0 x 10~® in the 
low mass region [113, 109]. 

The single event sensitivity quoted here depends on having minimal back- 
ground levels in the signal region. Part of this can be achieved by kinematic 
cuts which discriminate against neutrino interactions in the detector material. 
However, it will probably be necessary to reduce the amount of material in the 
fiducial region compared to NuTeV. We estimate that even if the decay region is 
composed only of helium gas, the number of neutrino interactions will approach 
a few thousand. The ideal detector, therefore, would consist of a long vacuum 
or quasi-vacuum pipe with appropriate segmentation for tracking. The decay 
pipe could be used in conjunction with larger neutrino detectors adapted for the 
muon storage ring. 

The muon storage ring would to be an ideal location to continue the search 
for neutral heavy leptons. The high intensity neutrino beam allows for a neutral 
heavy lepton search to be sensitive to the 10 - 100 MeV/c^ mass range. In 
addition, such a neutral heavy lepton program is very compatible with a neutrino 
detector which uses the same neutrino beam. It is also clear, however, that 
a neutral heavy lepton search would receive the most benefit at lower muon 
energies, and thus would yield best results at the earlier stages of the muon 
storage ring program. 

4.11 Neutrino Magnetic Moments 

Although neutrino oscillation searches focus on the mass differences between 
neutrino eigenstates, neutrinos can possess other observable properties such as 
a magnetic moment. A measurement of the neutrino magnetic moment (NMM) 
would not only have great impact in the field of cosmology, particularly the 
development of stellar models, but would also help constrain several Standard 
Model extensions. An important experimental advantage is that a NMM search 
can run parasitically as the front-end of a typical long baseline detector. 

Despite their lack of charge, neutrinos can possess a non-zero magnetic mo- 
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Neutral Heavy Lepton Search at the Muon Storage 
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Figure 60: Limits on iC/y^Lp as a function of Lq mass for one year of running. The 
curves show sensitivities for 20 GeV and 50 GeV muon energies. Sensitivities 
assume no background events in signal region. 
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ment that can arise through loop diagrams. In the Standard Model, extended to 
include a right-handed neutrino or with left-handed neutrinos which have mass, 
the expected magnitude of the [109] neutrino magnetic moment is given by 

//.c.3x10-^Vb-(^), (109) 

1 eV 

where is the Bohr magneton. Although quite minuscule, several extensions 
to the Standard Model could dramatically increase /Xj, : supersymmetric models 
can produce 10"^^//^ to 10"^^//^ [115] and calculations that invoke large extra 
dimensions easily yield IO'^^/ib or larger [116]. 

Relative to the Standard Model expectation, the excluded values of NMM are 
not at all stringent, being seven to nine orders of magnitiidc larger. The current 
limits on neutrino magnetic moment from laboratory experiments are fi,y < 1.5 to 
1.8 X 10^^°/iB for electron neutrinos [117] [118] and < 7.4 x 10~^^hb for muon 
neutrinos [119]. Astrophysical limits are stronger: the slow rate of plasmon decay 
in horizontal branch stars [120] implies iJ,„ < 10~^^iib, while the neutrino energy 
loss rate from supernova 1987a [121] yields fii, < 10~^'^fiB- Note, however, that 
several assumptions are implicit in the astrophysics limits, including the core 
temperature of the stars; if stellar models omit important processes, these limits 
might be overestimates. The supernova limit applies only to Dirac neutrinos and 
not to the Majorana case. 

Existing search schemes possess a weakness that sharply limits their ulti- 
mate sensitivity: the formulae for the hypothesized effect are quadratic in /Xj^ 
but linear in terms of the experimenter-controlled parameters. In contrast, the 
following scheme is quadratic in terms of the product of the NMM and a mag- 
netic field strength, /^i/-B; hence a carefully designed and executed experiment 
could improve the limits from current experiments and possibly the limits from 
astrophysics calculations, or actually detect a NMM. 

The energy £J of a neutrino with a magnetic moment in a magnetic field B 
gains a new term ^^-B. Consider a B field along the rc-axis, and a neutrino with 
momentum and helicity along the 2;-axis at t = 0. The eigenstates of the neutrino 
are projections along the x-axis, and the state of the neutrino is expressed as: 

p-i{E+tj,^B)t -i{E-^^B)t 

IT) = ^^M + ^^H>. (110) 

As the neutrino propagates, the relative phase of the two components changes, 
corresponding to a rotation to a sterile state in the case of a Dirac neutrino or 
to an antineutrino in the Majorana case. At a far detector, the signal would 
be a deficit in the number of neutrinos detected or increase in the number of 
antineutrinos detected with the B field in place compared to the sample detected 
with no B field turned on. 

In this phase rotation scheme, the energy splitting occurs as the neutrino 
passes through a field gradient and experiences a force F — V(/x,^-B). To preserve 
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this energy difference, which drives the phase difference in the absence of the B- 
field, the field must be turned off instead of allowing the neutrino to experience 
the reverse gradient as it exits the field region. (The principle of changing the 
energy of neutral dipolar molecules with time-varying electric fields has been 
demonstrated in the laboratory [122].) 

To be successful, there are two basic requirements for the magnetic field: 

1) The magnetic field must oscillate such that the neutrino experiences only 
one sign of the gradient. This study assumes that the neutrino bunch 
length is small compared to the oscillation length. If this assumption is 
not true, the effects discussed here will be diluted but the basic conclusions 
will still apply. 

2) The magnetic field must be as strong as possible. 

We have explored the possibility [123] of using two existing technologies for 
the B field: resonant cavities and kicker magnets. In both cases the maximum 
magnetic field is too small to yield improved magnetic moment limits given 
realistic equipment. We are currently exploring configurations involving a series 
of pulsed current sheets. 

The formula for the number of events lost to sterile states may be expressed 
very simply as: 

Niost = N^sin^n^Bt) (111) 

where t is the neutrino flight time from entering the magnetic fleld to detection. 
We note that, in contrast to an oscillation disappearance signal, this effect is 
explicitly independent of the neutrino energy. Fig. 61 compares the number of 
events that vanish because of phase rotations to the expected statistical fluctu- 
ations in the number of CC events (N) observed in a 50 kton [124] far detector. 
We see that for a cummulative 10 T fleld gradient and 10^^ muon decays we 
expect a > lOcr signiflcance for a NMM of 10~^^hb- With a 3T gradient, the 
limit drops below two a. The sensitivity can be greatly increased by increasing 
the fleld strength and more weakly by increasing the number of events in the far 
detector. Because the detector distance deterimines both t and N in Eqn. Ill, 
the "signiflcance" in the flgure is linearly dependent on distance. 

To conclude, we have discussed a novel neutrino magnetic moment search 
technique that uses oscillating magnetic flelds at the source of a long baseline 
detector's neutrino beam. This is the only technique we know of that is quadratic 
in both III, and a controllable parameter, and thus has the potential for improved 
sensitivity as we improve our ability to create oscillating magnetic fleld gradients. 
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Figure 61: Significance for several scenarios. 
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4.12 Exotic processes 



Neutrino factories offer the possibility of searching for exotic processes resulting 
in production of e~, /x"*", or r-lepton of either charge. While these searches are 
interesting in their own right, they are also useful in ruling out exotic contri- 
butions to long-baseline neutrino oscillation signals. The neutrino beam from a 
muon storage ring would consist of a virtually pure combination of v^. and (or 
charge-conjugate). At very short baselines the Pg and will not have had time 
to oscillate into other flavors: For a 20 GeV muon storage ring with a 700 m 
straight section, and neutrino oscillations with Am^ = 3.5 x 10~^ eV^, the 
oscillation probability is ~ 5 x 10^^. 

One could distinguish between exotic processes and the beginning of a neu- 
trino oscillation by exploiting their differing dependence on energy and distance. 
Specifically, these exotic processes would probably have a flat or rising depen- 
dence on the neutrino energy E-^,. In contrast, a neutrino oscillation would have 
al/El dependence. Also, if the distance L of the experiment changes, the rate 
of exotic events would decrease with the flux as In contrast, the neutrino 

oscillation probability would increase as (for L small compared to the os- 
cillation period), and so the rate of oscillated events would be independent of 
L. 

Current understanding of muon interactions allows for exotic processes in two 
forms. Anomalous lepton production could occur if muons decay to neutrino 
flavors other than those in the usual decay /i ePe^'^" anomalous 
neutrinos then interact in the target. Alternatively, they could be produced if a 
Pg or Vjj^ interacts with the target via an exotic process. 

The only direct experimental limit on exotic jJL — > ev^i^y decays is BR{^ 
eP^i^e) < 1.3% [125]. Indirect hmits are also very weak. The contribution of non- 
V — A interactions to the muon decay rate has been limited to 8% [125]. Also, 
the total muon decay rate is one of the main measurements used to constrain 
electroweak parameters [125]. To first order, 

1 Gprrvl 

— = ^ 112 
Tf, 1927r3 ^ ' 

Assuming the standard model, Gp is determined to 1 part in 10^ from muon 
lifetime measurements. If there are exotic contributions to the muon lifetime, 
the measured value of would be shifted from the true value. Since 

mw oc Gp^^"^, (113) 

the 0.1% uncertainty on mw corresponds to a 0.4% shift in the muon lifetime. 
Finally, the CKM matrix element V^d is determined from the rate of nuclear /3- 
decays relative to the muon hfetime. The assumption of unitarity on the CKM 
matrix gives us the following constraint on the first row: 

iKdr+iK.i'+iKdr = i. (114) 
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The experimental determination is[125]: 

\Vud? + + \Vud? = 0.991 ± 0.005. (115) 

The uncertainty on this constraint corresponds to a 0.5% shift in the muon hfe- 
time. Additional contributions to the muon decay rate would lead to a downward 
shift in the determined value of iKdp from the true value. We conclude that 
exotic decay modes of the muon with branching ratios totaling 0.5% are possible 
without contradicting current measurements or tests of the standard model. 

As a concrete example of such an exotic process we consider R-parity- violating 
super symmetric models. These models lead to Icpton-number-violating vertices 
with couplings A, and muon decay processes such as — > e^v^^r as shown in 
Fig. 1. The matrix element for these decays turns out to have the same form 
as for the standard W-exchange. The current constraints on the couplings A are 
reviewed in Ref. [126]. These constraints allow a branching ratio of 0.4% for the 
process in Fig. 62. 

Similar processes are allowed for anomalous lepton production as shown for 
example in Fig. 63. Estimates for allowed rates are in progress [127]. These 
diagrams involve the A' couplings. Currently, the best hmit on one of these 
couplings, A231, is from i/^ deep- inelastic scattering, so existing neutrino data is 
already providing constraints! The search for these types of effects at the muon 
storage ring could be input into a decision on whether to build a muon-proton 
collider where they could be studied in more detail. 

As a start on estimating the capabilities of an experiment at the neutrino 
source, we consider the detector concept illustrated in Fig. 64. This concept 
consists of a repeating sequence of 1.5 mm-thick Tungsten sheets with Silicon 
tracking, separated by 4 mm. Tungsten, being dense, provides a high target 
mass while being thin enough for a produced r to have a high probability of 
hitting the Silicon. The impact parameter of the r decay products is typically 
90 microns with a broad distribution, so we would like a hit resolution of 5 
microns or better. Although there is a lot of multiple scattering in the tungsten, 
the short extrapolation distance provides for a good impact parameter resolution 
on the r decay products. This configuration has been optimized for a 50 GeV 
muon beam. For lower energy beams, the planes should be spaced more closely, 
and the Tungsten thickness perhaps reduced. Studies of detectors with passive 
target mass and tracking with emulsion sheets [128] suggest that we can expect 
T reconstruction efficiencies as high as 30%. 

We would propose placing such a detector in a magnetic field, and measuring 
the momentum of muons and hadrons should be straightforward. However, each 
Tungsten sheet is 0.4 radiation lengths thick, and while we should obtain good 
energy resolution for electromagnetic showers, it will not be feasible to measure 
the charge of an electron before it showers. 

In summary, with a total mass of 6 tons of Tungsten, 200 m^ of Silicon 
tracking, located close to a muon storage ring withy 5 x lO^'^ muon decays at 50 
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Figure 62: Example of exotic muon decay in R-parity-violating SUSY. 
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Figure 63: Example of an exotic neutrino interaction in R-parity-violating SUSY. 



GeV, we expect a total of 35 billion neutrino interactions, 4 orders of magnitude 
above present neutrino interaction samples. Thus, there is much potential for 
detecting very rare exotic processes if we can adequately reduce backgrounds. 
Detailed simulations and studies of possible Silicon tracking technologies are 
needed to quantify this. 



4.13 Summary 

We have investigated possible conventional neutrino physics studies done at a 
detector located near a muon storage ring. We emphasized novel methods rather 
than extensions of existing experiments with additional statistics. 

For a reference machine with 50 GeV stored muons and 10^° muon decays 
per year we find that it is possible to: 

• Measure individual parton distributions within the proton for all light 
quarks and ant i- quarks. 

• Determine the effects of a nuclear environment on individual quark species. 
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• Measure the spin contributions of individual quark species, including strange 
quarks and do precision studies of the QCD evolution of spin effects with- 
out a need for data from beta decay measurements. 

• Measure charm production with raw event rates of up to 100 million charm 
events/year with ~ 2M double tagged events. 

• Measure the Weinberg angle in both hadronic and purely leptonic modes 
with a precision of 0.0001 to 0.0002. 

• Search for the production and decay of neutral heavy leptons with mixing 
angle sensitivity two orders of magnitude better than present limits in the 
30-80 MeV region. 

• Search for a neutrino magnetic moment which may be much larger than 
the Standard Model prediction in some supersymmetric theories. 

• Search for anomalous tau production predicted by some R-parity violating 
supersymmetric models. 

We note that the event rates at a near detector increase linearly with neutrino 
energy. In addition, the acceptance of small detectors is better for the narrower 
beam produced by higher energy machines. Almost all of the above measure- 
ments, with the exception of the neutral heavy lepton search, lose sensitivity if 
the beam energy is less than 50 GeV and gain if it is greater. 

If the storage ring beam energy is lowered to 20 GeV, the statistical power 
of almost all of the measurements considered here would drop a factor of 2.5 
or more. The number of deep-inelastic scattering events with high enough 
for perturbative QCD to be meaningful drops even further and the minimum 
X rises to 0.05. Measurements involving charm or tau production in the final 
state would be have lower statistics due to threshold effects, as would the inverse 
muon decay normalization for v — e scattering, which has a threshold of ~ 11 
GeV. 
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However, it should be remembered that a 50 GcV neutrino factory will pro- 
duce neutrino fluxes of order 10"^ higher than existing neutrino beams. At 20 
GeV the improvement for most physics processes is still greater than a factor of 
a thousand. 

5 Summary and Recommendations 

The main goal of the physics study presented in this report has been to answer 
the question: Is the physics program at a neutrino factory compelling ? The 

answer is a resounding yes, provided there are 10^^ or more muon decays per 
year in the beam forming straight section and the muon beam energy is ~ 20 GeV 
or greater. Based on our study, we believe that a neutrino factory in the next 
decade would be the right tool at the right time. 

The neutrino oscillation physics program 

The recent impressive atmospheric neutrino results from the Super-Kamiokande 
experiment have gone a long way towards establishing the existence of neutrino 
oscillations. This is arguably the most dramatic recent development in particle 
physics. Up to the present era, neutrino oscillation experiments at accelerators 
were searches for a phenomenon that might or might not be within experimental 
reach. The situation now is quite different. The atmospheric neutrino deficit 
defines the Sm^ and oscillation amplitude to which future long baseline oscillation 
experiments must be sensitive to, namely 5w? = 0(10^^) eV^ and sin^ 29 = 0(1). 
Experiments that achieve these sensitivities are guaranteed an excellent physics 
program that addresses fundamental physics questions. Furthermore, should all 
of the experimental indications for oscillations (LSND, atmospheric, and solar) 
be confirmed, we may be seeing evidence for the existence of sterile neutrinos. 
This would be a very exciting discovery which would raise many new questions 
requiring new experimental input. 

A neutrino factory would be a unique facility for oscillation physics, pro- 
viding beams containing high energy electron neutrinos (antineutrinos) as well 
as muon antineutrinos (neutrinos). These beams could be exploited to provide 
answers to the questions that we are likely to be asking after the next generation 
of accelerator based experiments. 

The oscillation physics that could be pursued at a neutrino factory is com- 
pelling. Experiments at a neutrino factory would be able to simultaneously mea- 
sure, or put stringent limits on, all of the appearance modes — > Ut-, — > u^, 
and Ufj^ — > Ut-. Comparing the sum of the appearance modes with the disappear- 
ance measurements would provide a unique basic check of candidate oscillation 
scenarios that cannot be made with a conventional neutrino beam. In addition, 
for all of the specific oscillation scenarios we have studied, the component 
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in the beam can be exploited to enable crucial issues to be addressed. These 
include: 

(i) A precise determination of (or stringent limits on) all of the leading oscilla- 

tion parameters, which in a three-flavor mixing scenario would be sin^ 2^13, 
sin^ 26*23, and Sml2- 

(ii) A determination of the pattern of neutrino masses. 

(iii) A quantitative test of the MSW effect. 

(iv) Stringent limits on, or the observation of, CP violation in the lepton sector. 

To be more quantitative in assessing the beam energy, intensity, and baseline 
required to accomplish a given set of physics goals it is necessary to consider two 
very different experimental possibihties: (a) the LSND oscillation results are not 
confirmed, or (b) the LSND results are confirmed. 

(a) LSND not confirmed. A 20 GeV neutrino factory providing 10^^ muon 
decays per year is a good candidate "entry-level" facility which would 
enable either (i) the first observation of — > t'^, oscillations, the first 
direct measurement of matter effects, and a determination of the sign of 
Sml2 and hence the pattern of neutrino masses, or (ii) a very stringent limit 
on sin^ 2^13 and a first comparison of the sum of all appearance modes with 
the disappearance measurements. The optimum baselines for this entry- 
level physics program appears to be of the order of 3000 km or greater, 
for which matter effects are substantial. Longer baselines also favor the 
precise determination of sin^ 2^13. A 20 GeV neutrino factory providing 
10^° muon decays per year is a good candidate "upgraded" neutrino factory 
(or alternatively a higher energy facihty providing a few xlO^^ decays per 
year). This would enable the first observation of, or meaningful limits on, 
z/g oscillations, and precision measurements of the leading oscillation 
parameters. In the more distant future, a candidate for a second (third 
?) generation neutrino factory might be a facility that provides 0(10^^) 
decays per year and enables the measurement of, or stringent limits on, 
CP violation in the lepton sector. 

(b) LSND confirmed. Less extensive studies have been made for the class of sce- 

narios that become of interest if the LSND oscillation results are confirmed. 
However, in the scenarios we have looked at we find that the z/g Vr rate 
is sensitive to the oscillation parameters and can be substantial. With a 
large leading Sni^ scale medium baselines (for example a few xlO km) are 
of interest, and the neutrino factory intensity required to effectively ex- 
ploit the z/g beam component might be quite modest (< lO-*^^ decays per 
year). If sterile neutrinos play a role in neutrino oscillations, we will have 
an exciting window on physics beyond the SM, and we anticipate that a 
neutrino factory would enable crucial measurements to be made exploiting 
the electron neutrino beam component. 
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The non-oscillation physics program 

A neutrino factory could provide beams that are a factor of 10^ — 10^ more in- 
tense than conventional neutrino beams. This would have an enormous impact 
on the detector technology that could be used for non-oscillation neutrino exper- 
iments. For example, the use of silicon pixel targets and hydrogen or deuterium 
polarized targets would become feasible. Hence, a neutrino factory would offer 
experimental opportunities that do not exist with lower intensity conventional 
beams. 

We have looked at a few explicit examples of interesting experiments that 
might be pursued at a neutrino factory: 

• Precise measurements of the detailed structure of the nucleon for each 
parton flavor, including the changes that occur in a nuclear environment. 

• A first measurement of the nucleon spin structure with neutrinos. 

• Charm physics with several million tagged particles. Note that charm 
production becomes significant for storage ring energies above 20 GeV. 

• Precise measurements of standard model parameters - cts, the weak mixing 
angle, and the CKM matrix elements. 

• Searches for exotic phenomena such as neutrino magnetic moments, anomolous 
couplings to the tau and additional neutral leptons. 

The non-oscillation measurements benefit from higher beam energies since 
event rates and the kinematic reach scale with energy, and perturbativc calcula- 
tions become more reliable in the kinematic regions accessed by higher energies. 

Recommendations 

The physics program we have explored for a neutrino factory is compelling. We 
recommend a sustained effort to study both the physics opportunities and the 
machine realities. 

(i) We encourage support for the R&D needed to learn whether a neutrino 

factory can be a real option in the next decade. 

(ii) We propose further studies of detector technologies optimized for a neutrino 

factory, including both novel low mass detectors for near experiments and 
very high mass detectors for long baselines. For long baseline experiments 
detectors should have masses of a few times 10 kt or more that are able 
to detect and measure wrong-sign muons, and detectors of a few kt or 
more able to observe tau-lepton appearance with high efficiency. It is also 
desirable to identify electrons, and if possible measure the sign of their 
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charge. Both the detector technologies themselves and the civil engineering 
issues associated with the construction of such massive detectors need to 
be addressed. 

(iii) We recommend continued studies to better compare the physics potential 
of upgraded conventional neutrino beams with the corresponding potential 
at a neutrino factory, and also studies to better understand the benefits of 
muon polarization. 

(iv) The present study concentrated on the muon storage ring as a neutrino 
source and did not cover the additional physics programs which would use 
the proton driver and the high intensity muon beams. We recommend a 
further study directed at these other facets of physics at a muon storage 
ring facility. 
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